Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,530,490 T→C 16.3% intergenic (+86/‑96) 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,4900TC16.3% 320.1 / 52.2 137intergenic (+86/‑96)ydcD/yncIputative immunity protein for RhsE/pseudogene
Reads supporting (aligned to +/- strand):  ref base T (50/64);  new base C (9/13);  total (59/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 141 total aligned reads displayed.

TAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGG  >  NC_000913/1530392‑1530586
                                                                                                  |                                                                                                
tAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATtg                                                                                                 <  2:1019843/100‑1 (MQ=255)
tAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATtg                                                                                                 <  2:101727/100‑1 (MQ=255)
  aaTGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATtgtg                                                                                               >  2:191383/1‑100 (MQ=255)
  aaTGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATtgtg                                                                                               <  2:160283/100‑1 (MQ=255)
   aTGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGa                                                                                              <  1:2140056/100‑1 (MQ=255)
    tGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGAt                                                                                             <  2:894685/100‑1 (MQ=255)
     ggTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATg                                                                                            <  2:2613304/100‑1 (MQ=255)
     ggTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATg                                                                                            <  2:2102566/100‑1 (MQ=255)
     ggTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGACAGGGCAAGATCATGAAAATTGTGAt                                                                                             >  1:3367642/1‑99 (MQ=255)
       tCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGAGATGTa                                                                                          >  2:1736193/1‑100 (MQ=255)
        cATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTaa                                                                                         >  2:2362165/1‑100 (MQ=255)
          tAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGATCAATGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAAt                                                                                       <  1:1921056/100‑1 (MQ=255)
          tAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAAt                                                                                       <  2:3041100/100‑1 (MQ=255)
          tAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAAt                                                                                       >  2:2023785/1‑100 (MQ=255)
                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:3731149/1‑100 (MQ=255)
                   acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                              <  1:2539721/100‑1 (MQ=255)
                   acacCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAAtttt                                                                              <  2:591824/100‑1 (MQ=38)
                   acaacaGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAtttt                                                                              <  1:1200495/100‑1 (MQ=255)
                     acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                            <  2:774146/100‑1 (MQ=38)
                     acCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCa                                                                            <  2:3600043/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)
                       cAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATCGTGATTTACATCACAATTTTCAtg                                                                          <  2:2024741/100‑2 (MQ=37)
                        aGAAAGAGTAAAATATAAAATACATACAGAAAGACATGAGAGCAAGAGACGGGGCAGGGCAAGATCATGAAATTCGTGATTTACATCACAATTTTCAt                                                                           <  2:1377457/98‑1 (MQ=37)
                        aGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCt                                                                         >  2:3104846/1‑100 (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)
                             gagaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCt                                                                    >  1:2096847/1‑100 (MQ=255)
                              agaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCtt                                                                   <  2:2763712/100‑1 (MQ=255)
                               gaAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCttt                                                                  >  2:2376982/1‑100 (MQ=255)
                                aaaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCTCGATTTTCATCTTTGCTTTa                                                                 <  1:2322871/100‑1 (MQ=255)
                                aaaaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTa                                                                 >  1:3066665/1‑100 (MQ=255)
                                  aaaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAAc                                                               <  1:460506/100‑1 (MQ=255)
                                   aaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGTGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACg                                                              <  1:1363601/100‑1 (MQ=255)
                                   aaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACg                                                              <  1:800251/100‑1 (MQ=255)
                                   aaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACg                                                              <  1:819521/100‑1 (MQ=255)
                                   aaTATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACg                                                              <  1:820567/100‑1 (MQ=255)
                                          aaTACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAg                                                       <  1:328133/100‑1 (MQ=255)
                                           atacatacAGAAAGACATGAGAGCAAGGGCCGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                      <  2:759769/100‑1 (MQ=255)
                                           atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                      >  2:3271242/1‑100 (MQ=255)
                                           atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                      <  2:2635560/100‑1 (MQ=255)
                                           atacatacAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAgg                                                      <  2:145528/100‑1 (MQ=255)
                                              catacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGa                                                   >  2:2519413/1‑100 (MQ=255)
                                              catacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGa                                                   <  2:158557/100‑1 (MQ=255)
                                               atacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAAGTGTGATGTAAATCACGATTTTCATCTTTGCTTTAAAGCCTAAAg                                                       >  1:37431/1‑95 (MQ=255)
                                                tacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                 >  2:1319559/1‑100 (MQ=255)
                                                tacaGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATc                                                 >  1:1475638/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)
                                                       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:1990206/100‑1 (MQ=255)
                                                          cATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                       <  2:1746661/100‑1 (MQ=255)
                                                          cATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAAc                                       >  2:3595869/1‑100 (MQ=255)
                                                            tGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCt                                     <  1:30149/100‑1 (MQ=255)
                                                              agagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCttt                                   >  1:3642527/1‑100 (MQ=255)
                                                               gagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTg                                  >  2:392389/1‑100 (MQ=255)
                                                               gagCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTg                                  <  2:3380170/100‑1 (MQ=255)
                                                                 gCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGa                                >  2:1731911/1‑100 (MQ=255)
                                                                  cAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaa                               >  1:1925510/1‑100 (MQ=255)
                                                                  cAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGaa                               <  2:2640824/100‑1 (MQ=255)
                                                                   aaGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAAt                              <  1:242107/100‑1 (MQ=255)
                                                                   aaGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAAt                              >  1:1836759/1‑100 (MQ=255)
                                                                   aaGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAAt                              <  1:3385376/100‑1 (MQ=255)
                                                                    aGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATa                             <  1:1447479/100‑1 (MQ=255)
                                                                      ggACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGAGTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAAc                           >  1:1108671/1‑100 (MQ=255)
                                                                       gACGGGGCAGGGCAAGATCATGCAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATa                             >  1:2904229/1‑97 (MQ=255)
                                                                           gggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaa                      >  2:2500822/1‑100 (MQ=255)
                                                                            ggcagggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTaaaaa                     >  2:2170580/1‑100 (MQ=255)
                                                                                nggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATa                 >  1:3007427/2‑100 (MQ=255)
                                                                                gggcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATa                 <  2:1088404/100‑1 (MQ=255)
                                                                                 ggcaAGATGATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGATTTAACGCCTACCGCTAATCAAAATTTAACCTTTGGAATAACTAAAAAGATaa                >  1:278311/1‑100 (MQ=255)
                                                                                  gcaAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaa               >  2:1181644/1‑100 (MQ=255)
                                                                                   caAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaa              <  1:3651172/100‑1 (MQ=255)
                                                                                     aGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATaaaaa             <  1:3438730/99‑1 (MQ=255)
                                                                                      gATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg           <  1:742722/100‑1 (MQ=255)
                                                                                      aaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg           >  1:3220820/5‑100 (MQ=25)
                                                                                      aaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg           >  2:3616846/5‑100 (MQ=25)
                                                                                      aaagaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAg           >  2:118858/5‑100 (MQ=37)
                                                                                         gaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc        >  1:706802/2‑100 (MQ=37)
                                                                                         gaTGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc        >  2:1887526/2‑100 (MQ=37)
                                                                                         cATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGAc        <  2:1481337/100‑1 (MQ=255)
                                                                                          aTGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACg       <  2:1066364/100‑1 (MQ=255)
                                                                                           tGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGc      <  1:2362154/100‑1 (MQ=255)
                                                                                           tGAAAATCGTGATTTACATCACAATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGc      >  2:3800986/1‑100 (MQ=38)
                                                                                            gAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGcc     >  1:1662958/1‑100 (MQ=255)
                                                                                             aaaaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCa    <  1:309014/100‑1 (MQ=255)
                                                                                               aaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAgg  <  2:937305/100‑1 (MQ=255)
                                                                                               aaTTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCCTTGGAATAACTAAAAAGATAAAAAAGGACGCCAgg  <  1:2102369/100‑1 (MQ=255)
                                                                                                  |                                                                                                
TAAATGGTCATAATGAACAACACCAGAAAGAGAAAAATATAAAATACATACAGAAAGACATGAGAGCAAGGGACGGGGCAGGGCAAGATCATGAAAATTGTGATGTAAATCACGATTTTCATCTTTGCTTTAACGCCTACAGGTGATCAAAATTTAACCTTTGGAATAACTAAAAAGATAAAAAAGGACGCCAGG  >  NC_000913/1530392‑1530586

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