Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 833,889 C→T L274L (CTA→TTA)  dinG → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913833,8890CT96.8% 313.3 / 0.1 93L274L (CTA→TTA) dinGATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (42/48);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

ATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAACGTTTGAATG  >  NC_000913/833794‑833985
                                                                                               |                                                                                                
aTTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACa                                                                                              <  2:539383/100‑1 (MQ=255)
 tttACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAg                                                                                             >  1:1384230/1‑100 (MQ=255)
  ttACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGc                                                                                            <  2:50770/100‑1 (MQ=255)
    aCTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTg                                                                                          <  1:1377661/100‑1 (MQ=255)
       gctggtgctggACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGAc                                                                                       >  2:1957100/1‑100 (MQ=255)
       gctggtgctggACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGAc                                                                                       <  1:2267970/100‑1 (MQ=255)
       gctggtgctggACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGAc                                                                                       >  2:478705/1‑100 (MQ=255)
            tgctggACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGtt                                                                                  >  2:100620/1‑100 (MQ=255)
             gctggACGTAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTc                                                                                 <  2:694632/100‑1 (MQ=255)
             gctggACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTc                                                                                 >  2:444588/1‑100 (MQ=255)
             gctggACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTc                                                                                 <  2:1384227/100‑1 (MQ=255)
             gctggACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTc                                                                                 <  1:841386/100‑1 (MQ=255)
              ctggACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCa                                                                                <  1:2000913/100‑1 (MQ=255)
              ctggACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCa                                                                                <  1:2000999/100‑1 (MQ=255)
               tggACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCAc                                                                               >  2:1657015/1‑100 (MQ=255)
               tggACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCAc                                                                               >  1:561821/1‑100 (MQ=255)
                ggACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACg                                                                              >  2:941672/1‑100 (MQ=255)
                ggACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACg                                                                              >  2:1928656/1‑100 (MQ=255)
                ggACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACg                                                                              >  2:1929363/1‑100 (MQ=255)
                    atctgCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAAc                                                                          <  2:118776/96‑1 (MQ=255)
                     aaGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACt                                                                         <  1:605384/100‑1 (MQ=255)
                     aaGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACt                                                                         >  1:1053531/1‑100 (MQ=255)
                     aaGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACt                                                                         >  2:1465120/1‑100 (MQ=255)
                     aaGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACt                                                                         >  1:404980/1‑100 (MQ=255)
                        gCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAAcagat                                                                      >  1:118778/1‑96 (MQ=255)
                        gCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGt                                                                      >  2:674387/1‑100 (MQ=255)
                        gCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGt                                                                      >  1:1698484/1‑100 (MQ=255)
                         ccATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGtn                                                                     <  1:2118147/100‑2 (MQ=255)
                         ccATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTc                                                                     <  1:2118034/100‑1 (MQ=255)
                         ccATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTc                                                                     >  2:226446/1‑100 (MQ=255)
                         ccATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTc                                                                     <  2:2251384/100‑1 (MQ=255)
                         ccATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTc                                                                     <  1:1477126/100‑1 (MQ=255)
                          cATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCg                                                                    <  2:89843/100‑1 (MQ=255)
                          cATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCg                                                                    <  1:1671734/100‑1 (MQ=255)
                          cATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCACGAAACTGGTCg                                                                    >  2:117901/1‑100 (MQ=255)
                               ccTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTAcc                                                               <  2:2080346/100‑1 (MQ=255)
                                 tGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTg                                                             <  1:477036/100‑1 (MQ=255)
                                 tGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTg                                                             <  2:294052/100‑1 (MQ=255)
                                  gCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGc                                                            >  2:1839552/1‑100 (MQ=255)
                                    cGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCAt                                                          >  1:1325588/1‑100 (MQ=255)
                                     ggATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGCCTTGTTCACGAAACTGGTCGCTACCTGCATg                                                         >  1:1369676/1‑100 (MQ=255)
                                      gctgtgGCGCGGGATGTGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGTTCGCTACCGGCATgg                                                        <  2:2182986/98‑1 (MQ=255)
                                       aTGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGa                                                       <  2:25082/100‑1 (MQ=255)
                                        tgtgGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAg                                                      >  1:644958/1‑100 (MQ=255)
                                         gtgGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAg                                                      <  2:644957/99‑1 (MQ=255)
                                         gtgGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGc                                                     <  2:1131126/100‑1 (MQ=255)
                                         gtgGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGc                                                     <  2:561820/100‑1 (MQ=255)
                                         gtgGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGc                                                     <  1:1957107/100‑1 (MQ=255)
                                         gtgGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGc                                                     <  2:2302992/100‑1 (MQ=255)
                                           ggCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAg                                                   >  1:2040798/1‑100 (MQ=255)
                                           ggCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAg                                                   >  2:2240450/1‑100 (MQ=255)
                                           ggCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAg                                                   >  1:1747001/1‑100 (MQ=255)
                                           ggCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAg                                                   <  1:226448/100‑1 (MQ=255)
                                            gcgcGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGt                                                  <  2:698862/100‑1 (MQ=255)
                                              gcgGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGttt                                                <  2:1683796/100‑1 (MQ=255)
                                                gggATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCg                                              >  1:597539/1‑100 (MQ=255)
                                                 ggATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGc                                             <  1:669654/100‑1 (MQ=255)
                                                   aTGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGccc                                           >  1:1219830/1‑100 (MQ=255)
                                                    tGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCg                                          <  2:372132/100‑1 (MQ=255)
                                                    tGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCg                                          <  2:486266/100‑1 (MQ=255)
                                                     gcgcTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGa                                         >  2:982020/1‑100 (MQ=255)
                                                      cgcTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGCAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCgca                                        >  1:1609825/1‑98 (MQ=255)
                                                       gcTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAg                                       <  2:450293/100‑1 (MQ=255)
                                                           gagaTGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCa                                   >  2:1885662/1‑100 (MQ=255)
                                                             gaTGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATc                                 >  1:535484/1‑100 (MQ=255)
                                                             gaTGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATc                                 >  2:975966/1‑100 (MQ=255)
                                                               tGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCAtgaa                               >  2:1776545/1‑97 (MQ=255)
                                                               tGAGCGCCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATccc                               <  1:2114829/100‑1 (MQ=255)
                                                                 ngcgcCGAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCAc                             >  2:1716345/2‑100 (MQ=255)
                                                                       gAAATCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAAACCATACCa                              >  2:2045843/1‑93 (MQ=255)
                                                                          aTCACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGa                    >  2:1685305/1‑100 (MQ=255)
                                                                            cACCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATc                  >  2:757858/1‑100 (MQ=255)
                                                                             aCCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATcc                 <  2:1369673/100‑1 (MQ=255)
                                                                             aCCGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATcc                 <  2:1325585/100‑1 (MQ=255)
                                                                              ccGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATccc                <  2:935854/100‑1 (MQ=255)
                                                                              ccGCGCCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATccc                <  1:478707/100‑1 (MQ=255)
                                                                                gcgcCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTg              >  1:793521/1‑100 (MQ=255)
                                                                                  gcCGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGaa            >  1:862447/1‑100 (MQ=255)
                                                                                  gcCGTGGTATCGGCTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGaa            >  2:920300/1‑100 (MQ=255)
                                                                                   ccGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAAc           <  1:444590/100‑1 (MQ=255)
                                                                                   ccGTGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCCGGCGATCCCTGAAc           <  2:1781122/100‑1 (MQ=255)
                                                                                      tGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAACGtt        <  2:404978/100‑1 (MQ=255)
                                                                                      tGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAACGtt        <  2:1219827/100‑1 (MQ=255)
                                                                                      tGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAACGtt        >  1:943939/1‑100 (MQ=255)
                                                                                      tGGTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAACGtt        <  2:869921/100‑1 (MQ=255)
                                                                                        gTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAACGTTTg      <  1:1465124/100‑1 (MQ=255)
                                                                                        gTATCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAACGTTTg      <  1:970390/100‑1 (MQ=255)
                                                                                          aTCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAACGTTTGaa    >  2:622496/1‑100 (MQ=255)
                                                                                          aTCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAACGTTTGaa    <  2:352076/100‑1 (MQ=255)
                                                                                          aTCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAACGTTTGaa    <  2:1559007/100‑1 (MQ=255)
                                                                                          aTCGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAACGTTTGaa    >  2:2164614/1‑100 (MQ=255)
                                                                                            cGGTTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAACGTTTGAATg  >  2:1551534/1‑100 (MQ=255)
                                                                                            cGGCTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAACGTTTGAATg  <  2:390229/100‑1 (MQ=255)
                                                                                               |                                                                                                
ATTTACTGCTGGTGCTGGACGAAGGCCATCACCTGCCGGATGTGGCGCGGGATGCGCTTGAGATGAGCGCCGAAATCACCGCGCCGTGGTATCGGCTACAGCTGGACTTGTTCACGAAACTGGTCGCTACCTGCATGGAGCAGTTTCGCCCGAAGACCATCCCACCGCTGGCGATCCCTGAACGTTTGAATG  >  NC_000913/833794‑833985

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