Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 324,875 T→A 6.2% N164Y (AAC→TAC) ‡ ykgH ← uncharacterized protein YkgH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913324,8750TA 6.2% 305.1 / 7.6 96N164Y (AAC→TAC) ‡ykgHuncharacterized protein YkgH
Reads supporting (aligned to +/- strand):  ref base T (41/49);  new base A (3/3);  total (44/52)
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.66e-01

TACCATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTCAGC  >  NC_000913/324726‑325022
                                                                                                                                                     |                                                                                                                                                   
tACCATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGtt                                                                                                                                                     >  1:1245548/1‑150 (MQ=255)
 aCCATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTc                                                                                                                                                    >  2:716209/1‑150 (MQ=255)
 aCCATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACAAGTTc                                                                                                                                                    >  1:1086527/1‑150 (MQ=255)
  ccATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTcc                                                                                                                                                   <  1:654918/150‑1 (MQ=255)
  ccATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTcc                                                                                                                                                   <  1:716209/150‑1 (MQ=255)
  ccATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTcc                                                                                                                                                   <  2:254047/150‑1 (MQ=255)
   cATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCt                                                                                                                                                  <  1:23280/150‑1 (MQ=255)
   cATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCt                                                                                                                                                  >  1:1249266/1‑150 (MQ=255)
    aTTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTg                                                                                                                                                 >  2:756342/1‑150 (MQ=255)
       tttttGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGTAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTc                                                                                                                                              >  2:231211/1‑150 (MQ=255)
       tttttGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTc                                                                                                                                              <  2:1330218/150‑1 (MQ=255)
         tttGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAg                                                                                                                                            <  2:1002386/150‑1 (MQ=255)
         tttGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGGAACTGGTCat                                                                                                                                            <  2:43451/150‑2 (MQ=255)
          ttGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAgg                                                                                                                                           <  1:422844/150‑1 (MQ=255)
           tGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGt                                                                                                                                          >  2:978311/1‑150 (MQ=255)
           tGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGt                                                                                                                                          >  1:390415/1‑150 (MQ=255)
            gACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGta                                                                                                                                         >  2:235840/1‑150 (MQ=255)
              cAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGtata                                                                                                                                       <  2:366418/150‑1 (MQ=255)
              cAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGGAACTGGTCattgtat                                                                                                                                       <  2:1229427/150‑7 (MQ=255)
               aaGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGGAACTGGTCattgtatc                                                                                                                                      <  1:125762/150‑8 (MQ=255)
                   cAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCAtt                                                                                                                                  >  1:67081/1‑150 (MQ=255)
                   cAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCAtt                                                                                                                                  >  1:373880/1‑150 (MQ=255)
                       aCATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTg                                                                                                                              <  2:1139365/150‑1 (MQ=255)
                        cATTAAAGAAATTACAATGCTTAGAAAAAATCCGAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGATCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGt                                                                                                                             <  1:496953/150‑1 (MQ=255)
                         aTTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGta                                                                                                                            <  1:337929/150‑1 (MQ=255)
                          ttACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGtat                                                                                                                           <  2:506326/150‑1 (MQ=255)
                          ttACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGtat                                                                                                                           <  2:381481/150‑1 (MQ=255)
                          ttACAGAAATTACAATGCTTAGAAAAAATCCCAGAGACAGGCAAATGGAAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTAGTCAGGTATATCATTTTTGtat                                                                                                                           >  2:1158509/1‑150 (MQ=255)
                           tACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGtata                                                                                                                          >  1:849178/1‑150 (MQ=255)
                            aCAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGCCCAGTTCCTGGTCAGGTATATCATTTTTGTATaa                                                                                                                         >  2:1486883/1‑150 (MQ=255)
                            aCAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATaa                                                                                                                         >  2:1271111/1‑150 (MQ=255)
                              aGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATCTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAg                                                                                                                       <  1:1069963/150‑1 (MQ=255)
                                aaaTTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCt                                                                                                                     >  1:59766/1‑150 (MQ=255)
                                  aTTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCa                                                                                                                   <  2:1479270/150‑1 (MQ=255)
                                    tACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACt                                                                                                                 <  1:93533/150‑1 (MQ=255)
                                    tACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACt                                                                                                                 >  2:236952/1‑150 (MQ=255)
                                           cTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAAtt                                                                                                          >  1:455195/1‑150 (MQ=255)
                                             tAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGc                                                                                                        <  1:946581/150‑1 (MQ=255)
                                              aGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCt                                                                                                       >  2:1269253/1‑150 (MQ=255)
                                              aGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCt                                                                                                       <  2:794891/150‑1 (MQ=255)
                                               gAAAAAATCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGa                                                                                                     <  1:336973/150‑1 (MQ=255)
                                                  aaaaTCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACAAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAg                                                                                                   >  2:1227954/1‑150 (MQ=255)
                                                     aTCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCg                                                                                                <  1:235840/150‑1 (MQ=255)
                                                       cccAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGtt                                                                                              <  2:1249266/150‑1 (MQ=255)
                                                       cccAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGtt                                                                                              >  1:364597/1‑150 (MQ=255)
                                                       cccAGAGCCAGGCAAATGGCAGAAGGGTATTTATTTTTTATTTCTCGTGAGCTTATAGTTACACACCAATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATTAAGCTCACTGTTAATTGCTGAAGGCGtt                                                                                              >  2:611279/1‑150 (MQ=255)
                                                        ccAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGttt                                                                                             >  2:153230/1‑150 (MQ=255)
                                                         cAGAGCCATGCAAATGGCAGCATGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTg                                                                                            <  1:1429939/150‑1 (MQ=255)
                                                         cAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTg                                                                                            <  1:1103791/150‑1 (MQ=255)
                                                         cAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTg                                                                                            <  1:1049200/150‑1 (MQ=255)
                                                           gagCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGGCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTa                                                                                          <  1:1319634/150‑1 (MQ=255)
                                                             gCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTAtc                                                                                        <  1:629200/150‑1 (MQ=255)
                                                              ccAGGCAAATGGCAGAAGGTTATTTATTTTTTATTTCTGGTGTTCTTATAGTTACTCAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTAtct                                                                                       <  2:806543/150‑1 (MQ=255)
                                                                aGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCa                                                                                     <  2:373880/150‑1 (MQ=255)
                                                                   cAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTAACTGTTAATTGCTGAAGGCGTTTGTATCTCATCg                                                                                  >  1:775764/1‑150 (MQ=255)
                                                                      aTGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCgaga                                                                               >  2:1396993/1‑150 (MQ=255)
                                                                        ggCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAAc                                                                             >  1:1512741/1‑150 (MQ=255)
                                                                           agcagGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTtg                                                                          <  1:1491887/150‑1 (MQ=255)
                                                                               gggTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTtgcctg                                                                      <  2:59766/150‑1 (MQ=255)
                                                                               gggTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTtgcctg                                                                      <  1:514852/150‑1 (MQ=255)
                                                                                  tatttattTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCa                                                                   <  2:59014/150‑1 (MQ=255)
                                                                                    tttattttttatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCAtt                                                                 <  1:1271111/150‑1 (MQ=255)
                                                                                     ttattttttatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTa                                                                <  2:748110/150‑1 (MQ=255)
                                                                                             tatttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTc                                                        <  2:83835/150‑1 (MQ=255)
                                                                                               tttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTg                                                      >  1:70661/1‑150 (MQ=255)
                                                                                                ttCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGa                                                     <  1:1486496/150‑1 (MQ=255)
                                                                                                       gAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTAGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATAc                                              <  1:1158509/150‑1 (MQ=255)
                                                                                                         gCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCa                                            >  2:279202/1‑150 (MQ=255)
                                                                                                          cTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCAt                                           <  1:978311/150‑1 (MQ=255)
                                                                                                           ttATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCAtt                                          <  2:455195/150‑1 (MQ=255)
                                                                                                            tataGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCAttt                                         >  1:1229774/1‑150 (MQ=255)
                                                                                                                 ttACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGatag                                    >  1:3134/1‑149 (MQ=255)
                                                                                                                   acacAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGataata                                  >  2:1292609/1‑150 (MQ=255)
                                                                                                                   acacAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGataata                                  <  2:67081/150‑1 (MQ=255)
                                                                                                                       aGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAgg                              <  2:1245548/150‑1 (MQ=255)
                                                                                                                        gCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCg                                                                                  <  2:99962/97‑1 (MQ=255)
                                                                                                                        gCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCg                                                                                  >  1:99962/1‑97 (MQ=255)
                                                                                                                         ccATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGcc                            <  1:564916/150‑1 (MQ=255)
                                                                                                                          cATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCt                           <  1:351197/150‑1 (MQ=255)
                                                                                                                           aTCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTg                          >  1:373322/1‑150 (MQ=255)
                                                                                                                           aTCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTAGCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTg                          >  1:639730/1‑150 (MQ=255)
                                                                                                                                cGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCAcgc                     >  1:133356/1‑150 (MQ=255)
                                                                                                                                cGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCAcgc                     >  2:724159/1‑150 (MQ=255)
                                                                                                                                 gTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCAcgcg                    <  1:756342/150‑1 (MQ=255)
                                                                                                                                   tCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCgag                  <  1:1396993/150‑1 (MQ=255)
                                                                                                                                     ggATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGat                >  1:135054/1‑150 (MQ=255)
                                                                                                                                      gATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGata               >  2:759872/1‑150 (MQ=255)
                                                                                                                                        atacctGACCAGGAACTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATAtt             >  1:478921/6‑150 (MQ=255)
                                                                                                                                           cctGACCAGGAACTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTcc          >  1:98807/3‑150 (MQ=255)
                                                                                                                                           cctGACCAGGAACTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTcc          >  1:505424/3‑150 (MQ=255)
                                                                                                                                               aCCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCtt      >  1:1120275/1‑150 (MQ=255)
                                                                                                                                               aCCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCtt      >  2:701796/1‑150 (MQ=255)
                                                                                                                                                ccAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTc     >  2:718498/1‑150 (MQ=255)
                                                                                                                                                ccAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTc     <  2:135054/150‑1 (MQ=255)
                                                                                                                                                  aGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTcag   >  2:875958/1‑150 (MQ=255)
                                                                                                                                                   gTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTcagc  <  1:488553/150‑1 (MQ=255)
                                                                                                                                                   gTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTACCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTcagc  <  2:639730/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TACCATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTCAGC  >  NC_000913/324726‑325022

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