Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,382,003 A→G 5.0% A19A (GCA→GCG ompG → outer membrane porin G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,382,0030AG 5.0% 440.3 / 7.6 133A19A (GCA→GCGompGouter membrane porin G
Reads supporting (aligned to +/- strand):  ref base A (52/73);  new base G (3/4);  total (55/78)
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.16e-01

Only 100 of 133 total aligned reads displayed.

AAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCT  >  NC_000913/1381858‑1382149
                                                                                                                                                 |                                                                                                                                                  
aaaTGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAg                                                                                                                                                 >  3:499875/1‑149 (MQ=255)
  aTGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGc                                                                                                                                               <  3:832153/149‑1 (MQ=255)
  aTGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGc                                                                                                                                               <  3:134008/149‑1 (MQ=255)
    gTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATAGCCTGTGCGCAGGcca                                                                                                                                             <  3:302523/149‑2 (MQ=255)
     tCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGa                                                                                                                                            >  1:1049197/1‑149 (MQ=255)
     tCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGa                                                                                                                                            >  4:755528/1‑149 (MQ=255)
     tCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGa                                                                                                                                            <  3:966761/149‑1 (MQ=255)
     tCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGa                                                                                                                                            <  3:795228/149‑1 (MQ=255)
     tCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGa                                                                                                                                            <  1:332456/149‑1 (MQ=255)
     tCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGa                                                                                                                                            <  2:52470/149‑1 (MQ=255)
     tCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGa                                                                                                                                            <  1:556766/149‑1 (MQ=255)
      cATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGccatt                                                                                                                                           <  2:7624/149‑4 (MQ=255)
       aTTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGccattc                                                                                                                                          <  2:552326/149‑5 (MQ=255)
       aTTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAg                                                                                                                                                 <  2:156732/142‑1 (MQ=255)
       aTTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGagg                                                                                                                                          >  3:640117/1‑149 (MQ=255)
       aTTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAAGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGcca                                                                                                                                             <  4:753626/146‑2 (MQ=255)
       aTTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAAGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAGGcca                                                                                                                                             >  3:753626/1‑145 (MQ=255)
        tttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAg                                                                                                                                                 <  4:979883/141‑1 (MQ=255)
        tttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGagga                                                                                                                                         >  1:337552/1‑149 (MQ=255)
         ttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaa                                                                                                                                        <  1:179952/149‑1 (MQ=255)
          tCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                       <  4:268008/149‑1 (MQ=255)
          tCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                       <  3:233098/149‑1 (MQ=255)
          tCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaa                                                                                                                                       >  2:447903/1‑149 (MQ=255)
          tCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAACGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACCCGCAGaggaaa                                                                                                                                       <  2:979089/149‑1 (MQ=255)
            tttGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaagg                                                                                                                                     >  1:902019/1‑149 (MQ=255)
            tttGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaagg                                                                                                                                     >  1:910270/1‑149 (MQ=255)
                   cAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACt                                                                                                                              <  2:1049197/149‑1 (MQ=255)
                    aGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaa                                                                                                                                        <  2:93388/138‑1 (MQ=255)
                          ggAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACtt                                                                                                                       <  2:32631/149‑1 (MQ=255)
                          ggAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACtt                                                                                                                       <  2:236039/149‑1 (MQ=255)
                              aTAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAat                                                                                                                   <  2:619505/149‑1 (MQ=255)
                                aGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAatat                                                                                                                 >  3:996372/1‑149 (MQ=255)
                                aGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAatat                                                                                                                 >  4:418309/1‑149 (MQ=255)
                                  cAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCg                                                                                                               <  1:1018747/149‑1 (MQ=255)
                                   aaTTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCgg                                                                                                              >  2:744687/1‑149 (MQ=255)
                                      tCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGgcgc                                                                                                           <  2:150284/149‑1 (MQ=255)
                                        gCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGa                                                                                                         <  1:744687/149‑1 (MQ=255)
                                         cTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGAt                                                                                                        <  3:458309/149‑1 (MQ=255)
                                         cTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGAt                                                                                                        <  3:912500/149‑1 (MQ=255)
                                          tAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATg                                                                                                       >  1:64244/1‑149 (MQ=255)
                                          tAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATg                                                                                                       >  3:979110/1‑149 (MQ=255)
                                            tatcctgggggAATGCCTTCCCCCTGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTa                                                                                                     >  1:838557/7‑149 (MQ=255)
                                             aTACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTAc                                                                                                    <  3:679925/149‑1 (MQ=255)
                                               aCAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGa                                                                                                  >  4:966343/1‑149 (MQ=255)
                                                ctgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGaa                                                                                                 >  3:435087/3‑149 (MQ=255)
                                                 tgggggAATGCCTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGaaa                                                                                                >  4:325496/2‑149 (MQ=255)
                                                 aGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGaaa                                                                                                >  2:515288/1‑149 (MQ=255)
                                                 aGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGaaa                                                                                                <  1:94881/149‑1 (MQ=255)
                                                    gggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAg                                                                                             >  1:285001/1‑149 (MQ=255)
                                                     ggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGa                                                                                            >  3:64258/1‑149 (MQ=255)
                                                            aTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTc                                                                                     >  4:205425/1‑149 (MQ=255)
                                                            aTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTc                                                                                     >  1:465521/1‑149 (MQ=255)
                                                            aTTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTc                                                                                     >  1:965937/1‑149 (MQ=255)
                                                             ttCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCg                                                                                    <  2:886858/149‑1 (MQ=255)
                                                              tCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGa                                                                                   <  1:85720/149‑1 (MQ=255)
                                                                ccccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAgg                                                                                 <  1:915179/149‑1 (MQ=255)
                                                                  ccAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGt                                                                               >  1:562718/1‑149 (MQ=255)
                                                                      gATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATg                                                                           >  3:393131/1‑149 (MQ=255)
                                                                      gATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATg                                                                           <  1:52120/149‑1 (MQ=255)
                                                                       ataataCAAGGAACAATAAAGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATgg                                                                          >  1:1086198/1‑149 (MQ=255)
                                                                        taataCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGc                                                                         <  2:105263/149‑1 (MQ=255)
                                                                           taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGaa                                                                      <  1:680774/149‑1 (MQ=255)
                                                                           taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGaa                                                                      <  1:1057763/149‑1 (MQ=255)
                                                                           taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGaa                                                                      <  4:64258/149‑1 (MQ=255)
                                                                           taCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGaa                                                                      <  4:52553/149‑1 (MQ=255)
                                                                            aCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTaa                                                                                                                    <  4:599825/102‑1 (MQ=255)
                                                                            aCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTaa                                                                                                                    >  3:599825/1‑102 (MQ=255)
                                                                            aCAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAg                                                                     >  4:495716/1‑149 (MQ=255)
                                                                            aCAAGAAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGCACTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAg                                                                     >  4:922030/1‑149 (MQ=255)
                                                                             cAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGa                                                                    <  1:637827/149‑1 (MQ=255)
                                                                              aaGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaaag                                                                                                                                      <  1:715051/82‑1 (MQ=255)
                                                                                ggAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGatat                                                                 >  1:239175/1‑149 (MQ=255)
                                                                                  aaCAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatgg                                                               <  1:994406/149‑1 (MQ=255)
                                                                                     aataatGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggatg                                                            <  1:838246/149‑1 (MQ=255)
                                                                                     aataatGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggatg                                                            >  2:807795/1‑149 (MQ=255)
                                                                                     aataatGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggatg                                                            <  2:399807/149‑1 (MQ=255)
                                                                                      ataatGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggatgg                                                           <  1:629813/149‑1 (MQ=255)
                                                                                      ataatGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATatggatgg                                                           >  3:994744/1‑149 (MQ=255)
                                                                                        aatGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGc                                                         <  1:236004/149‑1 (MQ=255)
                                                                                         atGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCt                                                        >  2:162911/1‑149 (MQ=255)
                                                                                              aaaGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGa                                                                    >  4:458300/1‑132 (MQ=255)
                                                                                              aaaGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGa                                                                    <  3:458300/132‑1 (MQ=255)
                                                                                               aaGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCTCAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGa                                                  <  2:910270/149‑1 (MQ=255)
                                                                                                 gTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGc                                                <  1:432826/149‑1 (MQ=255)
                                                                                                    attaCCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCtt                                             <  2:838557/149‑1 (MQ=255)
                                                                                                     ttaCCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTc                                            >  2:876964/1‑149 (MQ=255)
                                                                                                      taCCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGTTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCa                                           <  2:366308/149‑1 (MQ=255)
                                                                                                         ccTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTc                                        >  2:636578/1‑149 (MQ=255)
                                                                                                         ccTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTc                                        >  2:714455/1‑149 (MQ=255)
                                                                                                         ccTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTc                                        >  3:755203/1‑149 (MQ=255)
                                                                                                             tACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTAtt                                    <  4:37348/149‑1 (MQ=255)
                                                                                                               ccGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTAtttt                                  <  1:447903/149‑1 (MQ=255)
                                                                                                                 gCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTaa                                <  4:393131/149‑1 (MQ=255)
                                                                                                                 gCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTaa                                >  4:444595/1‑149 (MQ=255)
                                                                                                                   aCTGGTGAAGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATg                              >  2:1044995/1‑149 (MQ=255)
                                                                                                                         gATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCa                        >  3:1058789/1‑149 (MQ=255)
                                                                                                                                   ggAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTg              >  2:243031/1‑149 (MQ=255)
                                                                                                                                     aaTGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGga            <  4:1058789/149‑1 (MQ=255)
                                                                                                                                              cgcACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGGACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGctc   <  4:505775/149‑1 (MQ=255)
                                                                                                                                               gcACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGctct  <  3:799928/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCGCGATGTACGAAATAGAAAACGTCGAGGGTTATGGCGAAGATATGGATGGGCTGGCGGAGCCTTCAGTCTATTTTAATGCCGCCAACGGGCCGTGGAGAATTGCTCT  >  NC_000913/1381858‑1382149

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