Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 324,876 C→A 6.6% R163S (AGG→AGT ykgH ← uncharacterized protein YkgH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913324,8760CA 6.6% 304.2 / 8.4 91R163S (AGG→AGTykgHuncharacterized protein YkgH
Reads supporting (aligned to +/- strand):  ref base C (39/46);  new base A (3/3);  total (42/49)
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.62e-01

ACCATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTCAGCAGT  >  NC_000913/324727‑325025
                                                                                                                                                     |                                                                                                                                                     
aCCATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTc                                                                                                                                                       >  2:716209/1‑150 (MQ=255)
aCCATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACAAGTTc                                                                                                                                                       >  1:1086527/1‑150 (MQ=255)
 ccATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTcc                                                                                                                                                      <  1:716209/150‑1 (MQ=255)
 ccATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTcc                                                                                                                                                      <  2:254047/150‑1 (MQ=255)
 ccATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTcc                                                                                                                                                      <  1:654918/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                                                                                                                                             >  1:390415/1‑150 (MQ=255)
          tGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGt                                                                                                                                             >  2:978311/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:373880/1‑150 (MQ=255)
                  cAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCAtt                                                                                                                                     >  1:67081/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:1049200/150‑1 (MQ=255)
                                                        cAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTg                                                                                               <  1:1103791/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                        >  2:724159/1‑150 (MQ=255)
                                                                                                                               cGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCAcgc                        >  1:133356/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         >  2:701796/1‑150 (MQ=255)
                                                                                                                                              aCCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCtt         >  1:1120275/1‑150 (MQ=255)
                                                                                                                                               ccAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTc        <  2:135054/150‑1 (MQ=255)
                                                                                                                                               ccAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTc        >  2:718498/1‑150 (MQ=255)
                                                                                                                                                 aGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTcag      >  2:875958/1‑150 (MQ=255)
                                                                                                                                                  gTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTcagc     <  1:488553/150‑1 (MQ=255)
                                                                                                                                                  gTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTACCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTcagc     <  2:639730/150‑1 (MQ=255)
                                                                                                                                                     ccTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAAATGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATAACATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTCAGCAGt  >  2:1249180/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCATTTTTTTGACAAGGCAAAACATTACAGAAATTACAATGCTTAGAAAAAATCCCAGAGCCAGGCAAATGGCAGCAGGGTATTTATTTTTTATTTCTGGTGAGCTTATAGTTACACAGCCATCACCGTTCGGATACAATGACCAGTTCCTGGTCAGGTATATCATTTTTGTATAAAGCTCACTGTTAATTGCTGAAGGCGTTTGTATCTCATCGAGAACTTGCCTGCCATTATGATATTCTGACGCATACCATTTGATAATAGAGGCCTGCACGCGAGATATTTCCTCTTCAGCAGT  >  NC_000913/324727‑325025

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: