Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 137,057 G→A intergenic (‑140/‑26) yacC ← / → yacK hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913137,0570GA98.0% 690.5 / ‑3.0 203intergenic (‑140/‑26)yacC/yacKhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base A (71/128);  total (71/132)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 215 total aligned reads displayed.

AGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGT  >  NC_000913/136967‑137146
                                                                                          |                                                                                         
agGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCAcc                                                                                         <  2:4692562/93‑1 (MQ=255)
agGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCAcc                                                                                         <  2:181503/93‑1 (MQ=255)
agGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCAcc                                                                                         <  1:552133/93‑1 (MQ=255)
  gCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTTACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTTTTTCACCtc                                                                                       >  2:1059208/1‑92 (MQ=255)
  gCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTg                                                                                       <  1:876598/93‑1 (MQ=255)
   cgacgaCCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGc                                                                                      <  1:1341043/93‑1 (MQ=255)
    gacgacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCt                                                                                     <  1:2802029/93‑1 (MQ=255)
       gacCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAAGGGAAGGACTATGCTCAACCTTTGATTTTGTTTCACCCGCTTaa                                                                                  >  2:2145373/1‑93 (MQ=255)
        acCATTATTGTCGATTAAATTGTGTCTGCGGCTTTTCCCTCCCCTGAGGGGAAGGGCCATGGTCAACGTTTGGATTTGGTTCATCTGCGTAAg                                                                                 >  2:693862/1‑93 (MQ=255)
         ccATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGa                                                                                <  1:2598973/93‑1 (MQ=255)
           attattGTCGATTAAATTGTGTCTGCGGCTTTAACTTCCCGTAAGGGGGAGGGCTATGCTCAACGTTTGATTTTTTTTCCACTGGTTAaaaaa                                                                              >  2:4354914/1‑92 (MQ=255)
           attattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAAt                                                                              <  1:1352051/93‑1 (MQ=255)
           attattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAAt                                                                              <  2:3141353/93‑1 (MQ=255)
             tattGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATaa                                                                            >  2:159622/1‑93 (MQ=255)
               ttGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAgg                                                                          >  2:1558276/1‑93 (MQ=255)
               ttGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAgg                                                                          <  1:947719/93‑1 (MQ=255)
                tGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGa                                                                         >  1:746516/1‑93 (MQ=255)
                tGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGa                                                                         <  2:3598288/93‑1 (MQ=255)
                tGTCGATTAAAGTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGa                                                                         <  1:4693001/93‑1 (MQ=255)
                  tCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGaaa                                                                       <  1:3136779/93‑1 (MQ=255)
                   cGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAAt                                                                      <  2:1752531/93‑1 (MQ=255)
                   cGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAAt                                                                      <  1:66169/93‑1 (MQ=255)
                    gATTAAATTGTGTCTGCGGCTTGGCCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATa                                                                     <  2:1442171/93‑1 (MQ=255)
                    gATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATa                                                                     <  1:737066/93‑1 (MQ=255)
                     aTTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATaa                                                                    <  2:3908047/93‑1 (MQ=255)
                     aTTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATaa                                                                    <  1:4318024/93‑1 (MQ=255)
                     aTTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATaa                                                                    <  2:3451314/93‑1 (MQ=255)
                       tAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACt                                                                  >  1:2864910/1‑93 (MQ=255)
                        aaaTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTa                                                                 <  1:2989629/93‑1 (MQ=255)
                        aaaTTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTa                                                                 <  1:2536003/93‑1 (MQ=255)
                           ttGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGc                                                              <  1:4432388/93‑1 (MQ=255)
                            tgtgTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCa                                                             <  2:364511/93‑1 (MQ=255)
                             gtgtCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCaa                                                            <  1:491981/93‑1 (MQ=255)
                               gtCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAAcg                                                          <  1:1674223/93‑1 (MQ=255)
                                tCTGCGGTTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAAcgt                                                         <  1:4000112/93‑1 (MQ=255)
                                tCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAAcgt                                                         <  2:2014424/93‑1 (MQ=255)
                                tCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAAcgt                                                         >  1:4543482/1‑93 (MQ=255)
                                tCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAAcgt                                                         <  2:3542072/93‑1 (MQ=255)
                                tCTGCGGATTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAAcgt                                                         <  1:2918857/93‑1 (MQ=255)
                                   gCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAAcgtcgt                                                      <  2:4204358/93‑1 (MQ=255)
                                     ggCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAAcgtcgtta                                                    >  2:2377595/1‑91 (MQ=255)
                                      gCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGAt                                                   <  1:1709239/93‑1 (MQ=255)
                                         tGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTc                                                >  2:1213228/1‑93 (MQ=255)
                                          gACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCt                                               >  2:2921367/1‑93 (MQ=255)
                                          gACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCt                                               >  2:3126676/1‑93 (MQ=255)
                                             cTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaa                                            >  1:2030318/1‑93 (MQ=255)
                                              ttCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGCAATAACTATGCAACGTCGTGCTTTCTTaaa                                           >  1:401785/1‑93 (MQ=255)
                                              ttCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaaa                                           >  1:3713675/1‑93 (MQ=255)
                                              ttCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaaa                                           >  1:2224833/1‑93 (MQ=255)
                                              ttCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaaa                                           >  1:267118/1‑93 (MQ=255)
                                               tCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaaaa                                          <  2:2990143/93‑1 (MQ=255)
                                               tCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaaaa                                          >  1:4090876/1‑93 (MQ=255)
                                               tCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTaaaa                                          >  2:2061001/1‑93 (MQ=255)
                                                 ccGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAata                                        >  2:2605072/1‑93 (MQ=255)
                                                 ccGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAata                                        <  1:163554/93‑1 (MQ=255)
                                                   gTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATAtt                                      <  2:4273434/93‑1 (MQ=255)
                                                   gTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATAtt                                      <  2:415259/93‑1 (MQ=255)
                                                   gTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATAtt                                      >  2:454691/1‑93 (MQ=255)
                                                   gTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATAtt                                      <  2:919929/93‑1 (MQ=255)
                                                   gTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATAtt                                      >  2:1537145/1‑93 (MQ=255)
                                                    tAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTc                                     <  2:1412636/93‑1 (MQ=255)
                                                    tAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTc                                     <  2:2351920/93‑1 (MQ=255)
                                                     aaGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTcc                                    <  2:1953491/93‑1 (MQ=255)
                                                     aaGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTcc                                    <  1:124766/93‑1 (MQ=255)
                                                        gggAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTc                                 <  1:4380175/93‑1 (MQ=255)
                                                        gggAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTc                                 <  1:1873622/93‑1 (MQ=255)
                                                        gggAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTc                                 <  2:444317/93‑1 (MQ=255)
                                                        gggAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTc                                 <  1:706613/93‑1 (MQ=255)
                                                          gAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTcgc                               <  2:412286/93‑1 (MQ=255)
                                                          gAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTcgc                               <  2:4041451/93‑1 (MQ=255)
                                                          gAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTcgc                               <  2:254654/93‑1 (MQ=255)
                                                          gAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTcgc                               <  1:3332568/93‑1 (MQ=255)
                                                          gAAGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTcgc                               >  1:3979653/1‑93 (MQ=255)
                                                           aaGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTcgcg                              <  2:1400381/93‑1 (MQ=255)
                                                            aGGACTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCgcgc                             <  1:1229921/93‑1 (MQ=255)
                                                               aCTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTgg                          <  2:2477975/93‑1 (MQ=255)
                                                                cTATGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTggg                         >  1:4508933/1‑93 (MQ=255)
                                                                  aTGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGCCGTGATTTCTTAAAATATTCCGTCGCGCTGGgtg                       >  2:1290424/1‑93 (MQ=255)
                                                                   tGCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGgtgt                      <  1:848460/93‑1 (MQ=255)
                                                                    gCTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGtgtg                     <  2:3329106/93‑1 (MQ=255)
                                                                     cTCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTgg                    >  1:2903617/1‑93 (MQ=255)
                                                                      tCAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGc                   <  2:3766044/93‑1 (MQ=255)
                                                                       cAACGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCt                  <  1:4747960/93‑1 (MQ=255)
                                                                        aaCGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCtt                 >  1:3548782/1‑93 (MQ=255)
                                                                        aaCGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCtt                 <  2:2289656/93‑1 (MQ=255)
                                                                        aaCGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCtt                 >  2:4639147/1‑93 (MQ=255)
                                                                         aCGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTc                >  1:2943562/1‑93 (MQ=255)
                                                                         aCGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTc                <  1:2404829/93‑1 (MQ=255)
                                                                         aCGTTTGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTc                <  1:2063613/93‑1 (MQ=255)
                                                                            tttGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCggt             >  1:1558648/1‑92 (MQ=255)
                                                                            tttGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGc             >  2:3193537/1‑93 (MQ=255)
                                                                            tttGATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGc             >  2:3086088/1‑93 (MQ=255)
                                                                              tGATTTTGTTTCCCCTGCTTAAGAAGAGAGAAATAACTATGCAACGTCGTGAGTTCTGAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCtt           <  1:2274051/93‑1 (MQ=255)
                                                                               gATTTTGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCttt          <  2:2405113/93‑1 (MQ=255)
                                                                                 ttttGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGc        >  1:2707824/1‑93 (MQ=255)
                                                                                 ttttGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGc        >  1:4669198/1‑93 (MQ=255)
                                                                                 ttttGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGc        >  2:852293/1‑93 (MQ=255)
                                                                                  tttGTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGcc       >  1:1064400/1‑93 (MQ=255)
                                                                                     gTTTCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCt    >  2:1871615/1‑93 (MQ=255)
                                                                                       ttCACCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCtgt  <  1:1459179/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGGCGACGACCATTATTGTCGATTAAATTGTGTCTGCGGCTTGACCTTCCCGTAAGGGGAAGGACTATGCTCAACGTTTGATTTTGTTTCGCCTGCTTAAGAATAAGGAAATAACTATGCAACGTCGTGATTTCTTAAAATATTCCGTCGCGCTGGGTGTGGCTTCGGCTTTGCCGCTGT  >  NC_000913/136967‑137146

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