Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 577,085 G→A intergenic (‑260/‑313) nmpC ← / → essD DLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage/DLP12 prophage; putative phage lysis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913577,0850GA99.1% 459.9 / NA 115intergenic (‑260/‑313)nmpC/essDDLP12 prophage; truncated outer membrane porin (pseudogene);IS, phage, Tn; Phage or Prophage Related; outer membrane porin protein; locus of qsr prophage/DLP12 prophage; putative phage lysis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (50/64);  minor base T (0/1);  total (50/65)
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.99e-01

Only 100 of 115 total aligned reads displayed.

GAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGGCTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTG  >  NC_000913/576939‑577232
                                                                                                                                                  |                                                                                                                                                   
gtaatAGATCTCTATTTTTATAAGAACTTTTTTCTATCGAAACCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGAATAAATTTAAACCGTCTa                                                                                                                                                  <  1:1014063/146‑1 (MQ=255)
gAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTaa                                                                                                                                                 <  1:633715/151‑1 (MQ=255)
   gTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTaaaaa                                                                                                                                              <  2:1583059/151‑1 (MQ=255)
     aGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATa                                                                                                                                            <  2:903805/151‑1 (MQ=255)
           ctATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATaa                                                                                                                                           <  2:116958/146‑1 (MQ=255)
             aTTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTgg                                                                                                                                     <  2:41980/150‑1 (MQ=255)
             aTTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTg                                                                                                                                      <  2:675995/149‑1 (MQ=255)
               ttttATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGaa                                                                                                                                   <  1:1418078/150‑1 (MQ=255)
               ttttATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGaa                                                                                                                                   <  1:602015/150‑1 (MQ=255)
                   aTCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAatta                                                                                                                               <  1:1189105/150‑1 (MQ=255)
                       gAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAattatt                                                                                                                             <  1:1122940/148‑1 (MQ=255)
                       gAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAattat                                                                                                                              <  1:1716925/147‑1 (MQ=255)
                       gAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTAtttt                                                                                                                           <  2:342733/150‑1 (MQ=255)
                        aaCTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTAtttttt                                                                                                                         <  2:1138016/151‑1 (MQ=255)
                        aaCTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTAtttttt                                                                                                                         <  2:1125011/151‑1 (MQ=255)
                        aaCTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTAtttttt                                                                                                                         <  1:159058/151‑1 (MQ=255)
                         aCTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTg                                                                                                                        <  1:866451/151‑1 (MQ=255)
                          cTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTgg                                                                                                                       >  1:1353298/1‑151 (MQ=255)
                            tttttCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTgg                                                                                                                       <  2:294652/149‑1 (MQ=255)
                            tttttCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGtt                                                                                                                     >  1:94463/1‑151 (MQ=255)
                             ttttCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGtt                                                                                                                     <  1:1559777/150‑1 (MQ=255)
                             ttttCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGtt                                                                                                                     <  2:483869/150‑1 (MQ=255)
                              tttCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTa                                                                                                                    <  1:1049411/150‑1 (MQ=255)
                              tttCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTAt                                                                                                                   >  2:1583347/1‑151 (MQ=255)
                                tCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTAttt                                                                                                                 >  1:1022268/1‑151 (MQ=255)
                                tCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTAtt                                                                                                                  <  2:151658/150‑1 (MQ=255)
                                 cTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATttg                                                                                                                >  1:566012/1‑151 (MQ=255)
                                  tATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATttgt                                                                                                               >  1:905901/1‑151 (MQ=255)
                                   aTCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTAttta                                                                                                                <  1:774328/149‑2 (MQ=255)
                                   aTCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATttgtt                                                                                                              <  1:1618051/151‑1 (MQ=255)
                                   aTCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATttgtt                                                                                                              <  1:1266936/151‑1 (MQ=255)
                                   aTCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATttgt                                                                                                               <  2:167303/150‑1 (MQ=255)
                                    tCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATttgttg                                                                                                             >  2:1209503/1‑151 (MQ=255)
                                     cGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATttgttg                                                                                                             <  1:360490/150‑1 (MQ=255)
                                        aTCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATttgttg                                                                                                             >  2:970353/1‑147 (MQ=255)
                                         tCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGAtt                                                                                                        >  1:1515927/1‑151 (MQ=255)
                                          ccTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGAttt                                                                                                       >  1:1234170/1‑151 (MQ=255)
                                             aTTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGc                                                                                                     <  1:540649/150‑1 (MQ=255)
                                             aTTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCt                                                                                                    >  2:1013001/1‑151 (MQ=255)
                                                cATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTAt                                                                                                 <  1:1085499/151‑1 (MQ=255)
                                                     ctctTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTAtt                                                                                            >  1:1093023/1‑151 (MQ=255)
                                                      tctTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTAttt                                                                                           >  1:1067486/1‑151 (MQ=255)
                                                         tGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTg                                                                                          <  1:429569/149‑1 (MQ=255)
                                                             tGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTGTTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAgtgg                                                                                    >  2:1660070/1‑151 (MQ=255)
                                                             tGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAgtgg                                                                                    >  2:739062/1‑151 (MQ=255)
                                                             tGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAAACGACTAAAAATAACACTGGTAATTATTTTTTGGTTATTTGTTGAGATTTGATTATGTTTTTGTAgtg                                                                                     >  1:994783/1‑150 (MQ=255)
                                                              gAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAgtggt                                                                                   <  1:465119/151‑1 (MQ=255)
                                                              gAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAgtggt                                                                                   <  2:473879/151‑1 (MQ=255)
                                                                ataaaaataaaTCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGgtgt                                                                                 >  1:982759/1‑151 (MQ=255)
                                                                  aaaaataaaTCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTAtt                                                                                            <  1:1073277/138‑1 (MQ=255)
                                                                  aaaaataaaTCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTAtt                                                                                            >  2:1073363/1‑138 (MQ=255)
                                                                    aaataaatCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTc                                                                             <  1:1663960/151‑1 (MQ=255)
                                                                    aaataaatCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTc                                                                             <  1:1644298/151‑1 (MQ=255)
                                                                     aataaatCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCa                                                                            <  2:1495512/151‑1 (MQ=255)
                                                                      ataaatCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCa                                                                            >  1:1166990/1‑150 (MQ=255)
                                                                        aaatCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAAt                                                                          <  2:1203237/150‑1 (MQ=255)
                                                                         aatCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATAc                                                                        <  1:737172/151‑1 (MQ=255)
                                                                                tAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCGATACTCGGTa                                                                  <  1:791794/150‑1 (MQ=255)
                                                                                  gCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGc                                                                >  2:1545063/1‑150 (MQ=255)
                                                                                  gCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGc                                                                >  1:1458092/1‑150 (MQ=255)
                                                                                   ccAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATGACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCAt                                                              <  2:963568/151‑1 (MQ=255)
                                                                                    cAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCAtt                                                             <  2:154764/151‑1 (MQ=255)
                                                                                     aaTTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCAtt                                                             >  1:1627600/1‑150 (MQ=255)
                                                                                       tttATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATtctc                                                          >  2:1682489/1‑151 (MQ=255)
                                                                                       tttATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATtctc                                                          >  2:9140/1‑151 (MQ=255)
                                                                                        ttATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATtctct                                                         >  1:932980/1‑151 (MQ=255)
                                                                                          atatTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCa                                                       <  1:535147/151‑1 (MQ=255)
                                                                                          atatTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCa                                                       <  2:1303556/151‑1 (MQ=255)
                                                                                               aaCGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAatat                                                  >  1:34762/1‑151 (MQ=255)
                                                                                                 cGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCa                                                <  1:1209410/151‑1 (MQ=255)
                                                                                                   gCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCAtt                                              <  2:596686/151‑1 (MQ=255)
                                                                                                   gCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCAtt                                              <  2:1412582/151‑1 (MQ=255)
                                                                                                   gCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCAtt                                              >  1:1633085/1‑151 (MQ=255)
                                                                                                    cTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCAttt                                             >  2:1157143/1‑151 (MQ=255)
                                                                                                       ttatttatAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGt                                          >  2:246846/1‑151 (MQ=255)
                                                                                                         atttatAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTg                                         <  2:1562360/150‑1 (MQ=255)
                                                                                                         atttatAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTg                                         <  1:866006/150‑1 (MQ=255)
                                                                                                          tttatAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGt                                       >  1:325606/1‑151 (MQ=255)
                                                                                                             ataAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGttt                                     >  2:536422/1‑150 (MQ=255)
                                                                                                               aaGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTtacg                                  <  1:759921/151‑1 (MQ=255)
                                                                                                               aaGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTtac                                   <  2:326861/150‑1 (MQ=255)
                                                                                                                aGTGCTCTATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTtacgt                                 <  2:327342/151‑1 (MQ=255)
                                                                                                                      ctATAATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTAcgtacgta                            >  2:248060/1‑150 (MQ=255)
                                                                                                                        ataATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTaaaa                         <  2:975087/151‑1 (MQ=255)
                                                                                                                        ataATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTaaa                          >  2:1290334/1‑150 (MQ=255)
                                                                                                                        ataATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTaaa                          >  1:1660264/1‑150 (MQ=255)
                                                                                                                        ataATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTaa                           <  2:1627701/149‑1 (MQ=255)
                                                                                                                        ataATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTaa                           <  1:827642/149‑1 (MQ=255)
                                                                                                                         taATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTaaaa                         >  1:1684578/1‑150 (MQ=255)
                                                                                                                         taATTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTaaaa                         >  1:1668591/1‑150 (MQ=255)
                                                                                                                          aaTTTGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTaaaaa                        >  1:179404/1‑150 (MQ=255)
                                                                                                                              tGAAGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTg                    <  1:725123/150‑1 (MQ=255)
                                                                                                                                aaGGTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGtt                 >  1:1427237/1‑151 (MQ=255)
                                                                                                                                   gTTCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATg              <  2:1239306/151‑1 (MQ=255)
                                                                                                                                    ttCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGc             >  1:702256/1‑151 (MQ=255)
                                                                                                                                     tCAATTTAAACCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGc             <  2:1548129/150‑1 (MQ=255)
                                                                                                                                             aaCCGACTAAAAGTAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTaa       >  1:240557/1‑148 (MQ=255)
                                                                                                                                             aaCCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAag      >  2:58006/1‑149 (MQ=255)
                                                                                                                                             aaCCGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAag      >  1:478071/1‑149 (MQ=255)
                                                                                                                                               ccGACTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTg  <  2:589007/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GAAGTAGATCTCTATTTTTATCTGAACTTTTTTCTATCGAATCCTATTCATGGCTCTTGGCTGAATAAAAATAAATCTATTAGCCAATTTATATTAACGGCTGTTATTTATAAGTGCTCTATAATTTGAAGGTTCAATTTAAACCGGCTAAAAATAACACTGGAAATTATTTTTTGGTTATTTGTTGAGATTTGCTTATGTATTTGTAGTGGTGTTTTCAATACTCGGTAGCATTCTCTCAAATATCATTTAGTGGTTTACGTACGTAAAAAATTGGTTATGCTGTTAAGAGTG  >  NC_000913/576939‑577232

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