Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,178,668 A→G intergenic (+165/‑155) rplL → / → rpoB 50S ribosomal subunit protein L7/L12/RNA polymerase beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,178,6680AG99.6% 819.6 / ‑4.6 260intergenic (+165/‑155)rplL/rpoB50S ribosomal subunit protein L7/L12/RNA polymerase beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (137/122);  minor base C (1/0);  total (138/122)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 263 total aligned reads displayed.

TGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACAAAA  >  NC_000913/4178578‑4178756
                                                                                          |                                                                                        
tGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCg                                                                                        >  1:1122991/1‑93 (MQ=255)
tGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCg                                                                                        <  1:1347917/93‑1 (MQ=255)
 gTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGt                                                                                       <  1:2620343/93‑1 (MQ=255)
 gTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGt                                                                                       >  1:1539684/1‑93 (MQ=255)
  tAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTc                                                                                      <  1:2973817/93‑1 (MQ=255)
      gcgcCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGtt                                                                                  >  2:3937253/1‑93 (MQ=255)
       cgcCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTc                                                                                 <  1:948356/93‑1 (MQ=255)
       cgcCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTc                                                                                 <  1:3874703/93‑1 (MQ=255)
         ccAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCtg                                                                               >  2:2542611/1‑93 (MQ=255)
          cAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCtgt                                                                              <  2:2871524/93‑1 (MQ=255)
          cAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCtgt                                                                              >  2:643191/1‑93 (MQ=255)
           aGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCtgtg                                                                             <  1:4278313/93‑1 (MQ=255)
            gTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTgtgt                                                                            >  1:1928437/1‑93 (MQ=255)
            gTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTgtgt                                                                            >  1:1270402/1‑93 (MQ=255)
             tAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTa                                                                           >  2:2063524/1‑93 (MQ=255)
             tAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTa                                                                           >  1:3076421/1‑93 (MQ=255)
              aGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTaa                                                                          <  1:1145027/93‑1 (MQ=255)
               gCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTaaa                                                                         >  1:4550647/1‑93 (MQ=255)
                cGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAAt                                                                        <  1:1584906/93‑1 (MQ=255)
                cGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAAt                                                                        <  2:2242480/93‑1 (MQ=255)
                 gTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATc                                                                       >  1:200197/1‑93 (MQ=255)
                 gTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATc                                                                       >  2:4279143/1‑93 (MQ=255)
                 gTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATc                                                                       <  1:3748007/93‑1 (MQ=255)
                  tttCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCg                                                                      >  1:4429947/1‑93 (MQ=255)
                    tCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCa                                                                    <  1:1716272/93‑1 (MQ=255)
                     cacaCTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCaa                                                                   >  2:583041/1‑93 (MQ=255)
                        acTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGa                                                                <  2:111533/93‑1 (MQ=255)
                           gTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAAt                                                             >  2:1089904/1‑93 (MQ=255)
                           gTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAAt                                                             >  1:1440833/1‑93 (MQ=255)
                           gTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAAt                                                             >  2:1308553/1‑93 (MQ=255)
                            tttGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATg                                                            >  1:3272684/1‑93 (MQ=255)
                            tttGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATg                                                            >  1:3349235/1‑93 (MQ=255)
                            tttGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATg                                                            >  1:352076/1‑93 (MQ=255)
                              tGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGt                                                          <  1:4075153/93‑1 (MQ=255)
                                actactGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGttt                                                        <  2:3715600/93‑1 (MQ=255)
                                 ctactGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCTGTTCTGTGTAAATCGCAATGAAATGGTTTa                                                       >  2:4459562/1‑93 (MQ=255)
                                 ctactGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTa                                                       >  1:2027347/1‑93 (MQ=255)
                                  tactGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTaa                                                      <  2:4072504/93‑1 (MQ=255)
                                   actGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAg                                                     >  2:3645303/1‑93 (MQ=255)
                                    ctgctgTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGc                                                    <  2:3597267/93‑1 (MQ=255)
                                    ctgctgTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGc                                                    >  2:4025516/1‑93 (MQ=255)
                                      gctgTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGt                                                  <  2:1613014/93‑1 (MQ=255)
                                        tgtgCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGa                                                >  1:3370861/1‑93 (MQ=255)
                                        tgtgCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGa                                                >  2:2002558/1‑93 (MQ=255)
                                        tgtgCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGa                                                >  2:1204762/1‑93 (MQ=255)
                                        tgtgCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGa                                                >  1:3557220/1‑93 (MQ=255)
                                         gtgCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGAt                                               <  2:3451702/93‑1 (MQ=255)
                                          tgCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATa                                              <  2:4365711/93‑1 (MQ=255)
                                          tgCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATa                                              <  1:2566364/93‑1 (MQ=255)
                                            ccTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGc                                            >  2:4167846/1‑93 (MQ=255)
                                            ccTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGc                                            <  1:684364/93‑1 (MQ=255)
                                            ccTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGc                                            <  2:1783522/93‑1 (MQ=255)
                                             cTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCa                                           <  1:207681/93‑1 (MQ=255)
                                              tttCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCaa                                          >  1:4474409/1‑93 (MQ=255)
                                               ttCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAAc                                         >  2:790540/1‑93 (MQ=255)
                                               ttCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAAc                                         <  2:341053/93‑1 (MQ=255)
                                               ttCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAAc                                         >  2:4556293/1‑93 (MQ=255)
                                                tCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACa                                        <  1:378147/93‑1 (MQ=255)
                                                tCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACa                                        <  2:3760521/93‑1 (MQ=255)
                                                tCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACa                                        <  2:2628155/93‑1 (MQ=255)
                                                tCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACa                                        >  2:3018141/1‑93 (MQ=255)
                                                 cAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAg                                       <  2:4524155/93‑1 (MQ=255)
                                                   aTGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGc                                     <  1:4466968/93‑1 (MQ=255)
                                                    tGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCa                                    >  1:1476614/1‑93 (MQ=255)
                                                     gCTTGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCAt                                   <  1:559008/93‑1 (MQ=255)
                                                       ttGTTTCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTg                                 <  2:3129028/93‑1 (MQ=255)
                                                           ttCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGa                             >  1:1037053/1‑93 (MQ=255)
                                                           ttCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGa                             >  1:1963666/1‑93 (MQ=255)
                                                            tCTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGaa                            <  2:3780104/93‑1 (MQ=255)
                                                             cTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGaaa                           <  1:1203457/93‑1 (MQ=255)
                                                             cTATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGaaa                           >  1:2667075/1‑93 (MQ=255)
                                                              tATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAg                          <  2:273371/93‑1 (MQ=255)
                                                              tATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAg                          >  1:1808781/1‑93 (MQ=255)
                                                              tATCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAg                          >  1:1172898/1‑93 (MQ=255)
                                                               aTCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAgt                         <  2:1700953/93‑1 (MQ=255)
                                                                tCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAgtg                        >  1:4365733/1‑93 (MQ=255)
                                                                tCGACGACTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAgtg                        <  2:3488820/93‑1 (MQ=255)
                                                                  gacgacTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATGGCGGAAAGTGtt                      <  2:2914114/93‑1 (MQ=255)
                                                                   acgacTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTc                     >  1:855800/1‑93 (MQ=255)
                                                                   acgacTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTc                     <  2:3476842/93‑1 (MQ=255)
                                                                   acgacTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTc                     >  2:4242939/1‑93 (MQ=255)
                                                                   acgacTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTc                     >  1:2843578/1‑93 (MQ=255)
                                                                    cgacTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTcc                    <  2:2560837/93‑1 (MQ=255)
                                                                      acTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAt                  <  1:616365/93‑1 (MQ=255)
                                                                      acTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAt                  <  2:1458623/93‑1 (MQ=255)
                                                                      acTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAt                  >  1:3970727/1‑93 (MQ=255)
                                                                      acTTAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAt                  <  1:2478153/93‑1 (MQ=255)
                                                                        ttAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAttt                >  1:3619347/1‑93 (MQ=255)
                                                                         tAATATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCAtttt               >  1:1319444/1‑93 (MQ=255)
                                                                          aaTATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTc              >  2:138859/1‑93 (MQ=255)
                                                                          aaTATACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTc              >  1:2354779/1‑93 (MQ=255)
                                                                           atatACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTcc             <  1:4644694/93‑1 (MQ=255)
                                                                           atatACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTcc             >  2:95953/1‑93 (MQ=255)
                                                                           atatACTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGGTCCATTTTcc             >  2:1478438/1‑93 (MQ=255)
                                                                            tataCTGCGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCg            >  1:1700424/1‑93 (MQ=255)
                                                                              taCTGCGACAGGCCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGt          >  1:2416764/1‑93 (MQ=255)
                                                                                   cGACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTcaaca     <  2:2783031/93‑1 (MQ=255)
                                                                                    gACAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTcaacaa    >  1:2179203/1‑93 (MQ=255)
                                                                                     aCAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACaaa   >  1:303394/1‑93 (MQ=255)
                                                                                      cAGGGCGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACaaaa  >  2:3279937/1‑93 (MQ=255)
                                                                                          |                                                                                        
TGTAAGGCGCCAGTAGCGTTTCACACTGTTTGACTACTGCTGTGCCTTTCAATGCTTGTTTCTATCGACGACTTAATATACTGCGACAGGACGTCCGTTCTGTGTAAATCGCAATGAAATGGTTTAAGCGTGATAGCAACAGGCATTGCGGAAAGTGTTCCATTTTCCGGTCAACAAAA  >  NC_000913/4178578‑4178756

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