Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,200,772 A→G V12A (GTC→GCC)  intE ← e14 prophage; putative integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,200,7720AG98.7% 288.7 / NA 79V12A (GTC→GCC) intEe14 prophage; putative integrase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (40/38);  minor base C (0/1);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCAGCCAGGCGAGTGTTGGCCTCAGTGGCAATAGCGATAGCTTCAGCTTCGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTC  >  NC_000913/1200625‑1200916
                                                                                                                                                   |                                                                                                                                                
tCAGCCAGGCGAGTGTTGGCCTCAGTGGCAATAGCGATAGCTTCAGCTTCGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTg                                                                                                                                               >  1:906465/1‑151 (MQ=255)
  aGCCAGGCGAGTGTTGGCCTCAGTGGCAATAGCGATAGCTTCAGCTTCGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGa                                                                                                                                              >  1:969135/1‑150 (MQ=255)
   gCCAGGCGGGTGTTGGCCTCAGTGGCAATAGCGATAGCTTCAGTTTCTTTTGTCCCAAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGaa                                                                                                                                             <  2:386722/150‑1 (MQ=255)
    ccgggCGAGTGTTGGCCTCAGTGGCAATAGCTATAGCTTCAGCTTCGTTTTTTCCAAAAGAATGAATTTTTCCTTTCATTGAATTCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTg                                                                                                                                               <  2:1202455/144‑1 (MQ=255)
     cAGGCGAGTGTTGGCCTCAGTGGCAATAGCGATAGCTTCAGCTTCGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAAc                                                                                                                                           >  2:399045/1‑150 (MQ=255)
      aGGCGAGTGTTGGCCTCAGTGGCAATAGCGATAGCTTCAGCTTCGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACat                                                                                                                                         >  2:882794/1‑151 (MQ=255)
         cGAGTGTTGGCCTCAGTGGCAATAGCGATAGCTTCAGCTTCGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACattat                                                                                                                                      >  2:1110455/1‑151 (MQ=255)
          gAGTGTTGGCCTCAGTGGCAATAGCGATAGCTTCAGCTTCGTTTGTCCCCAAAGCATGGAATTTTCCTGTCACTGGATGTTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACattat                                                                                                                                      <  1:501206/150‑1 (MQ=255)
           aGTTTTGGCCTCAGTGGCAATAGGAATAGTTTCAGCTTCTTTTGTCCCAAAAGCATGGATTTTTCCTGTCATTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACatta                                                                                                                                       <  2:850250/148‑1 (MQ=255)
            gtgtTGGCCTCAGGGGCAATAGCGATAGCTTCAGCTTCGTTTGTTCCCAAAGTATGGATTTTTCTTATCACTGGAGGATTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACatta                                                                                                                                       <  2:431010/147‑1 (MQ=255)
             tgtTGGCCTCAGTGGCAATAGCGATATCTTCAGCTTCGTTTGTTCCAAAAGCATGAAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCt                                                                                                                                   <  1:1082021/150‑1 (MQ=255)
                 ggCCTCAGTGGCAATAGCGATAGCTTCAGCTTCGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTAcgc                                                                                                                              >  1:721052/1‑151 (MQ=255)
                 ggCCACAGTGGTAATAGCGATTCCTTCTGCTTCTTTTGTTCCCAAAGATGGGAATTTTCTTTTCGCTGCAGGTTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTa                                                                                                                                 <  1:638374/148‑1 (MQ=255)
                     tCAGTGGCAATAGCGATAGCTTCAGCTTCGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTc                                                                                                                          >  1:280347/1‑151 (MQ=255)
                        gTGGAAATAGCGATAGCTTCAGCTTCGTTTTTCCCCAAAGCATGAAATTTTCCTGTCACTGGATGTTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTg                                                                                                                        <  2:959029/150‑1 (MQ=255)
                          ggCAATAGCGATAGCTTCAGCTTCGTTTTTTCCCAAAGCATGGATTTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCt                                                                                                                     <  1:7752/151‑1 (MQ=255)
                            cAATAGCGATAGCTTCAGCTTCGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGc                                                                                                                   >  1:1173871/1‑151 (MQ=255)
                             aaTAGCGATATTTTCAGCTTCTTTTGTCCCCAAAGCATGGATTTTTCTTGTCATTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGc                                                                                                                   <  2:1023700/150‑1 (MQ=255)
                             aaTAGCGATAGCTTGAGTTTTTTTTTTTCCAAAAGCTTGATTTTTTCTTGTAAATGGATGTTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCt                                                                                                                     <  2:986299/148‑1 (MQ=255)
                             aaTAGCGATAGCTTCAGCTTCGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGcc                                                                                                                  >  1:948667/1‑151 (MQ=255)
                              aTAGCGATCGCTTCAGCTTCGTTTGTCCCCAAAGCATGAAATTTTCCTGTCACTGGTTGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCa                                                                                                                 <  2:643485/151‑1 (MQ=255)
                               tAGCGATAGCTTCAGCTTCGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCAt                                                                                                                >  1:105485/1‑151 (MQ=255)
                                  cGATAGCTTCAGCTTCGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTAc                                                                                                             >  2:912647/1‑151 (MQ=255)
                                      agcttcagcttcGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCa                                                                                                          >  2:361483/1‑150 (MQ=255)
                                         ttcagcttcGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAAt                                                                                                      >  1:90488/1‑151 (MQ=255)
                                          tcagtttcGTTTTTTCCCAATGCTTGGAATTTTCCTGTCACTGAATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAAt                                                                                                      <  2:1102291/150‑1 (MQ=255)
                                            agcttcGTTTTTCCCAAAAGCTTGGATTTTTCCTGCCTCTGGATGATTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATc                                                                                                     <  2:577677/149‑1 (MQ=255)
                                                 cGTTTGTCCCCAAAGCATGGAATTTCCCGGTCACTGAATGTTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTc                                                                                                           <  1:101339/138‑1 (MQ=255)
                                                   ttttttCCAAAAGCATGGAATTTTCCTGTCACGGTATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGt                                                                                                  <  2:642340/145‑1 (MQ=255)
                                                   tttGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCaa                                                                                             >  2:786938/1‑150 (MQ=255)
                                                     tGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATACTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCaaaa                                                                                           <  1:877758/150‑1 (MQ=255)
                                                        tCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTACACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTa                                                                                        >  2:741249/1‑150 (MQ=255)
                                                         cccAAAGAATGAAATTTTCCTCTCCTTGGATCATTATCCCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAAt                                                                                      <  1:662631/151‑1 (MQ=255)
                                                            aaaGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATg                                                                                     >  2:100707/1‑149 (MQ=255)
                                                             aaGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAg                                                                                   >  2:450314/1‑150 (MQ=255)
                                                                 atgtatttttCCTGTCTCTGGATGCTTAATCCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCa                                                                                <  2:257695/143‑1 (MQ=255)
                                                                 aTGCAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCAtt                                                                              >  1:559431/1‑151 (MQ=255)
                                                                   ggATTTTTCCTGTCCTTGGATTTTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCAt                                                                               <  2:1158032/148‑1 (MQ=255)
                                                                       ttttCCTGTCACTGGATGTTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCt                                                                         <  2:136637/150‑1 (MQ=255)
                                                                           ccTGTCACTGAGTGTTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGAt                                                                     <  2:435565/150‑1 (MQ=255)
                                                                            cTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGAtt                                                                    >  1:12536/1‑150 (MQ=255)
                                                                                cACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAgg                                                               >  1:546960/1‑151 (MQ=255)
                                                                                 aCTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGgt                                                              <  1:340971/151‑1 (MQ=255)
                                                                                  cTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACa                                                                 <  2:940396/147‑1 (MQ=255)
                                                                                   tgtatccTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAgg                                                               <  2:868350/142‑1 (MQ=255)
                                                                                   gggATGCATATACCGCCAATAGTCTTTATTTACCTTCCGACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTa                                                                   <  2:892871/143‑1 (MQ=255)
                                                                                    gtatgttaataCCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGgt                                                              <  2:984974/140‑1 (MQ=255)
                                                                                    ggATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGtt                                                           >  2:557173/1‑151 (MQ=255)
                                                                                     gATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGtt                                                           >  1:1201432/1‑150 (MQ=255)
                                                                                      aTGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTAc                                                         >  2:569851/1‑151 (MQ=255)
                                                                                         cTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCa                                                       >  1:997672/1‑150 (MQ=255)
                                                                                            ataCCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACt                                                    <  2:111714/150‑1 (MQ=255)
                                                                                               ccGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTccc                                                 <  1:361181/150‑1 (MQ=255)
                                                                                               ccGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTccc                                                 <  1:989729/150‑1 (MQ=255)
                                                                                                cGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTcccc                                                >  1:948584/1‑150 (MQ=255)
                                                                                                    aaTAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCaact                                            >  2:426447/1‑150 (MQ=255)
                                                                                                     aTAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCaactaa                                          <  1:120901/151‑1 (MQ=255)
                                                                                                      tAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCaactaac                                         >  2:791473/1‑151 (MQ=255)
                                                                                                       aGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCaactaac                                         <  1:361495/150‑1 (MQ=255)
                                                                                                        gACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTc                                       <  2:505664/151‑1 (MQ=255)
                                                                                                           tttatttaCCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGc                                    >  1:231126/1‑151 (MQ=255)
                                                                                                            ttatttaCCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCgg                                     <  1:1079619/149‑1 (MQ=255)
                                                                                                             tatttaCCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGt                                  >  2:127355/1‑151 (MQ=255)
                                                                                                             tatttaCCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGt                                  >  2:1012102/1‑151 (MQ=255)
                                                                                                              atttaCCTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGt                                  <  2:948665/150‑1 (MQ=255)
                                                                                                                   ccTTCCTACTATAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCAc                            >  1:1175608/1‑151 (MQ=255)
                                                                                                                          actaTAAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACt                       <  1:882798/149‑1 (MQ=255)
                                                                                                                            tataAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGcc                   >  2:504911/1‑151 (MQ=255)
                                                                                                                             ataAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGcc                   >  1:1173881/1‑150 (MQ=255)
                                                                                                                             ataAAGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGcc                   <  2:144606/150‑1 (MQ=255)
                                                                                                                                aaGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATaa               <  1:88939/151‑1 (MQ=255)
                                                                                                                                 aGCGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAAc              >  2:36743/1‑151 (MQ=255)
                                                                                                                                   cGGATATAAGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCg            >  1:749795/1‑151 (MQ=255)
                                                                                                                                        ataATTTAGGGCCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTATTGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCGAACTAACCCGGCGTCCGCACGCACTCGCCATAACCGGCCta       <  2:628194/151‑2 (MQ=255)
                                                                                                                                           aGTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGtt    >  2:949197/1‑151 (MQ=255)
                                                                                                                                            gTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGttt   >  1:1124214/1‑151 (MQ=255)
                                                                                                                                            gTTAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGttt   >  1:57070/1‑151 (MQ=255)
                                                                                                                                             ttAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTc  >  1:655920/1‑151 (MQ=255)
                                                                                                                                              tAGGGGCTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTc  >  2:848493/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TCAGCCAGGCGAGTGTTGGCCTCAGTGGCAATAGCGATAGCTTCAGCTTCGTTTGTTCCCAAAGCATGGAATTTTCCTGTCACTGGATGCTTATACCGCCAATAGACTTTATTTACCTTCCTACTATAAAGCGGATATAAGTTAGGGACTGAAACATTATTCTTACGCGGTCTGGCTGCCATTACTCAAAATCCGTTGCAAAAGTAATGAGTCATTTTTCTTGATTACAGGTGTTACCAACTCCCCAACTAACTCGGCGTCCTCACGCACTCGCCATAACCGGCCTTGTTTC  >  NC_000913/1200625‑1200916

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