Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,538,185 G→T 6.1% L45I (CTT→ATT)  narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,538,1850GT 6.1% 833.1 / 20.1 228L45I (CTT→ATT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base G (107/107);  new base T (6/8);  total (113/115)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

Only 100 of 239 total aligned reads displayed.

CATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATG  >  CP009273/1538039‑1538332
                                                                                                                                                  |                                                                                                                                                   
cATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCaaga                                                                                                                                                    >  1:2071180/1‑148 (MQ=255)
  tAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGACGGCaagaag                                                                                                                                                  >  2:1881799/1‑148 (MQ=255)
   aaaGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGt                                                                                                                                                 <  2:1381695/148‑1 (MQ=255)
     aGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAg                                                                                                                                               <  1:1752186/148‑1 (MQ=255)
       gAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGAc                                                                                                                                             >  1:2509468/1‑148 (MQ=255)
         gTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACaa                                                                                                                                           <  1:1559783/148‑1 (MQ=255)
            gggAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAactg                                                                                                                                        <  2:261082/148‑1 (MQ=255)
              gAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAactgac                                                                                                                                      >  1:2773890/1‑148 (MQ=255)
              gAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAactgac                                                                                                                                      >  2:2227894/1‑148 (MQ=255)
                acacGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAactgactg                                                                                                                                    <  1:3284115/148‑1 (MQ=255)
                   cGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGata                                                                                                                                 <  1:434940/148‑1 (MQ=255)
                   cGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGata                                                                                                                                 >  2:1687989/1‑148 (MQ=255)
                    gCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAATAAGTAGACAACTGACTGatat                                                                                                                                <  2:2954091/148‑1 (MQ=255)
                    gCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGatat                                                                                                                                <  2:493259/148‑1 (MQ=255)
                     cAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATc                                                                                                                               >  1:2591518/1‑148 (MQ=255)
                        taatGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCag                                                                                                                            >  1:1505104/1‑148 (MQ=255)
                        taatGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCag                                                                                                                            <  2:1744847/148‑1 (MQ=255)
                          atGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCagag                                                                                                                          <  2:2793777/148‑1 (MQ=255)
                           tGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAgaga                                                                                                                         <  2:2025746/148‑1 (MQ=255)
                           tGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAgaga                                                                                                                         <  1:3122438/148‑1 (MQ=255)
                                  gAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATttcttc                                                                                                                  <  1:2537130/148‑1 (MQ=255)
                                   aaaCGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCg                                                                                                                 >  1:3040577/1‑148 (MQ=255)
                                   aaaCGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCg                                                                                                                 >  1:3266355/1‑148 (MQ=255)
                                     aCGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAg                                                                                                               <  1:39094/148‑1 (MQ=255)
                                     aCGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAg                                                                                                               <  1:929467/148‑1 (MQ=255)
                                        gAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCaa                                                                                                            >  1:2288789/1‑148 (MQ=255)
                                         aGGGTAATGCTGTTAATAAAAAGAGTTGATCGGTATTTAAATTACAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAat                                                                                                           <  2:1454557/148‑1 (MQ=255)
                                         aGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAat                                                                                                           <  2:361671/148‑1 (MQ=255)
                                         aGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAat                                                                                                           <  2:1574276/148‑1 (MQ=255)
                                          gggTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAata                                                                                                          <  1:2249498/148‑1 (MQ=255)
                                            gTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATg                                                                                                        >  1:2764292/1‑148 (MQ=255)
                                              aaTGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGtt                                                                                                      >  2:409387/1‑148 (MQ=255)
                                               aTGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGttt                                                                                                     >  1:2209229/1‑148 (MQ=255)
                                                tGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGtttt                                                                                                    <  2:2897478/148‑1 (MQ=255)
                                                tGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGtttt                                                                                                    >  2:3060764/1‑148 (MQ=255)
                                                 gCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTc                                                                                                   <  1:2765989/148‑1 (MQ=255)
                                                   ggTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCt                                                                                                 <  1:289361/148‑1 (MQ=255)
                                                      taataaAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTa                                                                                              <  1:866822/148‑1 (MQ=255)
                                                       aataaAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGATATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTAt                                                                                             <  1:2386838/148‑1 (MQ=255)
                                                        ataaAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTAtt                                                                                            >  1:272062/1‑148 (MQ=255)
                                                        ataaAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTAtt                                                                                            <  1:7006/148‑1 (MQ=255)
                                                         taaAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTAttt                                                                                           >  2:1452176/1‑148 (MQ=255)
                                                         taaAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTAttt                                                                                           >  2:2152597/1‑148 (MQ=255)
                                                                agTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGa                                                                                    >  1:365929/1‑148 (MQ=255)
                                                                 gTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGaa                                                                                   <  2:948401/148‑1 (MQ=255)
                                                                  ttGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAg                                                                                  >  1:195488/1‑148 (MQ=255)
                                                                   tGATCGGTAGTAAAATTAAAACCGATTTTATTGAGAGTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAgg                                                                                 <  1:2795545/148‑1 (MQ=255)
                                                                     aTCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAATAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCACAGAAGGcc                                                                               >  2:645999/1‑148 (MQ=255)
                                                                         gtagtaAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGAt                                                                           <  1:535721/148‑1 (MQ=255)
                                                                         gtagtaAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGAt                                                                           <  1:2198800/148‑1 (MQ=255)
                                                                          tagtaAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGTCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGAtt                                                                          <  1:1077885/148‑1 (MQ=255)
                                                                          tagtaAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGAtt                                                                          <  1:3175758/148‑1 (MQ=255)
                                                                            gtaAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAACGCAAGAAGTAGACAACTGACTGAGATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGAtttt                                                                        <  2:619760/148‑1 (MQ=255)
                                                                                aaTTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTgg                                                                    <  2:2511991/148‑1 (MQ=255)
                                                                                aaTTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTgg                                                                    >  1:1864638/1‑148 (MQ=255)
                                                                                     aaaCCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTcc                                                               >  1:773944/1‑148 (MQ=255)
                                                                                     aaaCCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTcc                                                               <  2:1270033/148‑1 (MQ=255)
                                                                                      aaCCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGGATCCAGACACAGAAGTCAAGAAGTAGACAACTGACTGATATCCAGATATTTCTTCGAGCAATATGTTTTCCTTTATTTGCCCAGAAGGCCGGATTTTCTGTTTTCCa                                                              <  2:2544457/148‑1 (MQ=255)
                                                                                      aaCCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCa                                                              >  1:1942700/1‑148 (MQ=255)
                                                                                      aaCCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCa                                                              <  2:2173278/148‑1 (MQ=255)
                                                                                            ttttATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTcgcg                                                        >  1:801484/1‑148 (MQ=255)
                                                                                                  tGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAg                                                  <  2:1864638/148‑1 (MQ=255)
                                                                                                   gagaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGa                                                 <  1:2664738/148‑1 (MQ=255)
                                                                                                     gaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATa                                               <  1:1587534/148‑1 (MQ=255)
                                                                                                       ttAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAAc                                             >  1:923556/1‑148 (MQ=255)
                                                                                                        tAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACg                                            >  2:752289/1‑148 (MQ=255)
                                                                                                         aaCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGa                                           <  2:1922253/148‑1 (MQ=255)
                                                                                                         aaCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGa                                           <  2:1922447/148‑1 (MQ=255)
                                                                                                          aCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGAc                                          >  2:2297718/1‑148 (MQ=255)
                                                                                                           cgttAACTGCGATAAATAGCATCCAGACACAGACGGCAAGAATTAGACAACTGACTGATATCCAGAGATTTATTCGAGAATTATGTTTTCCTTTATTTTCCCAGAAGGCCGTATTTTCTGGTTTCCTGTCGCGAAAAATATAACGACt                                         <  2:1074426/145‑1 (MQ=255)
                                                                                                            ggTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTACACAACTGACTTATATCCAGAGATTTCTTCGAGAAATATGTTTTCCTTTTTTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTa                                        <  1:1719647/148‑1 (MQ=255)
                                                                                                             gTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTAt                                       >  2:784987/1‑148 (MQ=255)
                                                                                                             gTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTAt                                       >  1:439195/1‑148 (MQ=255)
                                                                                                             gTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTAt                                       >  2:191381/1‑148 (MQ=255)
                                                                                                                aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTAtttt                                    >  2:3104012/1‑148 (MQ=255)
                                                                                                                aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTAtttt                                    >  1:735606/1‑148 (MQ=255)
                                                                                                                aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTAtttt                                    >  1:1238147/1‑148 (MQ=255)
                                                                                                                 cTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTAttttt                                   >  2:1535292/1‑148 (MQ=255)
                                                                                                                  tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTtc                                  <  2:48075/148‑1 (MQ=255)
                                                                                                                    cgcTAAATAGCATCCAGACACAGAAGGCAATAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTtctc                                >  1:2336749/1‑148 (MQ=255)
                                                                                                                     gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCa                               <  1:2152597/148‑1 (MQ=255)
                                                                                                                     gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCa                               <  1:1416062/148‑1 (MQ=255)
                                                                                                                     gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCa                               <  2:2459434/148‑1 (MQ=255)
                                                                                                                       tAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCAtt                             >  1:2320747/1‑148 (MQ=255)
                                                                                                                        aaaTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCAttt                            >  2:90873/1‑148 (MQ=255)
                                                                                                                        aaaTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCAttt                            <  1:1499514/148‑1 (MQ=255)
                                                                                                                           tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGc                         >  2:2739785/1‑148 (MQ=255)
                                                                                                                             gCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAg                       >  2:300080/1‑148 (MQ=255)
                                                                                                                              cATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGt                      >  2:768636/1‑148 (MQ=255)
                                                                                                                               aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTg                     >  1:1932058/1‑148 (MQ=255)
                                                                                                                                 ccAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGcc                   >  1:1655601/1‑148 (MQ=255)
                                                                                                                                   aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCat                 <  2:3280363/148‑1 (MQ=255)
                                                                                                                                    gACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGAATTGCCata                >  2:2118639/1‑148 (MQ=255)
                                                                                                                                         aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGtt           >  1:971136/1‑148 (MQ=255)
                                                                                                                                           aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc         >  1:2436576/1‑148 (MQ=255)
                                                                                                                                           aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc         >  1:2579102/1‑148 (MQ=255)
                                                                                                                                           aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc         >  1:2670776/1‑148 (MQ=255)
                                                                                                                                             ggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc       >  1:597282/1‑148 (MQ=255)
                                                                                                                                                aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca    >  1:1165421/1‑148 (MQ=255)
                                                                                                                                                  gaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg  >  2:277697/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CATAAAGGAGTAGGGAACACGCAATAATGCGCCAGAAACGGAGGGTAATGCGGTTAATAAAAAGAGTTGATCGGTAGTAAAATTAAAACCGATTTTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATG  >  CP009273/1538039‑1538332

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: