Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,530,505 G→A 15.8% intergenic (+101/‑81) ydcD → / → yncI putative immunity protein for RhsE/pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,530,5050GA15.8% 365.3 / 50.3 145intergenic (+101/‑81)ydcD/yncIputative immunity protein for RhsE/pseudogene
Reads supporting (aligned to +/- strand):  ref base G (57/65);  new base A (9/14);  total (66/79)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 145 total aligned reads displayed.

ACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTG  >  NC_000913/1530408‑1530603
                                                                                                 |                                                                                                  
acaacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAt                                                                                                  <  1:3747282/100‑1 (MQ=255)
 caacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtt                                                                                                 >  2:2816216/1‑100 (MQ=255)
 caacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtt                                                                                                 <  1:2691336/100‑1 (MQ=255)
  aacaCCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAttt                                                                                                <  1:1464405/100‑1 (MQ=255)
   acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                                               <  1:2539721/100‑1 (MQ=255)
   acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                                               >  1:3731149/1‑100 (MQ=255)
   acaacaGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                                               <  1:1200495/100‑1 (MQ=255)
     acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                                             <  2:308728/100‑1 (MQ=38)
     acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                                             <  2:774146/100‑1 (MQ=38)
     acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                                             <  2:3704836/100‑1 (MQ=38)
      ccAGAAAGAGAAAAATATAAAATACATGCAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                                            <  2:3081755/100‑1 (MQ=37)
      ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCAt                                                                                            >  2:656780/1‑100 (MQ=255)
      ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCAt                                                                                            >  2:3407229/1‑100 (MQ=255)
      ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                                            <  1:427668/100‑1 (MQ=38)
      ccAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAt                                                                                            <  2:1176875/100‑1 (MQ=38)
       cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                                           <  1:2364434/100‑2 (MQ=37)
       cAGAAAGAGAAAAATATAAAATACATACAGAAAGACAGGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATc                                                                                           >  1:2037538/1‑100 (MQ=255)
         gAAAGAGAAAAATATAAAATACATACAGAAATACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGAGGTAAATCACGATTTTCATCtt                                                                                         <  1:795597/100‑1 (MQ=255)
         gAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCtt                                                                                         <  2:1083878/100‑1 (MQ=255)
          aaaGAGAAGATTATAAAACACATATACAAAGACATGAGAGCAAGATACAGATCAGGTTAAGCTCAACAAAATAGTAATTTACAACTCAATTTTCAtg                                                                                           <  1:1792911/97‑2 (MQ=255)
          aaaGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGTACGGGGCAGGGCAAGATTATAAAAATCGTGATTTACATCACAATTTTCAtgat                                                                                         <  1:1158111/99‑4 (MQ=25)
           aaGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTg                                                                                       <  2:1496519/100‑1 (MQ=255)
            agagAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGc                                                                                      <  2:779907/100‑1 (MQ=255)
             gagaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCt                                                                                     >  1:2096847/1‑100 (MQ=255)
              agaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCtt                                                                                    <  2:2763712/100‑1 (MQ=255)
               gaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGAAGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATATTTGCttt                                                                                   >  1:2278334/1‑100 (MQ=255)
                aaaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCTCGATTTTCATCTTTGCTTTa                                                                                  <  1:2322871/100‑1 (MQ=255)
                  aaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAAc                                                                                <  1:460506/100‑1 (MQ=255)
                   aaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGTGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACg                                                                               <  1:1363601/100‑1 (MQ=255)
                   aaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACg                                                                               <  1:820567/100‑1 (MQ=255)
                   aaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACg                                                                               <  1:800251/100‑1 (MQ=255)
                    atatAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGc                                                                              <  2:2962896/100‑1 (MQ=255)
                          aaTACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAg                                                                        <  1:328133/100‑1 (MQ=255)
                           atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgn                                                                       <  2:1461873/100‑2 (MQ=255)
                           atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                       <  2:2635560/100‑1 (MQ=255)
                           atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                       >  2:3271242/1‑100 (MQ=255)
                           atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                                       <  2:145528/100‑1 (MQ=255)
                              catacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGa                                                                    <  2:158557/100‑1 (MQ=255)
                              catacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGa                                                                    >  2:2519413/1‑100 (MQ=255)
                               atacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAAGTGTGATGTAAATCACGATTTTCATCTTTGCTTTAAAGCCTAAAg                                                                        >  1:37431/1‑95 (MQ=255)
                                tacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                                  <  2:1311441/100‑1 (MQ=255)
                                tacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                                  >  2:1319559/1‑100 (MQ=255)
                                tacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                                  >  2:742730/1‑100 (MQ=255)
                                tacaGAAAGACATGAGAGCAAGGGACGGGGAAGGGCAAGATCATGAAAATTATGATGTAAATCACGATTTTCATCTTTGCTTTAACGCATACAGGTg                                                                     >  1:123532/1‑97 (MQ=255)
                                  caGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaa                                                                >  2:3280529/1‑100 (MQ=255)
                                   aGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaaa                                                               >  1:1378671/1‑100 (MQ=255)
                                    gAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCaaaa                                                              <  2:3812714/100‑1 (MQ=255)
                                       aGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAAttt                                                           <  1:3040660/100‑1 (MQ=255)
                                        gACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTa                                                          <  2:3788447/100‑1 (MQ=255)
                                         aCATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTaa                                                         <  2:2290261/100‑1 (MQ=255)
                                          cATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                                        <  1:1075710/100‑1 (MQ=255)
                                          cATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                                        <  2:1746661/100‑1 (MQ=255)
                                          cATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                                        <  2:1990206/100‑1 (MQ=255)
                                          cATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                                        >  2:3595869/1‑100 (MQ=255)
                                           aTGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAcc                                                       <  1:3184157/100‑1 (MQ=255)
                                            tGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCt                                                      >  1:3362433/1‑100 (MQ=255)
                                            tGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCt                                                      <  1:30149/100‑1 (MQ=255)
                                            tGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCt                                                      >  2:2533457/1‑100 (MQ=255)
                                              agagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCttt                                                    >  1:3642527/1‑100 (MQ=255)
                                               gagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTg                                                   <  2:3380170/100‑1 (MQ=255)
                                               gagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTg                                                   >  2:392389/1‑100 (MQ=255)
                                               gagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTg                                                   >  2:2181803/1‑100 (MQ=255)
                                                agCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTgg                                                  <  2:105579/100‑1 (MQ=255)
                                                 gCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGa                                                 >  2:1731911/1‑100 (MQ=255)
                                                  cAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaa                                                <  2:2640824/100‑1 (MQ=255)
                                                  cAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaa                                                >  1:1925510/1‑100 (MQ=255)
                                                   aaGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAAt                                               >  1:907942/1‑100 (MQ=255)
                                                   aaGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAAt                                               <  1:242107/100‑1 (MQ=255)
                                                      ggACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAGTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAAc                                            >  1:1108671/1‑100 (MQ=255)
                                                       gACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACt                                           >  1:1598609/1‑100 (MQ=255)
                                                            ggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaaa                                      >  2:2170580/1‑100 (MQ=255)
                                                                gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATa                                  <  2:1549271/100‑1 (MQ=255)
                                                                 ggcaAGATGATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGATTTAACGCCTACCGCTAATCAAAATTTAACCTTTGGAATAACTAAAAAGATaa                                 >  1:278311/1‑100 (MQ=255)
                                                                 aaagcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaa                                 >  1:3027808/10‑100 (MQ=21)
                                                                 aaagcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaa                                 >  1:1884626/10‑100 (MQ=21)
                                                                  gcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa                                >  2:1181644/1‑100 (MQ=255)
                                                                   ngcaaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaa                               >  2:1944319/8‑100 (MQ=25)
                                                                   caAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaa                               <  1:3651172/100‑1 (MQ=255)
                                                                     aGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaa                              <  1:3438730/99‑1 (MQ=255)
                                                                      aaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg                            >  2:3616846/5‑100 (MQ=25)
                                                                         gaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc                         >  1:706802/2‑100 (MQ=37)
                                                                         gaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc                         >  2:1887526/2‑100 (MQ=37)
                                                                          aTGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAgg                           >  1:2926395/1‑97 (MQ=255)
                                                                            gAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGcc                      >  1:1662958/1‑100 (MQ=255)
                                                                             aaaaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCa                     <  1:309014/100‑1 (MQ=255)
                                                                             aaaaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCa                     <  2:2460884/100‑1 (MQ=255)
                                                                             aaaaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCa                     <  1:1463338/100‑1 (MQ=255)
                                                                               aaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAgg                   <  2:937305/100‑1 (MQ=255)
                                                                                   gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAg               <  2:3562551/100‑1 (MQ=255)
                                                                                   gtgATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAg               <  2:2453678/100‑1 (MQ=255)
                                                                                     gATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTa             >  1:282269/1‑100 (MQ=255)
                                                                                      aTGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACATTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTAt            >  1:2447229/1‑100 (MQ=255)
                                                                                         tAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCa         >  1:2990837/1‑100 (MQ=255)
                                                                                         tAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCa         >  2:2017882/1‑100 (MQ=255)
                                                                                          aaaTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCaa        >  2:1086357/1‑100 (MQ=255)
                                                                                            aTCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAg      >  2:90436/1‑100 (MQ=255)
                                                                                             tCACGATGTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACATTTGTAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGt     >  1:156990/1‑100 (MQ=255)
                                                                                              cACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGtt    <  1:1817779/100‑1 (MQ=255)
                                                                                               aCGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGttt   <  1:683987/100‑1 (MQ=255)
                                                                                                cGATTTTCATCTTTGATTTAACGCCTACAGCTGATAAAAATTTAACCTTTGGAATAACCAAAAAGATAAAAAAGGAAGCCAGTTGAGCACTCAATGtt    >  1:3528044/1‑98 (MQ=255)
                                                                                                 |                                                                                                  
ACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGGTGAGTATTCAAAGTTTG  >  NC_000913/1530408‑1530603

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: