Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,478,094 A→G G311G (GGA→GGG b1410 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,478,0940AG99.6% 797.4 / NA 256G311G (GGA→GGGb1410hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (130/125);  minor base C (0/1);  total (130/126)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 271 total aligned reads displayed.

TTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCT  >  NC_000913/1478004‑1478183
                                                                                          |                                                                                         
tttATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGcc                                                                                         <  1:6316849/93‑1 (MQ=255)
tttATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGcc                                                                                         >  2:2957587/1‑93 (MQ=255)
  tataTGCTAACAAATCGCAAAAATTTGATTATCAGCATTAATACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCg                                                                                       <  1:3805357/93‑1 (MQ=255)
  tataTGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTAGTTTCTGGTGGGCCCg                                                                                       >  1:4213170/1‑93 (MQ=255)
  tataTGCTAACAAATCGCAAAAATCTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCg                                                                                       >  1:6462288/1‑93 (MQ=255)
   atatGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGt                                                                                      <  1:5744090/93‑1 (MQ=255)
   atatGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGt                                                                                      <  2:6263869/93‑1 (MQ=255)
    tatGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTg                                                                                     >  1:3951120/1‑93 (MQ=255)
          aaCAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTa                                                                               >  1:1881413/1‑93 (MQ=255)
           aCAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTAt                                                                              >  2:3569810/1‑93 (MQ=255)
           aCAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTAt                                                                              <  1:4696910/93‑1 (MQ=255)
              aaTCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGc                                                                           >  2:2863331/1‑93 (MQ=255)
              aaTCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGc                                                                           <  1:3728318/93‑1 (MQ=255)
              aaTCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGc                                                                           >  2:6013482/1‑93 (MQ=255)
              aaTCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGc                                                                           >  1:1329796/1‑93 (MQ=255)
               aTCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGcc                                                                          >  2:607713/1‑93 (MQ=255)
                  gCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGt                                                                       <  1:4569717/93‑1 (MQ=255)
                  gCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGt                                                                       <  1:1305841/93‑1 (MQ=255)
                   cAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGtt                                                                      >  2:4796002/1‑93 (MQ=255)
                   cAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGtt                                                                      <  2:6651163/93‑1 (MQ=255)
                    aaaaaTTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTg                                                                     >  2:2403653/1‑93 (MQ=255)
                     aaaaTTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGa                                                                    >  1:2328919/1‑93 (MQ=255)
                       aaTTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGGCCGTGCCATTATCGCCGGTTGAtt                                                                  >  2:5692804/1‑93 (MQ=255)
                       aaTTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGAtt                                                                  >  2:1438182/1‑93 (MQ=255)
                       aaTTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGAtt                                                                  <  1:1631026/93‑1 (MQ=255)
                        aTTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGAttt                                                                 <  2:1500568/93‑1 (MQ=255)
                        aTTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGAttt                                                                 <  2:813225/93‑1 (MQ=255)
                         tttGATTATCAGCATGAAGACTGCACAGGCCCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTa                                                                <  2:1046365/93‑1 (MQ=255)
                         tttGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTa                                                                >  1:6130519/1‑93 (MQ=255)
                         tttGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTa                                                                >  2:3408365/1‑93 (MQ=255)
                         tttGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTa                                                                >  1:6241736/1‑93 (MQ=255)
                           tGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGc                                                              <  2:1735429/93‑1 (MQ=255)
                            gATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCg                                                             >  2:6209618/1‑93 (MQ=255)
                             aTTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGt                                                            <  1:1562942/93‑1 (MQ=255)
                              ttATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTa                                                           <  2:3322281/93‑1 (MQ=255)
                              ttATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTa                                                           >  2:4690693/1‑93 (MQ=255)
                                aTCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATc                                                         >  2:5171907/1‑93 (MQ=255)
                                 tCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCg                                                        <  2:2199950/93‑1 (MQ=255)
                                  cAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCgt                                                       >  1:4593341/1‑92 (MQ=255)
                                   aGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCt                                                      <  1:2204914/93‑1 (MQ=255)
                                   aGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCt                                                      <  1:393410/93‑1 (MQ=255)
                                    gCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCtt                                                     >  1:5735510/1‑93 (MQ=255)
                                      aTGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTa                                                   >  1:6636831/1‑93 (MQ=255)
                                       tGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGGTACTTTCTCGTGGGCCCGTGCCATTATCGCCGGATGATTTAGCGTATCGCTTTAt                                                  >  1:1902302/1‑93 (MQ=255)
                                       tGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTAt                                                  >  2:5533638/1‑93 (MQ=255)
                                       tGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTAt                                                  <  1:6555239/93‑1 (MQ=255)
                                       tGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTAt                                                  >  1:2305187/1‑93 (MQ=255)
                                       tGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTAt                                                  >  1:2679249/1‑93 (MQ=255)
                                       tGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTAt                                                  >  1:472789/1‑93 (MQ=255)
                                        gAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATg                                                 >  1:2376659/1‑93 (MQ=255)
                                         aaGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGc                                                >  1:2680971/1‑93 (MQ=255)
                                         aaGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGc                                                >  2:4608357/1‑93 (MQ=255)
                                           gACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGa                                              >  1:1113861/1‑93 (MQ=255)
                                           gACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGa                                              >  1:335323/1‑93 (MQ=255)
                                             cTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGaaa                                            <  1:5490705/93‑1 (MQ=255)
                                              tGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAg                                           <  1:5856112/93‑1 (MQ=255)
                                              tGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAg                                           <  1:6010589/93‑1 (MQ=255)
                                               gCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAgg                                          >  2:1322363/1‑93 (MQ=255)
                                               gCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAgg                                          >  2:2740415/1‑93 (MQ=255)
                                                cacaGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGCGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGc                                         <  2:4201335/93‑1 (MQ=255)
                                                   aGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGAt                                      <  2:488911/93‑1 (MQ=255)
                                                   aGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGAt                                      >  2:250736/1‑93 (MQ=255)
                                                   aGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGAt                                      <  2:2551657/93‑1 (MQ=255)
                                                    ggACCATCCGCGGATCGATAGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATg                                     >  1:6513579/1‑93 (MQ=255)
                                                    ggACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATg                                     >  1:4314198/1‑93 (MQ=255)
                                                    ggACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATg                                     <  2:6145987/93‑1 (MQ=255)
                                                     gACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGa                                    >  2:5321708/1‑93 (MQ=255)
                                                     gACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGa                                    >  1:3178140/1‑93 (MQ=255)
                                                       ccATCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGaaa                                  <  2:2637382/93‑1 (MQ=255)
                                                          tCAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTg                               <  2:1902775/93‑1 (MQ=255)
                                                           cAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGt                              <  1:1697533/93‑1 (MQ=255)
                                                           cAGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGt                              >  2:2011281/1‑93 (MQ=255)
                                                            aGCGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGtt                             <  1:2410658/93‑1 (MQ=255)
                                                              cGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCg                           >  1:3283758/1‑93 (MQ=255)
                                                              cGGATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCg                           <  2:3695597/93‑1 (MQ=255)
                                                                gATCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCggg                         >  1:2098361/1‑93 (MQ=255)
                                                                 aTCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGa                        <  1:407570/93‑1 (MQ=255)
                                                                  tCGATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGAc                       >  1:4817480/1‑93 (MQ=255)
                                                                    gATGGCGGCTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGc                     <  2:1276017/93‑1 (MQ=255)
                                                                      taatggcTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCAc                   >  1:1381148/5‑93 (MQ=255)
                                                                       ggcggcTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCAct                  >  1:6679313/1‑93 (MQ=255)
                                                                       ggcggcTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCAct                  >  1:3870991/1‑93 (MQ=255)
                                                                         cggcTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCActct                <  2:4178182/93‑1 (MQ=255)
                                                                         cggcTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCActct                <  1:5158400/93‑1 (MQ=255)
                                                                         cggcTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCActct                <  1:3511979/93‑1 (MQ=255)
                                                                          ggcTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCtt               >  2:6487448/1‑93 (MQ=255)
                                                                          ggcTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCtt               >  1:1018483/1‑93 (MQ=255)
                                                                          ggcTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCtt               <  2:3120025/93‑1 (MQ=255)
                                                                          ggcTACTTTCTGGTGGCCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCtt               <  1:4387674/93‑1 (MQ=255)
                                                                           gcTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCCGTCTTc              >  1:2134287/1‑93 (MQ=255)
                                                                           gcTACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTc              >  2:3997116/1‑93 (MQ=255)
                                                                             tACTTTCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCgg            >  1:2995205/1‑93 (MQ=255)
                                                                                tttCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGcc         <  2:4492090/93‑1 (MQ=255)
                                                                                  tCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTg       <  1:4555280/93‑1 (MQ=255)
                                                                                  tCTGGTGGGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTg       >  2:3730252/1‑93 (MQ=255)
                                                                                    tggtggGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCa     <  1:804180/93‑1 (MQ=255)
                                                                                     ggtggGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCAt    <  2:5822231/93‑1 (MQ=255)
                                                                                      gtggGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCAtc   >  2:1008677/1‑93 (MQ=255)
                                                                                       tggGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCAtct  >  2:6698928/1‑93 (MQ=255)
                                                                                       tggGCCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCAtct  <  2:4288384/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTTATATGCTAACAAATCGCAAAAATTTGATTATCAGCATGAAGACTGCACAGGACCATCAGCGGATCGATGGCGGCTACTTTCTGGTGGACCCGTGCCATTATCGCCGGTTGATTTAGCGTATCGCTTTATGCGAAAGGCGATGAAATTGTTCGGGACGCACTCTTCGGGCCTGCATCT  >  NC_000913/1478004‑1478183

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