Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,534,711 G→T 7.2% R713S (CGC→AGC)  narZ ← nitrate reductase 2 (NRZ), alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,534,7110GT 7.2% 761.8 / 22.6 206R713S (CGC→AGC) narZnitrate reductase 2 (NRZ), alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (84/107);  new base T (6/9);  total (90/116)
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.73e-01

Only 100 of 211 total aligned reads displayed.

ATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGT  >  CP009273/1534566‑1534857
                                                                                                                                                 |                                                                                                                                                  
aTCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGcc                                                                                                                                                  >  1:2253430/1‑148 (MQ=255)
aTCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGGTAGAGCGcc                                                                                                                                                  >  1:466442/1‑148 (MQ=255)
  cGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAg                                                                                                                                                >  2:2406066/1‑148 (MQ=255)
  cGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAg                                                                                                                                                <  2:1381846/148‑1 (MQ=255)
      gCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGacaa                                                                                                                                            <  2:3507412/148‑1 (MQ=255)
      gCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGacaa                                                                                                                                            >  1:2763702/1‑148 (MQ=255)
       cAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGacaaa                                                                                                                                           <  1:3067884/148‑1 (MQ=255)
       cAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGacaaa                                                                                                                                           <  2:2669113/148‑1 (MQ=255)
          tttGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGacaaacaa                                                                                                                                        >  1:3407166/1‑148 (MQ=255)
          tttGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGacaaacaa                                                                                                                                        >  1:2409896/1‑148 (MQ=255)
            tGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAAt                                                                                                                                      >  1:595930/1‑148 (MQ=255)
                  tCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTAcgcg                                                                                                                                >  1:2077949/1‑148 (MQ=255)
                     aCTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGgat                                                                                                                             <  1:1175760/148‑1 (MQ=255)
                     aCTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCCCATTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGgat                                                                                                                             >  1:1873242/1‑148 (MQ=255)
                      cttcttCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCTCCAGACAAACAAATTACGCGgatg                                                                                                                            <  1:366695/148‑1 (MQ=255)
                      cttcttCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGgatg                                                                                                                            >  1:364525/1‑148 (MQ=255)
                      cttcttCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGgatg                                                                                                                            <  1:3237550/148‑1 (MQ=255)
                      cttcttCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGgatg                                                                                                                            <  2:3294038/148‑1 (MQ=255)
                       ttcttcCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGgatga                                                                                                                           <  2:3380598/148‑1 (MQ=255)
                       ttcttcCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGgatga                                                                                                                           <  2:3275873/148‑1 (MQ=255)
                          ttcCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATtg                                                                                                                        <  2:21832/148‑1 (MQ=255)
                             cGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATtgctg                                                                                                                     <  2:3465139/148‑1 (MQ=255)
                             cGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATtgctg                                                                                                                     <  2:1633428/148‑1 (MQ=255)
                                    aTCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTAt                                                                                                              >  1:183086/1‑148 (MQ=255)
                                     tCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGAGTATTGCTGCTGCTATc                                                                                                             <  2:1352315/148‑1 (MQ=255)
                                      cccGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCt                                                                                                            <  2:256432/148‑1 (MQ=255)
                                        cGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTgg                                                                                                          <  1:2026656/148‑1 (MQ=255)
                                         gTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGgt                                                                                                         <  2:1282230/148‑1 (MQ=255)
                                         gTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGgt                                                                                                         <  2:1689954/148‑1 (MQ=255)
                                         gTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGgt                                                                                                         <  2:1623191/148‑1 (MQ=255)
                                            gCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGgttgt                                                                                                      <  1:1303142/148‑1 (MQ=255)
                                            gCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGgttgt                                                                                                      <  1:2714332/148‑1 (MQ=255)
                                            gCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGgttgt                                                                                                      <  1:937854/148‑1 (MQ=255)
                                             cTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGgttgtt                                                                                                     <  1:611574/148‑1 (MQ=255)
                                             cTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGgttgtt                                                                                                     <  1:1101393/148‑1 (MQ=255)
                                              tGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTc                                                                                                    <  1:1344420/148‑1 (MQ=255)
                                               ggCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCg                                                                                                   >  2:3193947/1‑148 (MQ=255)
                                                 cACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCa                                                                                                 >  2:688952/1‑148 (MQ=255)
                                                     gAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGcc                                                                                             <  2:1325438/148‑1 (MQ=255)
                                                      aGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCa                                                                                            >  1:821605/1‑148 (MQ=255)
                                                            tCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTa                                                                                      >  2:795154/1‑148 (MQ=255)
                                                             cGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTaa                                                                                     >  1:2290973/1‑148 (MQ=255)
                                                              gCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCTCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGACCATACGTaaa                                                                                    >  2:1968286/1‑148 (MQ=255)
                                                               cccTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAAt                                                                                   <  2:1398457/148‑1 (MQ=255)
                                                                  tGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGc                                                                                >  2:413831/1‑148 (MQ=255)
                                                                          gCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCGGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTc                                                                        <  1:2563259/148‑1 (MQ=255)
                                                                          gCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCGGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTc                                                                        <  1:2562573/148‑1 (MQ=255)
                                                                            tttCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAg                                                                      <  2:443399/148‑1 (MQ=255)
                                                                             ttCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAgc                                                                     >  2:2855003/1‑148 (MQ=255)
                                                                               cGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAgcgc                                                                   <  1:758773/148‑1 (MQ=255)
                                                                                    cccAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTggg                                                              >  2:3072134/1‑148 (MQ=255)
                                                                                      cagcagATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCTCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCg                                                            >  1:170993/1‑148 (MQ=255)
                                                                                      cagcagATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCg                                                            >  2:1048452/1‑148 (MQ=255)
                                                                                       agcagATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCTCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCgg                                                           <  2:3548944/148‑1 (MQ=255)
                                                                                       agcagATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCTCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCgg                                                           <  2:1792723/148‑1 (MQ=255)
                                                                                       agcagATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCgg                                                           >  1:1838547/1‑148 (MQ=255)
                                                                                       agcagATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCgg                                                           >  1:2711360/1‑148 (MQ=255)
                                                                                       agcagATACTTCTGCATATACTCGGGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCTCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCgg                                                           <  2:2065537/148‑1 (MQ=255)
                                                                                          agATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTaa                                                        <  2:2301484/148‑1 (MQ=255)
                                                                                           gATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTaaa                                                       <  2:2539149/148‑1 (MQ=255)
                                                                                           gATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTaaa                                                       <  2:2539290/148‑1 (MQ=255)
                                                                                              aCTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTc                                                    <  2:211783/148‑1 (MQ=255)
                                                                                              aCTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTc                                                    >  1:1771352/1‑148 (MQ=255)
                                                                                               cTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCt                                                   >  1:3438094/1‑148 (MQ=255)
                                                                                               cTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCt                                                   <  1:2697508/148‑1 (MQ=255)
                                                                                                ttCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCtg                                                  <  2:2015888/148‑1 (MQ=255)
                                                                                                 tCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCtgt                                                 >  2:1506951/1‑148 (MQ=255)
                                                                                                 tCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCtgt                                                 >  1:266377/1‑148 (MQ=255)
                                                                                                  cTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCtgtg                                                >  2:2694007/1‑148 (MQ=255)
                                                                                                    gCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTggg                                              >  1:3491/1‑148 (MQ=255)
                                                                                                      atatACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTggggg                                            >  1:520298/1‑148 (MQ=255)
                                                                                                       tataCTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGa                                           >  2:1673347/1‑148 (MQ=255)
                                                                                                        ataCTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGAt                                          >  1:763179/1‑148 (MQ=255)
                                                                                                         taCTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCTCCAGACAAAAAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATa                                         >  2:1331951/1‑148 (MQ=255)
                                                                                                         taCTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATa                                         <  1:2749932/148‑1 (MQ=255)
                                                                                                         taCTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATa                                         <  2:2576954/148‑1 (MQ=255)
                                                                                                         taCTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATa                                         <  2:2599428/148‑1 (MQ=255)
                                                                                                          aCTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCTCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATaa                                        >  1:3283120/1‑148 (MQ=255)
                                                                                                               tGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCgg                                   >  1:943697/1‑148 (MQ=255)
                                                                                                                ggCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGGTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGc                                  <  1:3275439/148‑1 (MQ=255)
                                                                                                                ggCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGc                                  <  1:2132574/148‑1 (MQ=255)
                                                                                                                ggCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGc                                  <  1:3395990/148‑1 (MQ=255)
                                                                                                                 gCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGcc                                 >  1:841954/1‑148 (MQ=255)
                                                                                                                          ggaggaGCCAAGCAGGTTAGAGCGCCAGACAAACAAATGACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCt                        <  2:1048597/148‑1 (MQ=255)
                                                                                                                           gaggatCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGAATGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAACTTCTGTGGGGGATAATCCGGCCTTGTCGGCtt                       >  2:590515/1‑148 (MQ=255)
                                                                                                                            aggagCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTc                      >  1:172020/1‑148 (MQ=255)
                                                                                                                             ggagCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCTGCCTTGTCGGCTTCa                     >  2:813705/1‑148 (MQ=255)
                                                                                                                               agCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGc                   <  1:3543393/148‑1 (MQ=255)
                                                                                                                                gCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCt                  >  1:3056682/1‑148 (MQ=255)
                                                                                                                                 ccAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCtt                 >  2:1731070/1‑148 (MQ=255)
                                                                                                                                 ccAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCtt                 <  2:183891/148‑1 (MQ=255)
                                                                                                                                   aaGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTTCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTa               <  2:2763702/148‑1 (MQ=255)
                                                                                                                                     gCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAAt             <  1:1164719/148‑1 (MQ=255)
                                                                                                                                         gTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCg         >  1:1268312/1‑148 (MQ=255)
                                                                                                                                           tAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAg       >  2:3520773/1‑148 (MQ=255)
                                                                                                                                           tAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAg       >  2:3520980/1‑148 (MQ=255)
                                                                                                                                             gagCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCg     <  1:1615119/148‑1 (MQ=255)
                                                                                                                                                cgcCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGt  <  1:3193947/148‑1 (MQ=255)
                                                                                                                                                cgcCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGt  <  1:101548/148‑1 (MQ=255)
                                                                                                                                                cgcCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGt  <  2:287387/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATCGCTGCAGTTTGCCACTCGACTTCTTCCGGTTTGATCCCGTCGCTGGCACCGAGTTCCTCGCCCTGAATCCCGCTTTCGGTCCCCAGCAGATACTTCTGCATATACTCGTGGCCTTTGCCGGAGGAGCCAAGCAGGTTAGAGCGCCAGACAAACAAATTACGCGGATGATTGCTGCTGCTATCTGGTTGTTCGCAGGCCATACGTAAATCGCCCGATTTCAGCGCCTGGGCGGTAAATTCTGTGGGGGATAATCCGGCCTTGTCGGCTTCAGCTTTAATCCCGAGCGGGT  >  CP009273/1534566‑1534857

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: