Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 478,502 C→A 85.0% D342E (GAC→GAA USA300HOU_0457 → NADH dehydrogenase (ubiquinone)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730478,5020CA85.0% 126.1 / 14.8 68D342E (GAC→GAAUSA300HOU_0457NADH dehydrogenase (ubiquinone)
Reads supporting (aligned to +/- strand):  ref base C (6/5);  new base A (26/31);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

AGTTCAATGTGCATTGGGTGTATATTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGACGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTACTTGTATCTTG  >  CP000730/478364‑478633
                                                                                                                                          |                                                                                                                                   
aGTTCAATGTGCATTGGGTGTATATTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTCAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATCCCAAAACAAGCATCTGCTAAAGAAGc                                                                                                                                   >  2:460820/1‑141 (MQ=255)
 gTTCAATGTGCATTGGGTGTATATTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCt                                                                                                                                  >  1:246446/1‑141 (MQ=255)
  ttCAATGTGCATTGGGTGTATATTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCtt                                                                                                                                 <  1:280789/141‑1 (MQ=255)
  ttCAATGAGCATTGGGTGTATATTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCAAAAGGAGCtt                                                                                                                                 >  1:73663/1‑141 (MQ=255)
          gCATTGGGTGTATATTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGACgaaactggcagg                                                                                                                         >  1:314565/1‑130 (MQ=255)
          gCATTGGGTGTATATTCAGCAGCGATTATTCATGTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTgg                                                                                                                         >  2:571138/1‑141 (MQ=255)
           cATTGGGTGTATATTGAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATCTTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGa                                                                                                                        >  1:452623/1‑141 (MQ=255)
             ttGGGTGTATATTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGAtt                                                                                                                      <  2:19307/141‑1 (MQ=255)
             ttGGGTGTATATTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAACCATCTGCTAACGAAGCTTATGGCTGGAtt                                                                                                                      >  2:174722/1‑141 (MQ=255)
               gggTGTATATTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGt                                                                                                                    >  1:87244/1‑141 (MQ=255)
                      tatTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTTTTTTTTTTTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATgggg                                                                                                             >  2:535061/1‑140 (MQ=255)
                      tatTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGa                                                                                                             >  1:528353/1‑141 (MQ=255)
                        ttCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACg                                                                                                           <  1:487379/141‑1 (MQ=255)
                          cagcagCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACgtg                                                                                                         <  2:582209/141‑1 (MQ=255)
                                cGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGACGCTTATGGCTGGATTGTCATTGGACGTGTATGAg                                                                                                   >  2:86128/1‑141 (MQ=255)
                                 gATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGc                                                                                                  <  1:204105/141‑1 (MQ=255)
                                  attattCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGACGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCt                                                                                                 <  1:121432/141‑1 (MQ=255)
                                   ttattCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGACGCGTATGGCTGGATTGTCGTGGGACGGGTATTAGCta                                                                                                >  2:307422/1‑141 (MQ=255)
                                        taTTTAATATTGCACTGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATc                                                                                           <  2:252965/140‑1 (MQ=255)
                                          tttAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGt                                                                                         >  2:344032/1‑141 (MQ=255)
                                               tatTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATTTTTGAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCAt                                                                                    <  1:571138/141‑1 (MQ=255)
                                                tttGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAACCAAGCATCTGCCAAAGACGTTTATGGCTGGATTGTCATGGGACGTGTATCAGCTATTACCGTGGCAtt                                                                                   <  2:308332/140‑1 (MQ=255)
                                                atTGCACGGTATTTTTAAAGCAACATTATTTTTACATTCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGACGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCAtt                                                                                   <  2:308309/141‑1 (MQ=255)
                                                atTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTCAAGAAGCTTATGGCTGGATTGTCATGGGCCGCGTATTAGCTATTATCGTGGCAtt                                                                                   >  2:181196/1‑141 (MQ=255)
                                                atTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCAtt                                                                                   <  2:256229/141‑1 (MQ=255)
                                                   gCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGACGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCt                                                                                >  1:437870/1‑141 (MQ=255)
                                                   gCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCt                                                                                <  1:191914/141‑1 (MQ=255)
                                                     ccGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGACGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTAt                                                                              <  2:11377/140‑1 (MQ=255)
                                                     aCGGTTTTTCTAACATAACATTATTTTTACCATCGGGAACGAATGTGGCGCCATTTCATGTTCCAAAACAAGCATCTGCTAATGAAGCTTATGGCGGGGTAGTCGTTGTACGTGTATTAGCTATTTTCGGGGCAGTTCTAt                                                                              <  2:215065/141‑1 (MQ=255)
                                                                 aaGCAACATTATTTTTAGAATCAGGTTCTATCGTGAAGCGATTCAAGATTCAAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGagta                                                                  <  1:218484/141‑1 (MQ=255)
                                                                          tattTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGaa                                                         >  1:207025/1‑141 (MQ=255)
                                                                            tttttACAATCAGGTTCTATCGTGAAGCGATTCAATATTTCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGt                                                       >  1:574726/1‑141 (MQ=255)
                                                                            tttttACAAGCAGGCTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGt                                                       <  2:245022/141‑1 (MQ=255)
                                                                            tttctACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGt                                                       <  2:244977/141‑1 (MQ=255)
                                                                             ttttACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTg                                                      <  2:87244/141‑1 (MQ=255)
                                                                              tttACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGc                                                     <  2:47858/141‑1 (MQ=255)
                                                                                tACAATCAGGTACTATCGGGAAGCGATTCAATATTCCAAAACAAGGAGCTGCTAAAGAAGCTTATGGCGGGATTGTCATGGGACGAGTATTAGATATTATCGTGGCATTTCTATTTTGGCTGAGTATTgagagaggggaat                                                   >  2:68217/1‑134 (MQ=255)
                                                                                   caTCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATg                                                <  1:210985/140‑1 (MQ=255)
                                                                                      cAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTACCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAg                                             <  1:406892/141‑1 (MQ=255)
                                                                                      cAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCAGGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAg                                             <  1:153111/141‑1 (MQ=255)
                                                                                        ggTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAgtg                                           >  1:474158/1‑141 (MQ=255)
                                                                                        ggTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAgtg                                           >  1:578388/1‑141 (MQ=255)
                                                                                          ttCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCAGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGtt                                         <  2:402101/141‑1 (MQ=255)
                                                                                           tCTATCGTGAAGCGATTCAATATCCCAAAACAAGCATCTGCTAAAGAAGCTCATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCacttaacacttca                                        <  1:460820/141‑13 (MQ=255)
                                                                                           tCTATCGTGAAACGATTCAATATTCCAAAACAAGAATCTGCTAAGGCAGCATATGGCTGGTTTGTAATGGGACGTGTATTAGCTATTCTCGTGGCAATTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGCGTTa                                        >  2:153502/1‑141 (MQ=255)
                                                                                            cTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTaa                                       >  1:17592/1‑141 (MQ=255)
                                                                                            cTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGGGCATATGAAGTGTTaa                                       >  1:261482/1‑141 (MQ=255)
                                                                                             tATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGGGTATTAGCTATTATCGTGGCGTTTCTATTTTGGCGGAGTAGTGACAGAAGTGCATATGACGTGTTAAg                                      >  1:153749/1‑141 (MQ=255)
                                                                                               tCGGGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGGCAGAAGTGCATATGAAGTGTTAAGTg                                    <  2:380032/141‑1 (MQ=255)
                                                                                                       cGATTCAATATTCCAGAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGTCGTGTATTAGCTATTATCGTGGCAGTTATATTGTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCAtt                            >  1:451514/1‑141 (MQ=255)
                                                                                                        gATTCAATATTCCAAAACCCGCGCCACCGAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTc                           >  2:588117/1‑141 (MQ=255)
                                                                                                         gtttagtcttCCAGACCAAGCAGCGGCTAAAGGAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTAGTATCGTGGCATTGCTATTTTGGATTAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCt                          <  2:385398/133‑1 (MQ=255)
                                                                                                          ttCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGCCAGAAGTGCATATGAAGTGTTAAGTGGACTCATTCTa                         >  1:33012/1‑141 (MQ=255)
                                                                                                              atactgcAAAACAAGGATCTGCGAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGGATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGAt                     <  2:265573/135‑1 (MQ=255)
                                                                                                                    cAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCAt               >  1:239402/1‑141 (MQ=255)
                                                                                                                     aaaaCAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATTGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGTATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCAtt              <  2:593748/141‑1 (MQ=255)
                                                                                                                       aaCAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACTAGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTAc            <  2:161670/141‑1 (MQ=255)
                                                                                                                       aaCAAGCATCTGCTAAAGAAGCTTATGGCGGGATTGTCATGGGACGAGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTAc            <  2:364051/141‑1 (MQ=255)
                                                                                                                        aCAAGCATCTGCTAAAGACGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTTTAGGATGGTCATTACt           <  1:303837/141‑1 (MQ=255)
                                                                                                                        aCAAGCAGCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTACt           <  1:358369/141‑1 (MQ=255)
                                                                                                                          aaGCGTCTGCTAAAGCAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTACTTg         <  1:505882/141‑1 (MQ=255)
                                                                                                                          aaGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATTAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTACTTg         <  1:51079/141‑1 (MQ=255)
                                                                                                                          aaGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTACTTg         >  1:257842/1‑141 (MQ=255)
                                                                                                                            gCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTAGTAGCTATTATCGTGGCATTTCTATTTTGGATGGGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTACTTGTa       <  2:427127/141‑1 (MQ=255)
                                                                                                                             cATCTGCTAAAGACGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTACgtgtgt      >  1:151216/1‑139 (MQ=255)
                                                                                                                              aTCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTACTTGTATc     <  2:59593/141‑1 (MQ=255)
                                                                                                                               tCTGCTAAAGACGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTACTTGTATCt    >  2:105929/1‑141 (MQ=255)
                                                                                                                                 tGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTACTTGTATCTTg  >  2:242792/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                   
AGTTCAATGTGCATTGGGTGTATATTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGACGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTACTTGTATCTTG  >  CP000730/478364‑478633

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: