Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,165,015 T→A 18.7% D320V (GAT→GTT)  parC ← DNA topoisomerase IV, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,165,0150TA18.7% 223.7 / 43.5 98D320V (GAT→GTT) parCDNA topoisomerase IV, subunit A
Reads supporting (aligned to +/- strand):  ref base T (29/50);  new base A (8/10);  total (38/60)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

ACCAGCCATTCGGAGAGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCA  >  NC_000913/3164919‑3165113
                                                                                                |                                                                                                  
acCAGCCATTCGGAGAGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAt                                                                                                    >  2:890097/1‑97 (MQ=255)
acCAGCCATTCGGAGAGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCgg                                                                                                 <  2:491440/100‑1 (MQ=255)
  cAGCCATTCGGAGAGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCggt                                                                                                >  1:2033519/1‑99 (MQ=255)
  cAGCCATTCGGAGAGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAAATCg                                                                                                  >  2:338421/1‑97 (MQ=255)
   agacagtCGGAGAGTATTTCCAGCAGGTTTTGCAGCGCCGGACGACCATCCACACTGATCATATTAAGGTTAATATGATATCTCATTTCCAGATCggt                                                                                                <  2:184151/92‑1 (MQ=255)
   aGCCATTCGGAGAGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCggtgg                                                                                              >  1:869308/1‑100 (MQ=255)
   aGCCATTCGGAGAGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCgg                                                                                                 >  1:1267374/1‑97 (MQ=255)
   aGCCATTCGGAGAGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACggttg                                                                                              <  1:1626865/100‑3 (MQ=255)
    gCCATTCGGACAGGATTTCCAGCAGGTTTTTCACCGCCGGACGCCCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCggtggt                                                                                             <  2:2396561/100‑1 (MQ=255)
     ccATTCGGAGAGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACggtggt                                                                                             <  1:459175/99‑1 (MQ=255)
       aTTCGGAGAGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACGGTGGTAGc                                                                                          <  1:2031143/100‑1 (MQ=255)
        ttCGGAGAGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCg                                                                                         <  2:596510/100‑1 (MQ=255)
           ggAGAGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAag                                                                                      <  1:2108653/100‑1 (MQ=255)
           ggAGAGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACGGGCGTAGCGa                                                                                        >  2:508807/1‑98 (MQ=255)
            gagaGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAaga                                                                                     <  2:213488/100‑1 (MQ=255)
            gagaGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAaga                                                                                     <  2:1269174/100‑1 (MQ=255)
             agagGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAagag                                                                                    <  2:1086487/100‑1 (MQ=255)
             agagGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAagag                                                                                    >  2:11229/1‑100 (MQ=255)
             agagGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAagag                                                                                    <  1:336708/100‑1 (MQ=255)
              gagGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACGGTGGTAGCGAagag                                                                                    >  1:1578104/1‑99 (MQ=255)
                 gATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAggtgg                                                                                <  1:470250/100‑1 (MQ=255)
                  aTTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAggtggt                                                                               <  1:683094/100‑1 (MQ=255)
                  aTTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAggtgg                                                                                <  1:1353931/99‑1 (MQ=255)
                  aTTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACGGTGGTAGCGAAGAggtggt                                                                               <  2:929257/100‑1 (MQ=255)
                   tttCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAggtggt                                                                               >  2:1672242/1‑99 (MQ=255)
                   tttCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGtt                                                                              <  1:548452/100‑1 (MQ=255)
                    ttCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTtc                                                                             <  2:1630773/100‑1 (MQ=255)
                     tCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTtc                                                                             >  2:435157/1‑99 (MQ=255)
                       cagcagGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTtcatn                                                                          <  2:834536/100‑2 (MQ=255)
                       cagcagGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACGGTGGTAGCGAAGAGGTGGTtcatc                                                                          <  1:1942666/100‑1 (MQ=255)
                        agcagGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACGGTGGTAGCGAAGAGGTGGTtcatca                                                                         <  1:1985845/100‑1 (MQ=255)
                          cagGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCAcc                                                                       <  1:1050498/100‑1 (MQ=255)
                          cagGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCAcc                                                                       <  1:215529/100‑1 (MQ=255)
                          cagGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCAcc                                                                       <  1:11230/100‑1 (MQ=255)
                           agGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCt                                                                      <  2:2007560/100‑1 (MQ=255)
                            ggTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTg                                                                     <  1:1299328/100‑1 (MQ=255)
                                tttCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATcc                                                                 <  1:1030648/100‑1 (MQ=255)
                                tttCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATcc                                                                 <  2:129133/100‑1 (MQ=255)
                                   cACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCata                                                              <  2:2033497/100‑1 (MQ=255)
                                   cACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCat                                                               >  1:1296497/1‑99 (MQ=255)
                                   cACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACGGTGGTAGCGAAGAGGAGGTTCATCACCTAATCCata                                                              >  1:1826015/1‑100 (MQ=255)
                                    aCCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCatat                                                             >  1:1488344/1‑100 (MQ=255)
                                    aCCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGACCGGTGGTAGCGAAGAGGTGGTTCATCACCTgc                                                                    >  2:836027/1‑92 (MQ=255)
                                     ccgccgGACTACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATc                                                            <  2:1296483/100‑1 (MQ=255)
                                     ccgccgGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATc                                                            <  1:1774814/100‑1 (MQ=255)
                                     ccgccgGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATc                                                            <  2:1267419/100‑1 (MQ=255)
                                      cgccgGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATcc                                                           <  1:2199408/100‑1 (MQ=255)
                                      cgccgGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATcc                                                           <  2:2096791/100‑1 (MQ=255)
                                       gccgGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCa                                                          >  2:1227594/1‑100 (MQ=255)
                                        ccgGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCa                                                          <  1:15887/99‑1 (MQ=255)
                                           gacgacCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCAcgcg                                                      >  2:1388596/1‑100 (MQ=255)
                                             cgacCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGt                                                    >  1:2168460/1‑100 (MQ=255)
                                              gacCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGtt                                                   <  1:685176/100‑1 (MQ=255)
                                              gacCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGt                                                    <  1:220003/99‑1 (MQ=255)
                                                   tCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAAc                                              >  2:209571/1‑100 (MQ=255)
                                                   tCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAAc                                              <  1:1854133/100‑1 (MQ=255)
                                                   tCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAAc                                              >  2:1057215/1‑100 (MQ=255)
                                                     cAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAAcgc                                            <  2:859789/100‑1 (MQ=255)
                                                     cAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAAcgc                                            >  2:1932750/1‑100 (MQ=255)
                                                       gACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCgg                                          <  2:1007090/100‑1 (MQ=255)
                                                        aCCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGc                                         <  2:2348043/100‑1 (MQ=255)
                                                         ccGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGca                                        <  1:495171/100‑1 (MQ=255)
                                                           gATCATATTATGGTTAATATGATAGCTCTTTTCCAGATCGGTGGTAGCGGAGAGGTGGGTGATCACCTGATCCATATCCACGCGGTTGGAACGCGGc                                         <  1:701527/97‑1 (MQ=255)
                                                           gATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGcaca                                      >  2:1210893/1‑100 (MQ=255)
                                                           gATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGcaca                                      >  1:1864020/1‑100 (MQ=255)
                                                           gATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCAGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATACACGCGGTTGGAACGCGGcaca                                      >  2:212156/1‑100 (MQ=255)
                                                           gATCATATTAAGGTTAATACGATAGCTCTTTTCCAGAACGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGcaca                                      >  1:1294424/1‑100 (MQ=255)
                                                              cATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATc                                   >  1:1693841/1‑100 (MQ=255)
                                                              cATATTAAGGTTAATACGATAACTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAAt                                    >  2:1473978/1‑99 (MQ=255)
                                                               atatTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCa                                  <  1:2337081/100‑1 (MQ=255)
                                                               atatTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCa                                  <  2:661162/100‑1 (MQ=255)
                                                               atatTAAGGTTAATACGATAGCTCTTTTCCAGAACGGTGGTCGCGAAGAGGTGGTTCATCACCTGAACCATATCCACGCGGTTGGAACGCGGCACAATc                                   >  2:1174833/1‑99 (MQ=255)
                                                                tatTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCAc                                 >  1:917993/1‑100 (MQ=255)
                                                                 atTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCAcc                                >  2:703830/1‑100 (MQ=255)
                                                                      ggTTAATACGATAGCTCTTTTCCAGAACGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGcg                           <  1:1025776/100‑1 (MQ=255)
                                                                        ttAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGcggg                         <  2:1204661/100‑1 (MQ=255)
                                                                         tAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGATATGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATAACCAGGcggg                         >  2:2132642/1‑99 (MQ=255)
                                                                           aTACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGcgggtcg                      >  1:1954840/1‑100 (MQ=255)
                                                                            tACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGcgggtcgg                     >  1:1806353/1‑100 (MQ=255)
                                                                            tACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAAGCACCAGGcgggtcgg                     >  1:2370991/1‑100 (MQ=255)
                                                                              cGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGcgggtcgggt                   >  1:1723267/1‑100 (MQ=255)
                                                                              cGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGcgggtcgggt                   >  2:1098934/1‑100 (MQ=255)
                                                                               gATAGCTCTTTTCCAAATCGGTGGTAGCGAATAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCATATTCACCAGAcgggtcgggt                   <  2:318346/99‑1 (MQ=255)
                                                                                  aGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTtctc               <  2:1336732/100‑1 (MQ=255)
                                                                                  aGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTtctc               <  2:32349/100‑1 (MQ=255)
                                                                                  aGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTtct                >  1:1012502/1‑99 (MQ=255)
                                                                                   gCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACGTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCg              <  2:1682891/100‑1 (MQ=255)
                                                                                    ctctTTTCCATATCGGTGGTAGCGAAAAGGTAGTTCATCACCTGATCCATATTCACGCGGTTGGAACGCGGCACAATCACAAGGCGGGACGGGATCTCGt             >  2:905313/1‑100 (MQ=255)
                                                                                     tctTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTg            <  2:2202504/100‑1 (MQ=255)
                                                                                      ctTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTgg           >  1:826435/1‑100 (MQ=255)
                                                                                       ttttCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAAGGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGt          <  1:1691688/100‑1 (MQ=255)
                                                                                       ttttCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGt          <  1:34995/100‑1 (MQ=255)
                                                                                          tCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAg       <  1:2083510/100‑1 (MQ=255)
                                                                                          tCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAg       <  1:1847447/100‑1 (MQ=255)
                                                                                          tCCAGAACGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAg       <  1:767725/100‑1 (MQ=255)
                                                                                            cAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGAt     <  2:1648022/100‑1 (MQ=255)
                                                                                            cAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGAt     <  1:1458250/100‑1 (MQ=255)
                                                                                             aGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATGTCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGAt     <  2:1728195/99‑1 (MQ=255)
                                                                                               aTCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCa  >  2:1052961/1‑100 (MQ=255)
                                                                                                |                                                                                                  
ACCAGCCATTCGGAGAGGATTTCCAGCAGGTTTTTCACCGCCGGACGACCATCCAGACCGATCATATTAAGGTTAATACGATAGCTCTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCA  >  NC_000913/3164919‑3165113

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