Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 993,254 A→G 11.7% intergenic (+20/+23) pepN → / ← ssuB aminopeptidase N/aliphatic sulfonate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,2540AG11.7% 314.3 / 28.0 123intergenic (+20/+23)pepN/ssuBaminopeptidase N/aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (51/56);  new base G (6/8);  total (58/65)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-01

Only 100 of 126 total aligned reads displayed.

ACCGCTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCAA  >  NC_000913/993105‑993402
                                                                                                                                                     |                                                                                                                                                    
aCCGCTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCa                                                                                                                                                      >  1:364295/1‑150 (MQ=255)
 ccGCTGATTCGCCTGAAACGTTACGATGCCATACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGATAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCaa                                                                                                                                                     >  1:831761/1‑150 (MQ=255)
 ccGCTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCaa                                                                                                                                                     >  2:269239/1‑150 (MQ=255)
  cGCTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGGACGAGAAGATAACTAAAGCACTTGCTTGATAAATAACAGAATGGCGGCAAt                                                                                                                                                    >  2:284099/1‑150 (MQ=255)
   gCTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCCAATGGCGGCAATa                                                                                                                                                   >  1:257387/1‑150 (MQ=255)
    cTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGATAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAg                                                                                                                                                  >  1:396489/1‑150 (MQ=255)
     tGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAgc                                                                                                                                                 <  2:895646/150‑1 (MQ=255)
      gATTCGCCTGAAACTTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCAGGAACAGTTGAAAGGGCTGGAAAATCTCTATGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAgcg                                                                                                                                                >  2:488001/1‑150 (MQ=255)
              tGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCAtt                                                                                                                                        <  2:1225475/150‑1 (MQ=255)
                 aaCGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCgg                                                                                                                                     >  2:1112657/1‑150 (MQ=255)
                  aCGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTCggt                                                                                                                                    <  1:62074/150‑2 (MQ=255)
                  aCGTTACGAGGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCggg                                                                                                                                    <  1:1452829/150‑1 (MQ=255)
                   cGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCgggg                                                                                                                                   >  2:813043/1‑150 (MQ=255)
                   cGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGGTGAAAGTGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCgggg                                                                                                                                   <  1:1368714/150‑1 (MQ=255)
                    gTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTCggtta                                                                                                                                  <  1:1249022/150‑4 (MQ=255)
                      tACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAAt                                                                                                                                <  1:539855/150‑1 (MQ=255)
                       aCGATGCCAAACGGCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAAtt                                                                                                                               <  1:313976/150‑1 (MQ=255)
                           tGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAcc                                                                                                                           <  1:544083/150‑1 (MQ=255)
                           tGCCAAACATCAGGAGAAAATGCGCGCGGCGCTGGAACTGTTGAAAGTGCGGCAAAATCTCTCTGGCGATCTGTAGGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAcc                                                                                                                           <  2:1350413/150‑1 (MQ=255)
                                aaCGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGt                                                                                                                      >  2:1466791/1‑150 (MQ=255)
                                  cGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCATTAGCGCCGCTATTCGGGGAATTTACCCCGGttt                                                                                                                    >  1:843821/1‑150 (MQ=255)
                                  cGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGttt                                                                                                                    >  1:1289712/1‑150 (MQ=255)
                                    tctggAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTTAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTtc                                                                                                                  <  2:758034/147‑1 (MQ=255)
                                      aGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTtctc                                                                                                                <  2:686071/150‑1 (MQ=255)
                                       ggAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCa                                                                                                               <  2:968571/150‑1 (MQ=255)
                                        gagaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAg                                                                                                              <  1:215748/150‑1 (MQ=255)
                                         agaAAATGCGCTCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAgg                                                                                                             <  2:1718429/150‑1 (MQ=255)
                                          gaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGc                                                                                                            <  1:1366130/150‑1 (MQ=255)
                                           aaaaTGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCg                                                                                                           <  1:431894/150‑1 (MQ=255)
                                             aaTGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAg                                                                                                         >  1:688013/1‑150 (MQ=255)
                                              aTGCGCGCGGCGCTGGAACAGTTGATAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCAGTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGt                                                                                                        >  1:255927/1‑150 (MQ=255)
                                              aTGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGt                                                                                                        <  2:1477639/150‑1 (MQ=255)
                                              aTGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGt                                                                                                        <  1:632072/150‑1 (MQ=255)
                                              aTGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGt                                                                                                        <  1:1403926/150‑1 (MQ=255)
                                                   cgcggcgcggGAACAGTTGAAAGGGCTGGAAACTCTCTCTGGCGATCTGTACGAGAAGAAAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAg                                                                                                   <  2:272613/150‑1 (MQ=255)
                                                   cgcgGCGCTGGAACAGTTTAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAg                                                                                                   <  2:611567/150‑1 (MQ=255)
                                                   cgcgGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAg                                                                                                   <  2:1698278/150‑1 (MQ=255)
                                                    gcggcgCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTTTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATACCGCCGCCATTCGGGGAATTAACCCCTGATTTCTCAGGCGAGTGTCAGa                                                                                                  >  1:144098/1‑150 (MQ=255)
                                                      ggcgCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGAtt                                                                                                >  1:383399/1‑150 (MQ=255)
                                                        cgcTGGAACTGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTca                                                                                              <  2:477121/150‑1 (MQ=255)
                                                            ggAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcaccac                                                                                          <  1:1431715/150‑1 (MQ=255)
                                                             gAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACg                                                                                         >  1:44551/1‑150 (MQ=255)
                                                              aaCAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGt                                                                                        >  2:1516166/1‑150 (MQ=255)
                                                                  gTTGACAGGGCTGGAAATAGTCTCTTGCGTTCTGTACGAGAATATAACTAAAGCACTGTCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCa                                                                                    <  2:819551/150‑1 (MQ=255)
                                                                  gTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCa                                                                                    <  2:364295/150‑1 (MQ=255)
                                                                     gAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCa                                                                                 >  2:476544/1‑150 (MQ=255)
                                                                         ggtcTGGAAAATCTCTCAGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAAAAGAGCCGACATTCTGGGAATATACACCTGTTTTCTCAGGCGAGTTTCACATACAACCCGTTGCATCACCCg                                                                             >  2:697171/4‑150 (MQ=255)
                                                                              ggAAAATCTCTCTGGCGATCTGTACGTGAAGATAACTAAAGCCCTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCa                                                                        <  1:925779/150‑1 (MQ=255)
                                                                                  aaTCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAacac                                                                    >  1:716240/1‑150 (MQ=255)
                                                                                    tctctcTGGCGTTCTGTACGAGAATATAACTAAAGCACTTGCTTGATAAATAACCGAATGGCGGGAATAGCGCAGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCt                                                                  <  2:275060/150‑1 (MQ=255)
                                                                                      tctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTcc                                                                <  2:352462/150‑1 (MQ=255)
                                                                                       ctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCg                                                               >  1:1427878/1‑150 (MQ=255)
                                                                                          ttgCGATCTGTACGAGAAGATAACTAAAGCTATGGCTTGATAAATAACTGAATGGCGGCAATAGCGCAGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCtt                                                            <  1:284099/148‑1 (MQ=255)
                                                                                            gCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTcc                                                          <  2:568463/150‑1 (MQ=255)
                                                                                            gCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACAACGTTGCATCACCCGTTGCAACACCTCCGCTTcc                                                          <  2:911732/150‑1 (MQ=255)
                                                                                                tCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCTGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCt                                                      >  2:1219070/1‑150 (MQ=255)
                                                                                                 cTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTc                                                     <  2:831761/150‑1 (MQ=255)
                                                                                                  tGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGATGCAACACATCCGCTTCCAGCTcc                                                    >  2:767486/1‑150 (MQ=255)
                                                                                                   gTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCg                                                   >  2:30024/1‑150 (MQ=255)
                                                                                                   gTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCg                                                   >  1:1575915/1‑150 (MQ=255)
                                                                                                      cGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTa                                                <  2:990069/150‑1 (MQ=255)
                                                                                                      cGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCCAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTa                                                <  2:257387/150‑1 (MQ=255)
                                                                                                         gaagaTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATc                                             <  2:520406/150‑1 (MQ=255)
                                                                                                         gaagaTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGAGCCGACATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATc                                             <  2:656561/150‑1 (MQ=255)
                                                                                                           agaTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCCCCCGTTGCAACCCCTCCGCTTCCAGCTCCGCTAATCGc                                           >  2:955738/1‑150 (MQ=255)
                                                                                                            gaTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGcc                                          >  2:1715538/1‑149 (MQ=255)
                                                                                                            gaTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCAGCCATTCGGGGAATTAAAACATGTTTTCTCAGGCTAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCa                                          >  2:131926/1‑150 (MQ=255)
                                                                                                              tAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCAcc                                        <  2:1021659/150‑1 (MQ=255)
                                                                                                              aaaCTAACGCACTGGCTTGATAAATACCCGACTGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGTATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCAcc                                        <  1:276779/149‑1 (MQ=255)
                                                                                                                 cTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTTATCGCACCGAc                                     <  2:396489/150‑1 (MQ=255)
                                                                                                                 cTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAc                                     >  2:1703238/1‑150 (MQ=255)
                                                                                                                 cTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAc                                     <  1:645695/150‑1 (MQ=255)
                                                                                                                   aaaGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTAGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAccc                                   <  2:1288804/150‑1 (MQ=255)
                                                                                                                    aaGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCt                                  >  1:1027772/1‑150 (MQ=255)
                                                                                                                      gCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTaa                                >  1:659917/1‑150 (MQ=255)
                                                                                                                            gCTTGATAAATAACGGAATGGCGGCAATAGAGACGGCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAcg                          <  2:858749/150‑1 (MQ=255)
                                                                                                                                gATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGa                      >  1:1117104/1‑150 (MQ=255)
                                                                                                                                gATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGa                      >  2:1149200/1‑150 (MQ=255)
                                                                                                                                gATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATAGCACCGCCCCTAATCGACGCGGa                      >  2:509198/1‑150 (MQ=255)
                                                                                                                                gATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATGTACCCCTGGTTTCTCAGGCGAGTTTCAGAGTCACCACGTTGGATCACCCGTTGCAACACCTCCGCTGCCAGCTCCGCTAATCGCTCCGCCCCTGATCGACGCggc                      >  1:116991/1‑149 (MQ=255)
                                                                                                                                    aataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACggg                  <  1:1227977/150‑1 (MQ=255)
                                                                                                                                      taaCCGAATGGCGGCTATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGCCCCTAATCGACGCGGACggggg                >  1:1289414/1‑150 (MQ=255)
                                                                                                                                      taaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACggggg                <  2:1669622/150‑1 (MQ=255)
                                                                                                                                       ccccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACGCCTGATTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGCACGGGTGa               >  1:1553238/3‑150 (MQ=255)
                                                                                                                                        cccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCGAATCGGCGCGGACGGGGGat              >  1:1630228/2‑150 (MQ=255)
                                                                                                                                         ccGAATGGCGGCGCTATTGCCGCCATTCGGGGAATTTAACCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCTTCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGata             >  1:1141315/1‑150 (MQ=255)
                                                                                                                                         ccGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTTCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGata             <  1:400073/150‑1 (MQ=255)
                                                                                                                                            aaTGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATcc          >  2:456793/1‑150 (MQ=255)
                                                                                                                                             aTGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAc                           <  1:1008713/132‑1 (MQ=255)
                                                                                                                                             aTGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAc                           <  1:313740/132‑1 (MQ=255)
                                                                                                                                             aTGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAc                           >  2:313740/1‑132 (MQ=255)
                                                                                                                                             aTGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAc                           >  2:1008713/1‑132 (MQ=255)
                                                                                                                                              tGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCAc        >  1:1408325/1‑150 (MQ=255)
                                                                                                                                               ggcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACt       >  2:363889/1‑150 (MQ=255)
                                                                                                                                                gcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTg      <  1:1112657/150‑1 (MQ=255)
                                                                                                                                                 cggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGCTTCACATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGt     <  1:1078661/150‑1 (MQ=255)
                                                                                                                                                  ggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTc    >  2:1406514/1‑150 (MQ=255)
                                                                                                                                                   gcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCa   <  1:1516166/150‑1 (MQ=255)
                                                                                                                                                   gcAAAAGCGCCGCAATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCACATTCACCCCGTTGCATCACCCGATGCAACACCTCCGCGGCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCa   >  1:1053036/1‑150 (MQ=255)
                                                                                                                                                    cAATAGAGCCGCCATTAGGGGAATTTACACCTGTTTTCTCAGGCGAGTTTAATATTCACCACGTTGCATAACACGTTGCAACACCTCCGCTTCCAGCTCCGCTAAACGCACCGAAACTAAACGACGCGGACGGCGGGTATCCACtctcta  >  2:708660/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACCGCTGATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCAA  >  NC_000913/993105‑993402

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