Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,515,696 T→C Y73C (TAC→TGC)  dam ← DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,515,6960TC98.6% 243.0 / NA 69Y73C (TAC→TGC) damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (24/44);  total (24/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATT  >  NC_000913/3515559‑3515828
                                                                                                                                         |                                                                                                                                    
aaCCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGAACAAACAGCTCGCGTGCGGCCTGTACGCACt                                                                                                                                   <  2:688094/140‑1 (MQ=255)
 aCCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACt                                                                                                                                   <  2:87336/140‑1 (MQ=255)
 aCCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACt                                                                                                                                   <  1:873233/140‑1 (MQ=255)
   cGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACtca                                                                                                                                 <  2:696057/140‑1 (MQ=255)
      aGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTGGCGTGCGGCCTGTACGCACtcatca                                                                                                                              >  1:326536/1‑140 (MQ=255)
       gCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAg                                                                                                                             >  1:753245/1‑140 (MQ=255)
       gCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAg                                                                                                                             >  2:1005656/1‑140 (MQ=255)
       gCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAAAAAATAGCTCGCGTGCGGCCTGTACGCACTCATCAg                                                                                                                             >  2:378639/1‑140 (MQ=255)
           ttCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACg                                                                                                                         <  1:738712/140‑1 (MQ=255)
           ttCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGCTTCGGGAACAAAAACCCCGCGTGCGGCCTGCACGCACTCCTCAGTACg                                                                                                                         >  2:499189/1‑140 (MQ=255)
                atatAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCt                                                                                                                    >  1:654120/1‑140 (MQ=255)
                 tataAAAACAGTACCGCCCGACGGAACGGATCCTGTCTTCTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCtt                                                                                                                   <  2:353308/140‑1 (MQ=255)
                   taAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTca                                                                                                                 <  2:150752/140‑1 (MQ=255)
                   taAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTca                                                                                                                 <  1:953276/140‑1 (MQ=255)
                   taAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTca                                                                                                                 <  2:977820/140‑1 (MQ=255)
                     aaaaCAGTACCGCCCGACGGAAGGGATCCTGGCTTCTGTCGAACTCTTTGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTcaca                                                                                                               <  2:388986/140‑1 (MQ=255)
                     aaaaCAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTcaca                                                                                                               <  2:983904/140‑1 (MQ=255)
                     aaaaCAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTcaca                                                                                                               <  1:166507/140‑1 (MQ=255)
                       aaCAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAAt                                                                                                             <  1:12479/140‑1 (MQ=255)
                       aaCAGTACCGCCCGACGGAACGGATCCTGGCATTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAAt                                                                                                             <  2:450261/140‑1 (MQ=255)
                            tACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCGTGTACGCGCTCTTCAGTACGCCTCTTCACAATGTtat                                                                                                        >  1:313053/1‑140 (MQ=255)
                            tACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTtat                                                                                                        <  2:533779/140‑1 (MQ=255)
                            tACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTtat                                                                                                        >  2:902896/1‑140 (MQ=255)
                            tACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTtat                                                                                                        <  1:7861/140‑1 (MQ=255)
                            tACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTtat                                                                                                        >  1:944909/1‑140 (MQ=255)
                              ccGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATag                                                                                                      <  2:654120/140‑1 (MQ=255)
                                 cccGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTAGGCACTCATCAGTACGCATCTTCACAATGTTATAgaga                                                                                                   >  1:455314/1‑140 (MQ=255)
                                 cccGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTca                                                                                                                 >  1:675104/1‑126 (MQ=255)
                                 cccGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTca                                                                                                                 <  2:675104/126‑1 (MQ=255)
                                 cccGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGGAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAgaga                                                                                                   >  2:87475/1‑140 (MQ=255)
                                 cccGACGGAACGGATCCTGGCTATTGTTGTACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAgaga                                                                                                   <  1:499189/140‑1 (MQ=255)
                                         aacggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAg                                                                                           >  2:836886/1‑140 (MQ=255)
                                           cggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGt                                                                                         <  2:944909/140‑1 (MQ=255)
                                           cggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGt                                                                                         <  2:753245/140‑1 (MQ=255)
                                           cggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGt                                                                                         <  1:4463/140‑1 (MQ=255)
                                            ggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTc                                                                                        <  1:87475/140‑1 (MQ=255)
                                            ggaTCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTc                                                                                        <  2:756019/140‑1 (MQ=255)
                                                    gCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCGTCTGCACAATGTTATAGAGACAGATCAGGTCGCTATTGa                                                                                >  2:390538/1‑140 (MQ=255)
                                                     cTTTTGTTGAACTCTACGCGGAACTCATAGTAAACCTGGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGat                                                                               <  1:430785/140‑1 (MQ=255)
                                                      ttttGTTGAACTCTTCGCGGAACTGGTAGTAAACCTCTGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGata                                                                              <  1:359588/140‑1 (MQ=255)
                                                      ttttGTTGAACTCTTCGCGGAACTGATAGTACACCTCGGTGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGata                                                                              <  2:299732/140‑1 (MQ=255)
                                                        ttgttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCCCAATGTTATAGAGAGTGATCAGGTCGCTATTGATATc                                                                            >  1:421121/1‑140 (MQ=255)
                                                          gttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTTATATCgg                                                                          <  1:855945/140‑1 (MQ=255)
                                                          gttgAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCgg                                                                          <  1:833008/140‑1 (MQ=255)
                                                                   ttCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGatat                                                                 >  1:88270/1‑140 (MQ=255)
                                                                    tCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGAtata                                                                <  1:902896/140‑1 (MQ=255)
                                                                        ggAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACga                                                            >  1:760654/1‑140 (MQ=255)
                                                                          aaCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACgaga                                                          <  2:760654/140‑1 (MQ=255)
                                                                          aaCTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACgaga                                                          <  1:508958/140‑1 (MQ=255)
                                                                          aaCTGATAGTAAACCGCGGCGCGATATGTGTCGGGAATAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACgaga                                                          <  1:476151/140‑1 (MQ=255)
                                                                               aTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                     <  1:208100/140‑1 (MQ=255)
                                                                               aTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                     <  2:804962/140‑1 (MQ=255)
                                                                                          cGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGATAAGTCGGTGTTgaga                                          >  1:430147/1‑140 (MQ=255)
                                                                                           ggCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                     >  2:869043/1‑128 (MQ=255)
                                                                                           ggCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGt                                                     <  1:869043/128‑1 (MQ=255)
                                                                                           ggCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaa                                         >  1:713716/1‑140 (MQ=255)
                                                                                           ggCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaa                                         >  1:644125/1‑140 (MQ=255)
                                                                                           ggCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGaa                                         >  2:143259/1‑140 (MQ=255)
                                                                                             cgcAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAac                                       >  2:70508/1‑140 (MQ=255)
                                                                                                      tttCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCg                              <  2:539588/140‑1 (MQ=255)
                                                                                                      tttCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCg                              >  1:865652/1‑140 (MQ=255)
                                                                                                       ttCGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCgg                             <  1:858921/140‑1 (MQ=255)
                                                                                                         cGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCa                           <  1:70508/140‑1 (MQ=255)
                                                                                                         cGGGAACAAACAGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCa                           <  2:13017/140‑1 (MQ=255)
                                                                                                               caaacaGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCgg                             <  2:786854/132‑1 (MQ=255)
                                                                                                               caaacaGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCgg                             >  1:786854/1‑132 (MQ=255)
                                                                                                                aaacaGCTCGCGTGCGGCCTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACaa                    <  2:98287/140‑1 (MQ=255)
                                                                                                                    aGCTCGCGTGCGGCCTTTACGCCCTCATCAGGATGAATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAgg                <  2:455314/140‑1 (MQ=255)
                                                                                                                        cgcgTGCGGCCTGTACGTACTCATCAGTACGGATCTTCAGAATGTTATTGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCa            <  2:378210/140‑1 (MQ=255)
                                                                                                                           gTGCGGCCTGTACGCACTCATCCGTACGCATCTTCTCAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAAcc         <  2:326536/140‑1 (MQ=255)
                                                                                                                                  cTGTACGCACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACAtt  <  1:143259/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
AACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCATCTTCACAATGTTATAGAGACTGATCAGGTCGCTATTGATATCGGCAAGGATATAACGAGAAAAGTCGGTGTTGAGAAACACCGACCCGGCACCTACAAAAGGCTCAACCAGACATT  >  NC_000913/3515559‑3515828

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