Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,293,295 C→T 5.7% V149I (GTT→ATT)  napG ← ferredoxin‑type protein essential for electron transfer from ubiquinol to periplasmic nitrate reductase (NapAB)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,293,2950CT 5.7% 1191.5 / 38.8 332V149I (GTT→ATT) napGferredoxin‑type protein essential for electron transfer from ubiquinol to periplasmic nitrate reductase (NapAB)
Reads supporting (aligned to +/- strand):  ref base C (146/167);  new base T (13/6);  total (159/173)
Fisher's exact test for biased strand distribution p-value = 9.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 360 total aligned reads displayed.

CAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAATGTCCTCACACATTTCG  >  CP009273/2293146‑2293443
                                                                                                                                                     |                                                                                                                                                    
cAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATaacaaac                                                                                                                                                      >  2:2124371/1‑150 (MQ=255)
 aGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATaacaaaca                                                                                                                                                     <  2:1207381/150‑1 (MQ=255)
 aGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATaacaaaca                                                                                                                                                     >  2:105451/1‑150 (MQ=255)
   cACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATc                                                                                                                                                   <  2:2614558/150‑1 (MQ=255)
       cATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATcgcag                                                                                                                                               <  1:4629368/150‑1 (MQ=255)
       cATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATcgcag                                                                                                                                               >  1:778713/1‑150 (MQ=255)
       cATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATcgcag                                                                                                                                               >  1:4435909/1‑150 (MQ=255)
        aTACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATcgcagc                                                                                                                                              <  2:1108483/150‑1 (MQ=255)
         tccctttgCGCCCTGACCGCACCCGGTACAGGCGTCGCTGTGACCCGGCGGCCGACAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCCGGGGGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATcgcagcg                                                                                                                                             <  2:2828922/146‑1 (MQ=255)
           cTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATcgcagcgca                                                                                                                                           >  2:4135647/1‑150 (MQ=255)
            tttttCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATcgcagcgcag                                                                                                                                          >  1:253952/1‑150 (MQ=255)
            tttttCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATcgcagcgcag                                                                                                                                          >  1:254234/1‑150 (MQ=255)
              tttCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGcc                                                                                                                                        <  1:5593044/150‑1 (MQ=255)
              tttCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGcc                                                                                                                                        <  2:2197973/150‑1 (MQ=255)
               ttCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACTTCGCAGCGCAGccc                                                                                                                                       >  2:4456526/1‑150 (MQ=255)
                tCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGcccc                                                                                                                                      <  1:887468/150‑1 (MQ=255)
                tCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGcccc                                                                                                                                      >  1:3378925/1‑150 (MQ=255)
                 cGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCTATTTTCGGGCCTGCGCGTGAACATGCAGCGCAGCGCCGcccgg                                                                                                                                     >  2:3501074/1‑148 (MQ=255)
                  gCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTg                                                                                                                                    >  1:2396137/1‑150 (MQ=255)
                  gCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTg                                                                                                                                    >  1:4014306/1‑150 (MQ=255)
                  gCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTg                                                                                                                                    <  1:2785221/150‑1 (MQ=255)
                   cACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGa                                                                                                                                   >  1:1318079/1‑150 (MQ=255)
                    aCTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGaa                                                                                                                                  >  2:2307691/1‑150 (MQ=255)
                    aCTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGaa                                                                                                                                  >  2:330675/1‑150 (MQ=255)
                       tACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAATATCGCAGCGCAGCCCCTGAAAGt                                                                                                                               >  1:3462220/1‑150 (MQ=255)
                        aCCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGtt                                                                                                                              <  1:4502652/150‑1 (MQ=255)
                         ccGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTg                                                                                                                             <  1:496237/150‑1 (MQ=255)
                         ccGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTg                                                                                                                             >  2:5286009/1‑150 (MQ=255)
                            cAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAATATCGCAGCGCAGCCCCTGAAAGTTgaga                                                                                                                          >  1:5320443/1‑150 (MQ=255)
                                cGGGACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGt                                                                                                                      <  1:2711561/150‑1 (MQ=255)
                                     cAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTcc                                                                                                                 <  2:469171/150‑1 (MQ=255)
                                      aGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCt                                                                                                                <  2:1781519/150‑1 (MQ=255)
                                      aGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCt                                                                                                                <  1:2218220/150‑1 (MQ=255)
                                            cGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTcca                                                                                                          >  1:1427642/1‑150 (MQ=255)
                                                    aCCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTAcc                                                                                                  <  1:405774/150‑1 (MQ=255)
                                                    aCCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTAcc                                                                                                  <  1:514666/150‑1 (MQ=255)
                                                     ccGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCg                                                                                                 <  2:4069903/150‑1 (MQ=255)
                                                              gAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCa                                                                                        >  1:434465/1‑150 (MQ=255)
                                                              gAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATAAATTTTCGGGCATTCTCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCATGGTCCACCAGTACCTCCAGCCCCa                                                                                        >  1:1191061/1‑150 (MQ=255)
                                                                   cGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAATATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCcgc                                                                                   >  1:4492995/1‑150 (MQ=255)
                                                                      gCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCcgcgcg                                                                                <  1:1767377/150‑1 (MQ=255)
                                                                        gTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGcgtc                                                                              <  2:3363732/150‑1 (MQ=255)
                                                                         tGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGcgtcg                                                                             <  2:171723/150‑1 (MQ=255)
                                                                          gCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGcgtcgt                                                                            <  2:4946177/150‑1 (MQ=255)
                                                                           cTTACCGGTACGCGTGTTGCGCTCCAGGTCCTGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCgtcgtc                                                                           <  1:2960308/150‑1 (MQ=255)
                                                                           cTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCgtcgtc                                                                           >  1:623327/1‑150 (MQ=255)
                                                                           cTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCgtcgtc                                                                           <  2:4900347/150‑1 (MQ=255)
                                                                            ttACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGtcgtcg                                                                          >  2:2425018/1‑150 (MQ=255)
                                                                            ttACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCAACCAGTACCGCCAGCCCCATCCGCGCGtcgtcg                                                                          >  2:3225401/1‑150 (MQ=255)
                                                                             tACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGtcga                                                                         >  1:3656066/1‑150 (MQ=255)
                                                                              aCCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTGTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGtcgat                                                                        >  1:3037878/1‑150 (MQ=255)
                                                                               ccGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGtcgatc                                                                       <  1:473371/150‑1 (MQ=255)
                                                                               ccGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGtcgatc                                                                       <  2:4354308/150‑1 (MQ=255)
                                                                               ccGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGtcgatc                                                                       >  1:4137204/1‑150 (MQ=255)
                                                                                    aCGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGTCAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTc                                                                  >  2:2564102/1‑150 (MQ=255)
                                                                                    aCGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTc                                                                  <  2:1281723/150‑1 (MQ=255)
                                                                                     cgcgTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCa                                                                 >  2:4249656/1‑150 (MQ=255)
                                                                                      gcgTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCaa                                                                <  2:4869636/150‑1 (MQ=255)
                                                                                       cgTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAAt                                                               <  2:253952/150‑1 (MQ=255)
                                                                                         tgtTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAAtct                                                             >  1:3288415/1‑150 (MQ=255)
                                                                                          gtTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAAtctc                                                            <  2:3669459/150‑1 (MQ=255)
                                                                                           ttGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGATTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCa                                                           >  1:3115154/1‑150 (MQ=255)
                                                                                            tGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCAc                                                          <  1:5262892/150‑1 (MQ=255)
                                                                                                      ctccagGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAgcgc                                                >  1:1971496/1‑150 (MQ=255)
                                                                                                      ctccagGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAgcgc                                                >  2:759753/1‑150 (MQ=255)
                                                                                                       tccagGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCa                                               <  1:4473078/150‑1 (MQ=255)
                                                                                                        ccagGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCAc                                              <  2:1115980/150‑1 (MQ=255)
                                                                                                         cagGGTGATGGCCTCATCAATTTTCGGGCCTTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCAcc                                             >  2:5072532/1‑150 (MQ=255)
                                                                                                         cagGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCAcc                                             <  2:4845510/150‑1 (MQ=255)
                                                                                                          agGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCg                                            >  2:653288/1‑150 (MQ=255)
                                                                                                              tGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTg                                        >  2:4652701/1‑150 (MQ=255)
                                                                                                              tGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTg                                        <  2:852118/150‑1 (MQ=255)
                                                                                                               gATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCTCACCGcattg                                       >  2:3978258/1‑146 (MQ=255)
                                                                                                               gATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTgg                                       <  2:5452733/150‑1 (MQ=255)
                                                                                                                   gCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGcac                                   <  1:3886366/150‑1 (MQ=255)
                                                                                                                   gCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGcac                                   <  2:1616613/150‑1 (MQ=255)
                                                                                                                    ccTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCGCGATCCAGCGCACCGCTTGGGcaca                                  <  2:1312882/150‑1 (MQ=255)
                                                                                                                      tcatcaATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACt                                >  1:3427659/1‑150 (MQ=255)
                                                                                                                       catcaATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACtt                               >  2:2805587/1‑150 (MQ=255)
                                                                                                                        atcaATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACttt                              <  2:3663736/150‑1 (MQ=255)
                                                                                                                        atcaATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACttt                              <  2:2443841/150‑1 (MQ=255)
                                                                                                                        atcaATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACttt                              >  2:4789691/1‑150 (MQ=255)
                                                                                                                          caATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTgg                            <  1:4194688/150‑1 (MQ=255)
                                                                                                                            aTTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGgcg                          >  1:4218991/1‑150 (MQ=255)
                                                                                                                            aTTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGgcg                          >  1:439066/1‑150 (MQ=255)
                                                                                                                                tCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACg                      >  2:2305382/1‑150 (MQ=255)
                                                                                                                                    gCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAAt                  >  2:2237667/1‑150 (MQ=255)
                                                                                                                                        tCGCGATAACAAATATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAATGTcc              >  1:367926/1‑150 (MQ=255)
                                                                                                                                          gcgATAACAATCATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAATGTCCTc            >  2:3925635/1‑150 (MQ=255)
                                                                                                                                            gATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAATGTCCTcac          >  1:1640172/1‑150 (MQ=255)
                                                                                                                                             aTAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAATGTCCTcaca         <  2:3143934/150‑1 (MQ=255)
                                                                                                                                              tAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCTCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCAATCGCACCGCTTGGGCACACTTTGGCGAACGGAATGTCCTcacac        >  2:3547257/1‑150 (MQ=255)
                                                                                                                                              tAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAATGTCCTcacac        >  2:124422/1‑150 (MQ=255)
                                                                                                                                              tAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAATGTCCTcacac        <  1:4789747/150‑1 (MQ=255)
                                                                                                                                                acaaacaTCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAATGTCCTCACACAt      <  1:105451/150‑1 (MQ=255)
                                                                                                                                                 caaacaTCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAATGTCCTCACACAtt     >  2:3090984/1‑150 (MQ=255)
                                                                                                                                                  aaacaTCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAATGTCCTCACACAttt    <  1:1702880/150‑1 (MQ=255)
                                                                                                                                                  aaacaTCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAATGTCCTCACACAttt    >  1:288687/1‑150 (MQ=255)
                                                                                                                                                   aacaTCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAATGTCCTCACACATTTc   >  1:5150377/1‑150 (MQ=255)
                                                                                                                                                    acaTCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAATGTCCTCACACATTTCg  >  1:572932/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGCACGCATACTTTTTCGCACTTACCGCAACCGGTACAGGCGTCGCTGTGAACCGTCGGCAGAAAGCGGGCGTGCTTACCGGTACGCGTGTTGCGCTCCAGCTCCAGGGTGATGGCCTCATCAATTTTCGGGCATTCGCGATAACAAACATCGCAGCGCAGCCCCTGAAAGTTGAGACAGTTTTCCTGGTCCACCAGTACCGCCAGCCCCATCCGCGCGTCGTCGATCGATTCAATCTCACGATCCAGCGCACCGCTTGGGCACACTTTGGCGCACGGAATGTCCTCACACATTTCG  >  CP009273/2293146‑2293443

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