Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,303,518 A→C E24A (GAA→GCA)  yhbW → putative luciferase‑like monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,303,5180AC98.6% 496.8 / ‑2.1 143E24A (GAA→GCA) yhbWputative luciferase‑like monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (74/67);  total (74/69)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

Only 100 of 145 total aligned reads displayed.

TTAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTG  >  NC_000913/3303423‑3303613
                                                                                               |                                                                                               
ttAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAgcagca                                                                                             <  1:2570792/100‑1 (MQ=255)
ttAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAgcagca                                                                                             >  2:668353/1‑100 (MQ=255)
 tAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCAt                                                                                            >  2:2291633/1‑100 (MQ=255)
   aTCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATtc                                                                                          >  1:266537/1‑100 (MQ=255)
   aTCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATtc                                                                                          <  1:3585552/100‑1 (MQ=255)
   aTCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATtc                                                                                          >  2:1504711/1‑100 (MQ=255)
   aTCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATtc                                                                                          >  1:506838/1‑100 (MQ=255)
      ggTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGGCCCCATTCCCGAAGGTTCTTCAGCGCGATCAGCATTCTcc                                                                                       >  1:3602141/1‑100 (MQ=255)
      ggTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTcc                                                                                       >  1:1521680/1‑100 (MQ=255)
      ggTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTcc                                                                                       <  1:857597/100‑1 (MQ=255)
       gTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTccc                                                                                      >  2:2726929/1‑100 (MQ=255)
        tAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCa                                                                                     >  2:900843/1‑100 (MQ=255)
         aaCAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCAc                                                                                    <  2:1277578/100‑1 (MQ=255)
           cAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCActc                                                                                  >  2:1394647/1‑100 (MQ=255)
           cAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCActc                                                                                  <  2:3207005/100‑1 (MQ=255)
            aaaGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCActct                                                                                 >  2:2731451/1‑100 (MQ=255)
                tGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCg                                                                             >  2:377080/1‑100 (MQ=255)
                  aGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGAt                                                                           <  1:2773502/100‑1 (MQ=255)
                   gCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGAtc                                                                          >  1:1039528/1‑100 (MQ=255)
                   gCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGAtc                                                                          >  1:907834/1‑100 (MQ=255)
                   gCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGAtc                                                                          <  2:3222932/100‑1 (MQ=255)
                   gCTGTTATGACCGACAACACCATTGCGATTTCGCTACTCGGTATGGCCCCCGTTCCCGACGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCGCGAtc                                                                          <  2:118352/100‑1 (MQ=255)
                    cTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGAtct                                                                         <  1:1205497/100‑1 (MQ=255)
                      gTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCg                                                                       >  1:1041289/1‑100 (MQ=255)
                      gTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCg                                                                       >  1:2184432/1‑100 (MQ=255)
                        tATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGcc                                                                     <  2:1435293/100‑1 (MQ=255)
                        tATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGcc                                                                     >  2:883995/1‑100 (MQ=255)
                          tGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCg                                                                   >  1:2021174/1‑100 (MQ=255)
                          tGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATGCTCCCACTCTCTCGATCTCGCCCg                                                                   >  2:2775169/1‑100 (MQ=255)
                           gACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGt                                                                  >  2:1128150/1‑100 (MQ=255)
                             cTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCt                                                                >  2:695816/1‑100 (MQ=255)
                              tGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTg                                                               <  1:118606/100‑1 (MQ=255)
                              tGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTg                                                               <  2:3921431/100‑1 (MQ=255)
                              tGATAAAACCATTGCGTTTTCGCTACTCGATATGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTg                                                               >  1:802991/1‑100 (MQ=255)
                               gATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTgg                                                              >  1:3416098/1‑100 (MQ=255)
                               gATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCCCGAGCAGACTTCTCCCACTCTCTCGACCCCGCCCGTccgt                                                              >  2:3048709/1‑97 (MQ=255)
                                aTAAAATCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGCAGGTGCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCTTATGGc                                                             <  2:3760879/100‑1 (MQ=255)
                                aTAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGc                                                             <  1:903026/100‑1 (MQ=255)
                                 tAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGCAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCCGGCt                                                            >  2:2494700/1‑100 (MQ=255)
                                 tAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCt                                                            <  1:1690150/100‑1 (MQ=255)
                                  aaaaCCATTGCGTTTGCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTg                                                           <  1:3402766/100‑1 (MQ=255)
                                   aaaCCATTGCGTTTTCGCTACTCAATCTGGCCCCCATTCCCGCAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCCGGCtca                                                          >  2:730946/1‑98 (MQ=255)
                                    aaCCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGaa                                                         <  2:3859276/100‑1 (MQ=255)
                                        aTTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAgc                                                     >  2:1281940/1‑100 (MQ=255)
                                        aTTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAgc                                                     >  1:3688065/1‑100 (MQ=255)
                                          tGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAgcgc                                                   >  1:277913/1‑100 (MQ=255)
                                            cGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCgg                                                 >  2:1428299/1‑100 (MQ=255)
                                            cGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCgg                                                 >  2:1928700/1‑100 (MQ=255)
                                            cGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCgg                                                 <  1:1928407/100‑1 (MQ=255)
                                             gTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGc                                                >  1:2364058/1‑100 (MQ=255)
                                             gTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGc                                                <  2:3843021/100‑1 (MQ=255)
                                                  cGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTatca                                           <  2:1800681/100‑1 (MQ=255)
                                                  cGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTatca                                           >  1:2040054/1‑100 (MQ=255)
                                                    cTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTatcatc                                         >  2:1482124/1‑100 (MQ=255)
                                                      aCTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGc                                       >  2:348009/1‑100 (MQ=255)
                                                      aCTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCCCGAGCAGCATCCTCCCACTCTCTCGACCCCGCCCGTCTCGCTGAAAAGCGCGGCTATCAgcgc                                       >  2:1422440/1‑100 (MQ=255)
                                                        tCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTa                                     <  2:1030494/100‑1 (MQ=255)
                                                         cGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTAc                                    >  1:137780/1‑100 (MQ=255)
                                                         cGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTAc                                    >  1:2249998/1‑100 (MQ=255)
                                                          gATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCGTCGCTAct                                   <  2:3836266/100‑1 (MQ=255)
                                                           aTCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTActg                                  >  2:2024796/1‑100 (MQ=255)
                                                              tGGCCCCCATTCCCGAAGGTTCTTCAGCCCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTActggct                               <  2:2863420/100‑1 (MQ=255)
                                                               ggCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTActggctg                              <  1:3947622/100‑1 (MQ=255)
                                                                gCCCCCATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTActggctgg                             >  1:209700/1‑100 (MQ=255)
                                                                 cccccATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACtggctggc                            >  1:3915305/1‑100 (MQ=255)
                                                                 cccccATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACtggctggc                            >  1:582226/1‑100 (MQ=255)
                                                                 cccccATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACtggctggc                            <  1:1786186/100‑1 (MQ=255)
                                                                  ccccATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCa                           <  1:3376702/100‑1 (MQ=255)
                                                                  ccccATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCa                           <  1:1643342/100‑1 (MQ=255)
                                                                    ccATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGa                         >  1:682965/1‑100 (MQ=255)
                                                                    ccATTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGa                         <  2:78824/100‑1 (MQ=255)
                                                                      aTTCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAac                       <  1:1110078/100‑1 (MQ=255)
                                                                       ttCCCGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAaca                      >  1:498791/1‑100 (MQ=255)
                                                                         cccGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAAcacc                    >  1:1816701/1‑100 (MQ=255)
                                                                         cccGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAAcacc                    >  1:3990484/1‑100 (MQ=255)
                                                                          ccGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAAcacca                   >  1:3579093/1‑100 (MQ=255)
                                                                           cGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAAcaccac                  <  2:537158/100‑1 (MQ=255)
                                                                           cGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAAcaccac                  <  2:2327717/100‑1 (MQ=255)
                                                                           cGAAGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAAcaccac                  <  1:3331744/100‑1 (MQ=255)
                                                                              aGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAat               >  1:2926331/1‑100 (MQ=255)
                                                                              aGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAat               >  2:3163803/1‑100 (MQ=255)
                                                                              aGGTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAat               <  2:135478/100‑1 (MQ=255)
                                                                               ggTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACATCACAata              <  2:3628752/100‑1 (MQ=255)
                                                                               ggTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAata              >  2:2310434/1‑100 (MQ=255)
                                                                               ggTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAata              <  1:110370/100‑1 (MQ=255)
                                                                                gTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAatat             <  2:2507772/100‑1 (MQ=255)
                                                                                gTTCTTCAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCCGGCTGAACAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAatan             <  1:3264497/100‑2 (MQ=255)
                                                                                 ttcttcAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATg            >  2:3099074/1‑100 (MQ=255)
                                                                                 ttcctcAGCGCGAGCAGCATTCTCCCCCTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATg            >  1:3144177/1‑100 (MQ=255)
                                                                                   cttcAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGAc          >  1:2859930/1‑100 (MQ=255)
                                                                                   cttcAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGAc          <  1:251395/100‑1 (MQ=255)
                                                                                   cttcAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGAc          <  1:1906400/100‑1 (MQ=255)
                                                                                   cttcAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGAc          >  2:2801909/1‑100 (MQ=255)
                                                                                    ttcAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACt         <  1:1066237/100‑1 (MQ=255)
                                                                                     tcAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTg        <  1:1629620/100‑1 (MQ=255)
                                                                                     tcAGCGCGAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTg        <  2:1803889/100‑1 (MQ=255)
                                                                                         cgcgAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCAt    <  2:1647347/100‑1 (MQ=255)
                                                                                         cgcgAGCAGCATAATCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCCCGGGTATCCTCGCTACTGACTGGCAGAACACCACAATATGACTGGCAt    <  2:1541591/100‑1 (MQ=255)
                                                                                          gcgAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCAtt   <  2:2933921/100‑1 (MQ=255)
                                                                                           cgAGCAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTg  <  1:301141/100‑1 (MQ=255)
                                                                                               |                                                                                               
TTAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTG  >  NC_000913/3303423‑3303613

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