Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,537,806 A→T 100% L171* (TTA→TAA)  narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,537,8060AT98.5% 1274.1 / 4.6 389L171* (TTA→TAA) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base T (204/179);  total (207/182)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 390 total aligned reads displayed.

GTAGCGGTACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTT  >  CP009273/1537657‑1537952
                                                                                                                                                     |                                                                                                                                                  
  aGCGGTACCCACATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAACATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTAt                                                                                                                                                  <  1:4101527/150‑1 (MQ=255)
     ggTACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCACATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCg                                                                                                                                               >  2:4426950/1‑150 (MQ=255)
       tACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGcc                                                                                                                                             >  2:988502/1‑150 (MQ=255)
       tACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGcc                                                                                                                                             >  2:2754657/1‑150 (MQ=255)
          ccAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCAtt                                                                                                                                          >  1:3216124/1‑150 (MQ=255)
           cAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTa                                                                                                                                         <  1:689340/150‑1 (MQ=255)
            aaaTCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTaa                                                                                                                                        >  2:2122573/1‑150 (MQ=255)
            aaaTCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTaa                                                                                                                                        >  1:3562757/1‑150 (MQ=255)
             aaTCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAAt                                                                                                                                       <  1:5109294/150‑1 (MQ=255)
                ccATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATccc                                                                                                                                    >  1:2771403/1‑150 (MQ=255)
                 cATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCa                                                                                                                                   >  1:4031167/1‑150 (MQ=255)
                            cGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTccc                                                                                                                        >  1:1149636/1‑150 (MQ=255)
                            cGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTccc                                                                                                                        >  2:2840226/1‑150 (MQ=255)
                             gCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCAGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCt                                                                                                                       <  2:3236411/150‑1 (MQ=255)
                                   gACATCACCGATCCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATTCAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTg                                                                                                                 <  1:3341170/150‑1 (MQ=255)
                                   gACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTg                                                                                                                 >  1:4251384/1‑150 (MQ=255)
                                    aCATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTgg                                                                                                                <  1:3280108/150‑1 (MQ=255)
                                      aTCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCt                                                                                                              >  2:3874203/1‑150 (MQ=255)
                                         aCCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGTCATATTCAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTTTCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCtttt                                                                                                           <  2:1473314/150‑1 (MQ=255)
                                          ccGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTg                                                                                                          <  1:4086617/150‑1 (MQ=255)
                                               ccGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAaag                                                                                                     <  2:5175544/150‑1 (MQ=255)
                                                 gTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAaagaa                                                                                                   <  1:2550366/150‑1 (MQ=255)
                                                 gTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAaagaa                                                                                                   >  1:4381459/1‑150 (MQ=255)
                                                      ccTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCGACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAAc                                                                                              <  1:3902389/150‑1 (MQ=255)
                                                      ccTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAAc                                                                                              <  2:3955895/150‑1 (MQ=255)
                                                        tGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTg                                                                                            <  2:108455/150‑1 (MQ=255)
                                                          cGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGat                                                                                          <  2:3504227/150‑1 (MQ=255)
                                                          cGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGat                                                                                          >  1:3799909/1‑150 (MQ=255)
                                                           ggTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGata                                                                                         >  2:206775/1‑150 (MQ=255)
                                                           ggTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGata                                                                                         <  1:2977709/150‑1 (MQ=255)
                                                                 ccATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTAAACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCACCCGTGTTTGGCTTTTTGAAAGAAGAAACTGATATTGccc                                                                                   >  1:847328/1‑150 (MQ=255)
                                                                  cATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCa                                                                                  <  1:976499/150‑1 (MQ=255)
                                                                   aTTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCAt                                                                                 >  1:4728629/1‑150 (MQ=255)
                                                                   aTTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCAt                                                                                 >  1:1553548/1‑150 (MQ=255)
                                                                   aTTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGGTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCAt                                                                                 >  2:406705/1‑150 (MQ=255)
                                                                   aTTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAAACAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCAt                                                                                 <  2:2382009/150‑1 (MQ=255)
                                                                   aTTGAAGCCGAGAAATGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTAAACCTAATTTTACTTATCAGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCAt                                                                                 <  1:1419220/150‑1 (MQ=255)
                                                                    ttGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATg                                                                                <  2:90053/150‑1 (MQ=255)
                                                                     tGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGc                                                                               <  1:3860010/150‑1 (MQ=255)
                                                                            gagaAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCAGCCATTAAACCCAAGAGCGCTCCCTTGTTTGGCTGTTGGAAAGAAGAAACTGATATAGACCATGCTCTAAGc                                                                        >  1:2403380/1‑150 (MQ=255)
                                                                             agaAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCa                                                                       <  1:1709332/150‑1 (MQ=255)
                                                                             agaAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCa                                                                       >  2:3800628/1‑150 (MQ=255)
                                                                             aaaaaGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCa                                                                       >  1:4758378/3‑150 (MQ=255)
                                                                                aaGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTGTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAg                                                                    >  2:41830/1‑150 (MQ=255)
                                                                                 aGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCACATGCTCGAAGCAAAGt                                                                   >  2:4535405/1‑150 (MQ=255)
                                                                                  ggCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCTCCATTAATCCCAAGAGCGCTCCCTTTTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGtt                                                                  <  2:4281445/150‑1 (MQ=255)
                                                                                   gCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGttt                                                                 <  1:4884452/150‑1 (MQ=255)
                                                                                     aaaTACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTGTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGc                                                               >  2:2992825/1‑150 (MQ=255)
                                                                                       aTACAGGTACAAAAATGACCAGAGGTGCAACCAGCTGCATTACACATACACCTAAGTTTCCATATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCGAAGTTTGCAc                                                             >  1:1293471/1‑150 (MQ=255)
                                                                                          cAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTg                                                          <  1:4142213/150‑1 (MQ=255)
                                                                                           aGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGc                                                         >  2:3342915/1‑150 (MQ=255)
                                                                                             gTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCaa                                                       <  1:3055092/150‑1 (MQ=255)
                                                                                             gTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTGTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCaa                                                       >  1:617641/1‑150 (MQ=255)
                                                                                                cAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCCAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAAc                                                    >  2:3726814/1‑150 (MQ=255)
                                                                                                 aaaaaTGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAAcc                                                   >  2:3233770/1‑150 (MQ=255)
                                                                                                 aaaaaTGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAAcc                                                   <  2:1390608/150‑1 (MQ=255)
                                                                                                 aaaaaTGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAAcc                                                   >  1:2230664/1‑150 (MQ=255)
                                                                                                 aaaaaTGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTGTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAAcc                                                   >  1:1081591/1‑150 (MQ=255)
                                                                                                     aTGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCAt                                               <  1:4559717/150‑1 (MQ=255)
                                                                                                       gACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAg                                             <  2:3134153/150‑1 (MQ=255)
                                                                                                        aCCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGc                                            >  1:3159546/1‑150 (MQ=255)
                                                                                                        aCCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGc                                            <  1:2606653/150‑1 (MQ=255)
                                                                                                            gCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAg                                        <  1:2585419/150‑1 (MQ=255)
                                                                                                              ggTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcg                                      >  2:558104/1‑150 (MQ=255)
                                                                                                               gTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcga                                     <  1:216401/150‑1 (MQ=255)
                                                                                                                 gCAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgata                                   <  2:2516900/150‑1 (MQ=255)
                                                                                                                  cAACCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataa                                  >  2:3060454/1‑150 (MQ=255)
                                                                                                                   aaCCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataac                                 >  2:3653952/1‑150 (MQ=255)
                                                                                                                   aaCCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataac                                 >  2:1230279/1‑150 (MQ=255)
                                                                                                                   aaCCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataac                                 >  2:5036210/1‑150 (MQ=255)
                                                                                                                   aaCCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataac                                 <  1:2654812/150‑1 (MQ=255)
                                                                                                                   aaCCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataac                                 >  1:5003743/1‑150 (MQ=255)
                                                                                                                    aCCAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataacg                                <  1:888256/150‑1 (MQ=255)
                                                                                                                     ccAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataacga                               >  2:2007905/1‑150 (MQ=255)
                                                                                                                     ccAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataacga                               >  2:4727838/1‑150 (MQ=255)
                                                                                                                      cAGCTGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataacgat                              <  2:4675179/150‑1 (MQ=255)
                                                                                                                          tGCATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAat                          <  1:2864686/150‑1 (MQ=255)
                                                                                                                            cATTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAatat                        <  2:3216124/150‑1 (MQ=255)
                                                                                                                             aTTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGcgataacgataatatac                       >  2:5104512/1‑145 (MQ=255)
                                                                                                                             aTTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATc                       >  2:2968690/1‑150 (MQ=255)
                                                                                                                             aTTACACTTACACCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATc                       >  2:741234/1‑150 (MQ=255)
                                                                                                                              ttacacttacacCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATcc                      <  2:1018704/150‑1 (MQ=255)
                                                                                                                              ttacacttacacCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCAATGCTCGAAGCTAAGTTTGCACCTGCATAACCGCATAGAAAAGCGATAACGATAAATTTcc                      <  2:1675219/150‑1 (MQ=255)
                                                                                                                                  acttacacCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCaaa                  >  2:1914554/1‑150 (MQ=255)
                                                                                                                                    ttacacCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCACAGCGAGAACGAGAAATATCCCAAAAg                >  1:232544/1‑150 (MQ=255)
                                                                                                                                     tacacCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAgg               <  1:5220336/150‑1 (MQ=255)
                                                                                                                                      acacCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGa              >  2:1761634/1‑150 (MQ=255)
                                                                                                                                        acCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGt            <  1:1819591/150‑1 (MQ=255)
                                                                                                                                        acCTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGt            <  2:2771403/150‑1 (MQ=255)
                                                                                                                                         ccTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTa           >  1:3646444/1‑150 (MQ=255)
                                                                                                                                         ccTAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTa           >  1:248174/1‑150 (MQ=255)
                                                                                                                                           tAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTAtt         >  2:3646945/1‑150 (MQ=255)
                                                                                                                                           tAAGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTAtt         <  2:1553548/150‑1 (MQ=255)
                                                                                                                                            aaGTTTCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTc        >  1:3620501/1‑150 (MQ=255)
                                                                                                                                               tttCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCACAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGa     >  2:1334519/1‑150 (MQ=255)
                                                                                                                                                ttCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGAt    <  2:3562757/150‑1 (MQ=255)
                                                                                                                                                 tCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATTCTCGAAGCAAAGTTTGCACCTGCAAAACCGCAAAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGAtt   >  1:1930143/1‑150 (MQ=255)
                                                                                                                                                 tCCTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGAtt   <  1:4613703/150‑1 (MQ=255)
                                                                                                                                                  ccTTATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGAttt  <  1:2754657/150‑1 (MQ=255)
                                                                                                                                                  ccTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGAttt  <  2:2269573/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GTAGCGGTACCCAAATCCATGCGGCATTCGCCAGCGACATCACCGAACCGTCGGCCTGCGGTACGCCATTGACGCCGAGAAAGGCAAATACAGGTACAAAAATGACCAGCGGTGCAACCAGCTGCATTACACTTACACCTAAGTTTCCTAATCCGCCATTAATCCCAAGAGCGCTCCCTTGTTTGGCTTTTGGAAAGAAGAAACTGATATTGCCCATGCTCGAAGCAAAGTTTGCACCTGCAAAACCGCATAGCAAAGCGATAACGATAAATATCCCAAAAGGAGTATTCGGATTT  >  CP009273/1537657‑1537952

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