Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 740,004 G→A M166I (ATG→ATA phr → deoxyribodipyrimidine photolyase, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913740,0040GA98.0% 347.0 / ‑2.2 101M166I (ATG→ATAphrdeoxyribodipyrimidine photolyase, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (53/46);  total (53/48)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

Only 100 of 101 total aligned reads displayed.

TCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACT  >  NC_000913/739907‑740100
                                                                                                 |                                                                                                
tCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATAcc                                                                                                <  2:389081/100‑1 (MQ=255)
 ccTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCg                                                                                               >  2:2558274/1‑100 (MQ=255)
 ccTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCg                                                                                               >  1:2045465/1‑100 (MQ=255)
  cTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCgg                                                                                              <  1:796255/100‑1 (MQ=255)
    gccgccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAg                                                                                            >  2:460910/1‑100 (MQ=255)
    gccgccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGGCTTTACGCCTTTTACGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAg                                                                                            >  1:1364776/1‑100 (MQ=255)
     ccgccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGt                                                                                           >  2:2116002/1‑100 (MQ=255)
        ccTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCg                                                                                        <  1:2608473/100‑1 (MQ=255)
         cTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGt                                                                                       <  2:2334481/100‑1 (MQ=255)
           ggCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCg                                                                                     >  1:260020/1‑100 (MQ=255)
            gcgcGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGc                                                                                    <  1:361374/100‑1 (MQ=255)
               cgGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTgc                                                                                 <  1:2634112/100‑1 (MQ=255)
                   gatgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCa                                                                             >  2:948953/1‑100 (MQ=255)
                     tgaCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCaaa                                                                           >  1:2531469/1‑100 (MQ=255)
                       aCCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAg                                                                         <  2:2657373/100‑1 (MQ=255)
                        gcgttAAGCACGAGAGGTACAAAGTGTTTACGCCTTTTAAGAATGCCGGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGt                                                                        <  2:1877560/96‑1 (MQ=255)
                        ccGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGt                                                                        >  2:2320305/1‑100 (MQ=255)
                         cGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGtt                                                                       <  2:1842698/100‑1 (MQ=255)
                          ggTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTc                                                                      >  1:2301030/1‑100 (MQ=255)
                            tAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgt                                                                    <  1:263764/100‑1 (MQ=255)
                            tAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgt                                                                    <  2:77236/100‑1 (MQ=255)
                            tAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgt                                                                    <  1:1695141/100‑1 (MQ=255)
                            tAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgt                                                                    >  1:2029848/1‑100 (MQ=255)
                            tAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgt                                                                    <  1:1918157/100‑1 (MQ=255)
                              aTCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGATGCGCCAAAAGTTCgtag                                                                  >  2:1317404/1‑100 (MQ=255)
                               tCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgtagt                                                                 <  1:2182805/100‑1 (MQ=255)
                               tCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCgtagt                                                                 >  2:2159988/1‑100 (MQ=255)
                                 aCGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGtagtag                                                               >  1:2325637/1‑100 (MQ=255)
                                  cGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGc                                                              <  2:1847599/100‑1 (MQ=255)
                                   gagaTGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCg                                                             <  2:1877526/100‑1 (MQ=255)
                                     gaTGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGa                                                           <  1:2752987/100‑1 (MQ=255)
                                       tGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                                         >  1:1244185/1‑100 (MQ=255)
                                       tGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                                         >  2:517785/1‑100 (MQ=255)
                                       tGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                                         >  2:516823/1‑100 (MQ=255)
                                       tGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATc                                                         <  2:1538654/100‑1 (MQ=255)
                                        gTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCg                                                        <  1:1239932/100‑1 (MQ=255)
                                         tACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGa                                                       >  2:288031/1‑100 (MQ=255)
                                          aCAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGAt                                                      >  2:2464454/1‑100 (MQ=255)
                                          aCAAAGTCTTTACGCCTTGTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGAt                                                      <  2:1364767/100‑1 (MQ=255)
                                           cAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATa                                                     >  2:2722510/1‑100 (MQ=255)
                                            aaaGTCTTTACGCTTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATag                                                    <  1:1609393/100‑1 (MQ=255)
                                              aGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATagag                                                  <  1:2167762/100‑1 (MQ=255)
                                               gTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGc                                                 <  1:443242/100‑1 (MQ=255)
                                               gTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGc                                                 <  1:2320323/100‑1 (MQ=255)
                                               gTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCAGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGc                                                 >  2:1457910/1‑100 (MQ=255)
                                                tCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGcc                                                <  1:2160006/100‑1 (MQ=255)
                                                 cTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGccc                                               <  1:372309/100‑1 (MQ=255)
                                                 cTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGccc                                               <  2:2531450/100‑1 (MQ=255)
                                                 cTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGccc                                               >  1:2555800/1‑100 (MQ=255)
                                                   ttACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTc                                             >  2:1961418/1‑100 (MQ=255)
                                                    tACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCg                                            >  2:478557/1‑100 (MQ=255)
                                                    tACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCg                                            >  1:2393611/1‑100 (MQ=255)
                                                     aCGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGc                                           <  1:1317412/100‑1 (MQ=255)
                                                      cGCCTTGTAAGAATGCCTGGCTGAAACGGCTGCGGGAATGGATACCGGAGGGCGTCGCGGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGcc                                          <  2:2507323/100‑1 (MQ=255)
                                                       gCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGcca                                         >  2:960969/1‑100 (MQ=255)
                                                       gCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGcca                                         <  2:260017/100‑1 (MQ=255)
                                                       gCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGcca                                         >  1:189085/1‑100 (MQ=255)
                                                          ctttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatcc                                      <  2:648637/99‑1 (MQ=255)
                                                           tttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatcca                                     <  2:325897/100‑1 (MQ=255)
                                                           tttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatcca                                     <  2:880696/100‑1 (MQ=255)
                                                            ttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                                    >  2:2649805/1‑100 (MQ=255)
                                                            ttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                                    >  2:1803201/1‑100 (MQ=255)
                                                            ttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                                    <  1:581158/100‑1 (MQ=255)
                                                            ttAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGccatccat                                    <  1:1406420/100‑1 (MQ=255)
                                                              aaGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTa                                  <  2:2029831/100‑1 (MQ=255)
                                                               aGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTAc                                 >  1:1404527/1‑100 (MQ=255)
                                                               aGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTAc                                 <  1:1853120/100‑1 (MQ=255)
                                                                 aaTGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGc                               <  2:1477973/100‑1 (MQ=255)
                                                                  aTGCCTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTAGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACgtt                              >  1:1215456/1‑98 (MQ=255)
                                                                     ccTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCTATCCATTACGCTGaa                           >  1:2725928/1‑100 (MQ=255)
                                                                     ccTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGaa                           >  2:1103002/1‑100 (MQ=255)
                                                                      cTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAAt                          <  1:2549371/100‑1 (MQ=255)
                                                                      cTGGCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAAt                          <  1:2645210/100‑1 (MQ=255)
                                                                        ngCTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTa                        >  1:2163379/2‑100 (MQ=255)
                                                                          cTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATc                      >  2:567125/1‑100 (MQ=255)
                                                                          cTGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATc                      >  1:569080/1‑100 (MQ=255)
                                                                           tGAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATc                      >  1:1301195/1‑99 (MQ=255)
                                                                            gAAACGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCt                    >  2:780960/1‑100 (MQ=255)
                                                                              aaCGGCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCg                  >  1:967691/1‑100 (MQ=255)
                                                                                  gCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAg              >  2:2099547/1‑100 (MQ=255)
                                                                                  gCTGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAg              <  1:2204396/100‑1 (MQ=255)
                                                                                    tGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTc            >  2:2173523/1‑100 (MQ=255)
                                                                                    tGCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTc            >  2:1957372/1‑100 (MQ=255)
                                                                                     gCGGGAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCt           >  1:708999/1‑100 (MQ=255)
                                                                                        ggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc        <  1:41541/100‑1 (MQ=255)
                                                                                        ggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc        <  2:731544/100‑1 (MQ=255)
                                                                                        ggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc        >  2:2453677/1‑100 (MQ=255)
                                                                                        ggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc        <  2:1328634/100‑1 (MQ=255)
                                                                                        ggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc        >  1:942154/1‑100 (MQ=255)
                                                                                        ggAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGCAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTc        <  2:1937679/100‑1 (MQ=255)
                                                                                         gAAGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCg       <  1:2174673/100‑1 (MQ=255)
                                                                                           aGGGATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGAt     <  1:2074199/100‑1 (MQ=255)
                                                                                            gggATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATa    >  2:2642546/1‑100 (MQ=255)
                                                                                             ggATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATAc   >  1:764448/1‑100 (MQ=255)
                                                                                             ggATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATAc   >  1:2606206/1‑100 (MQ=255)
                                                                                             ggATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATAc   >  2:411352/1‑100 (MQ=255)
                                                                                             ggATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATAc   <  1:659740/100‑1 (MQ=255)
                                                                                              gATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACt  >  2:813520/1‑100 (MQ=255)
                                                                                              gATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACt  >  1:1637173/1‑100 (MQ=255)
                                                                                              gATACCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTAGTCAGTCTCTCGATACt  >  1:1194943/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TCCTGCCGCCTGGCGCGGTGATGACCGGTAATCACGAGATGTACAAAGTCTTTACGCCTTTTAAGAATGCCTGGCTGAAACGGCTGCGGGAAGGGATGCCGGAGTGCGTCGCTGCGCCAAAAGTTCGTAGTAGCGGATCGATAGAGCCCTCGCCATCCATTACGCTGAATTATCCTCGTCAGTCTTTCGATACT  >  NC_000913/739907‑740100

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