Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,001,421 2 bp→GC coding (381‑382/903 nt) rarD ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,001,4210CG100.0% 748.4 / NA 207G128R (GGA→CGA) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (102/105);  total (102/105)
*NC_0009134,001,4220GC99.5% 732.1 / NA 207D127E (GAC→GAGrarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (101/104);  minor base T (0/1);  total (102/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

TAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTT  >  NC_000913/4001334‑4001508
                                                                                       ||                                                                                      
 aataGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGcc                                                                                   <  2:4919690/93‑1 (MQ=255)
  ataGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCg                                                                                  <  2:3213908/93‑1 (MQ=255)
  ataGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCg                                                                                  <  2:1959899/93‑1 (MQ=255)
  ataGGTAGCGAACCAAAAGTCCACAGCTGGAATAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCg                                                                                  >  2:5884364/1‑93 (MQ=255)
   taGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGa                                                                                 <  1:4598849/93‑1 (MQ=255)
     ggTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAgg                                                                               >  1:11580/1‑93 (MQ=255)
     ggTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAgg                                                                               <  1:62683/93‑1 (MQ=255)
     ggTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAgg                                                                               <  1:3578871/93‑1 (MQ=255)
       tAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaa                                                                             <  2:5802821/93‑1 (MQ=255)
         gCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaaaa                                                                           <  2:3756626/93‑1 (MQ=255)
          cGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGCAAat                                                                          >  1:2717425/1‑93 (MQ=255)
          cGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAat                                                                          <  2:5097259/93‑1 (MQ=255)
           gacccAAAAGTCCACAGCTGGACTAACACGCCACATATCTCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatc                                                                         <  1:188335/90‑1 (MQ=255)
           gAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatc                                                                         <  2:3832326/93‑1 (MQ=255)
              ccAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatcatc                                                                      >  2:850101/1‑93 (MQ=255)
              ccAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatcatc                                                                      >  1:2875781/1‑93 (MQ=255)
              acAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatcatc                                                                      <  1:4977155/92‑1 (MQ=255)
               cAAAAGTCCACAGCTGTACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcc                                                                     <  1:698502/93‑1 (MQ=255)
               cAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatcatac                                                                     >  1:320132/1‑91 (MQ=255)
               cAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcc                                                                     >  1:4753205/1‑93 (MQ=255)
                aaaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATccc                                                                    >  1:1956998/1‑93 (MQ=255)
                aaaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATccc                                                                    >  2:4401759/1‑93 (MQ=255)
                aaaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATAGCTCGCCGAGGAAAATCATccc                                                                    >  1:5174534/1‑93 (MQ=255)
                  aaGTCCCCAGCTGGACTAACACGCCACATATCGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCa                                                                  <  1:347895/93‑1 (MQ=255)
                        ccaGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACaa                                                            <  2:1440839/92‑1 (MQ=255)
                        acaGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACaa                                                            >  2:360685/1‑93 (MQ=255)
                        acaGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACaa                                                            >  1:4909774/1‑93 (MQ=255)
                         caGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAAt                                                           >  1:4368788/1‑93 (MQ=255)
                          aGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATg                                                          >  2:5402296/1‑93 (MQ=255)
                           gCTGGACTAACACGCCACATATCGCCAGAATCAACGCCAGCCATTGCATCCGGCGGAATCGATCGCCGAGGAAAATCATCCCCAGCACAATGt                                                         >  2:323203/1‑93 (MQ=255)
                           acTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGt                                                         <  1:6123509/92‑1 (MQ=255)
                            cTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGtt                                                        <  1:346242/93‑1 (MQ=255)
                            cTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGtt                                                        <  2:6003558/93‑1 (MQ=255)
                            cTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGtt                                                        <  1:2922318/93‑1 (MQ=255)
                              ggACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCa                                                      >  2:457272/1‑93 (MQ=255)
                              ggACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCa                                                      >  2:4966885/1‑93 (MQ=255)
                              ggACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCa                                                      >  1:3657336/1‑93 (MQ=255)
                               gACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAc                                                     >  2:2015610/1‑93 (MQ=255)
                                aCTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAcc                                                    <  1:4142407/93‑1 (MQ=255)
                                  tAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAg                                                  <  2:5934542/93‑1 (MQ=255)
                                  tAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAg                                                  >  2:727250/1‑93 (MQ=255)
                                  tAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAg                                                  >  1:774843/1‑93 (MQ=255)
                                  tAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCACCACAATGTTCACCAg                                                  <  2:5834679/93‑1 (MQ=255)
                                   aaCACGCCTCATATCGCTAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGc                                                 <  2:238062/93‑1 (MQ=255)
                                   aaCACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCTCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGc                                                 <  2:4688028/93‑1 (MQ=255)
                                     cacGCCAGATATCGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCgg                                               <  1:963603/93‑1 (MQ=255)
                                     cacGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCgg                                               >  1:149547/1‑93 (MQ=255)
                                        gCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGtt                                            <  2:5641768/93‑1 (MQ=255)
                                        cccACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGtt                                            <  2:2172110/92‑1 (MQ=255)
                                         ccACATATCGCGAGAATCGCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCACCCCCAGCACAATGTTCACCAGCGGGTta                                           <  2:3568118/93‑1 (MQ=255)
                                         ccACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGATta                                           <  1:6021191/93‑1 (MQ=255)
                                           acaTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaat                                         >  1:815517/1‑93 (MQ=255)
                                           acaTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaat                                         >  1:883898/1‑93 (MQ=255)
                                           acaTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaat                                         >  1:2303260/1‑93 (MQ=255)
                                            caTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaata                                        <  2:112794/93‑1 (MQ=255)
                                             atatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaataa                                       <  2:1093837/93‑1 (MQ=255)
                                             atatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaataa                                       <  1:1891929/93‑1 (MQ=255)
                                             atacCGCAAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaataa                                       <  1:4419068/93‑1 (MQ=255)
                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCGTCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                                      <  2:1661176/93‑1 (MQ=255)
                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                                      >  2:1281148/1‑93 (MQ=255)
                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                                      <  2:2834462/93‑1 (MQ=255)
                                                tCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGt                                    <  1:5943444/93‑1 (MQ=255)
                                                tCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATAATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGt                                    <  2:3018721/93‑1 (MQ=255)
                                                 cGCCAGGATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTa                                   <  2:2034206/93‑1 (MQ=255)
                                                  gCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTaa                                  >  1:1910587/1‑93 (MQ=255)
                                                  cccAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTaa                                  <  1:1857229/92‑1 (MQ=255)
                                                   ccAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAc                                 <  2:3977596/93‑1 (MQ=255)
                                                    cAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAcc                                <  2:4625536/93‑1 (MQ=255)
                                                    cAGAATCACCGCCAGCCATTGCATCCGGCTGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAcc                                >  1:3886360/1‑93 (MQ=255)
                                                    cAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAcc                                <  2:315160/93‑1 (MQ=255)
                                                      gAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCaa                              >  1:4300379/1‑93 (MQ=255)
                                                      gAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCaa                              <  1:5277161/93‑1 (MQ=255)
                                                      gAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCaa                              <  1:56293/93‑1 (MQ=255)
                                                      gAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCaa                              >  1:2707085/1‑93 (MQ=255)
                                                       aaTCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAg                             >  1:1032102/1‑93 (MQ=255)
                                                        aTCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAgg                            >  2:2295916/1‑93 (MQ=255)
                                                         tCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGc                           >  1:167990/1‑93 (MQ=255)
                                                         tCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGc                           >  1:557464/1‑93 (MQ=255)
                                                          ccccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCt                          <  1:3838055/91‑1 (MQ=255)
                                                          cACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCt                          <  2:5059843/93‑1 (MQ=255)
                                                          cACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCt                          <  2:487945/93‑1 (MQ=255)
                                                            ccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcg                        <  1:5760719/93‑2 (MQ=255)
                                                             cGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcgc                       >  2:6059368/1‑91 (MQ=255)
                                                                 agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                   >  2:5216623/1‑93 (MQ=255)
                                                                  gccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCa                  >  1:266380/1‑93 (MQ=255)
                                                                  gccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCa                  >  1:847098/1‑93 (MQ=255)
                                                                   ccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAg                 >  1:2355647/1‑93 (MQ=255)
                                                                   ccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAg                 >  1:2883739/1‑93 (MQ=255)
                                                                     aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa               >  1:5273974/1‑93 (MQ=255)
                                                                      ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat              >  2:1359258/1‑93 (MQ=255)
                                                                      ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat              <  2:1562238/93‑1 (MQ=255)
                                                                      ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat              <  1:1240198/93‑1 (MQ=255)
                                                                      ttGCATCCGGAGGAATCGCTCGCCGAGTAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat              <  1:1483689/93‑1 (MQ=255)
                                                                       tGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCata             <  1:129761/93‑1 (MQ=255)
                                                                        gCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCatat            >  1:4080229/1‑93 (MQ=255)
                                                                          aTCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATACAGGCACCAAGGCTCGCTTCCAGCATATgt          >  2:3296401/1‑92 (MQ=255)
                                                                          aTCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg          <  1:3342344/93‑1 (MQ=255)
                                                                           tCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGt         <  2:4271262/93‑1 (MQ=255)
                                                                                 ggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGt   <  1:72612/93‑1 (MQ=255)
                                                                                 ggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGt   >  2:274199/1‑93 (MQ=255)
                                                                                       ||                                                                                      
TAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTT  >  NC_000913/4001334‑4001508

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