Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 291 A→T 5.0% intergenic (+36/‑46) thrL → / → thrA thr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132910AT 5.0% 372.4 / 6.9 130intergenic (+36/‑46)thrL/thrAthr operon leader peptide/Bifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (63/58);  new base T (2/5);  total (67/63)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.99e-01

Only 100 of 132 total aligned reads displayed.

GGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAC  >  NC_000913/145‑440
                                                                                                                                                  |                                                                                                                                                      
    tAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCg                                                                                                                                                >  1:1443797/1‑150 (MQ=255)
     aGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCgg                                                                                                                                               <  2:419974/150‑1 (MQ=255)
      gcgcTCAGACAGATAAAAATTAAAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTAACA‑GTGCggg                                                                                                                                              <  1:1793090/150‑1 (MQ=255)
       cgcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGc                                                                                                                                             <  2:1219737/150‑1 (MQ=255)
        gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCt                                                                                                                                            >  2:2010942/1‑150 (MQ=255)
        gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCt                                                                                                                                            >  1:704359/1‑150 (MQ=255)
        gcACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCt                                                                                                                                            >  2:1112913/1‑150 (MQ=255)
          acagacagATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCttt                                                                                                                                          <  2:633070/150‑1 (MQ=255)
          acagacagATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGCttt                                                                                                                                          <  1:1677215/150‑1 (MQ=255)
          acagacagATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGGGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCttt                                                                                                                                          <  1:156458/150‑1 (MQ=255)
           cagacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGCCGCGTCCAGGAAACACAGAAAAAAGCCCGCACCTGACC‑GTGCGGGCtttt                                                                                                                                         >  2:1004121/1‑150 (MQ=255)
            agacagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCCCAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGCttttt                                                                                                                                        <  1:1548930/150‑1 (MQ=255)
              acagaTAAAAATTACACAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCttttttt                                                                                                                                      <  2:1243508/150‑1 (MQ=255)
               cagaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCtttttttt                                                                                                                                     <  2:1615167/150‑1 (MQ=255)
                 taTAAAAATTACAGAGAACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTc                                                                                                                                   <  1:521657/149‑1 (MQ=255)
                 gaTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTc                                                                                                                                   >  1:1409295/1‑150 (MQ=255)
                  aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCCCCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAACCACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCg                                                                                                                                  <  1:384387/150‑1 (MQ=255)
                  aTAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACTGGTTACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCTCCAGACA‑GTGCGGGCtttttttttag                                                                                                                                  >  2:182607/1‑148 (MQ=255)
                    aaaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCA‑CTGTCAGGTGCGGGCttttttctgtgtt                                                                                                                                <  1:1663288/150‑8 (MQ=38)
                     aaaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGAcc                                                                                                                               <  1:1281506/150‑1 (MQ=255)
                      aaaTTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGACTACAGGAAACACAGAAAAAAGCCCGCACGAGCCA‑GTGCGGGCGTTTTTTTTCGAcct                                                                                                                              >  1:346480/1‑149 (MQ=255)
                         ttACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACATAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAg                                                                                                                           <  1:591399/150‑1 (MQ=255)
                         ttACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGAAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAg                                                                                                                           <  1:393461/150‑1 (MQ=255)
                            cAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACATGTAACGGTGCGGGCTGCCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGTCCCAAGGTa                                                                                                                        <  1:1176116/150‑1 (MQ=255)
                               agTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACg                                                                                                                     <  2:473257/150‑1 (MQ=255)
                                     caacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGCCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaa                                                                                                               >  2:121554/1‑150 (MQ=255)
                                     caacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTTCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaa                                                                                                               >  2:1718975/1‑150 (MQ=255)
                                     caacaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaa                                                                                                               >  1:478478/1‑150 (MQ=255)
                                       acaTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaaca                                                                                                             >  2:602880/1‑150 (MQ=255)
                                        caTCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTaacaa                                                                                                            >  1:396385/1‑150 (MQ=255)
                                             tttaaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACACAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                       <  1:693134/147‑1 (MQ=255)
                                             aTGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATg                                                                                                       >  1:1909037/1‑150 (MQ=255)
                                                aaaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGa                                                                                                    >  2:1399930/1‑150 (MQ=255)
                                                 aaCGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAg                                                                                                   >  2:669956/1‑150 (MQ=255)
                                                   cGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACCACCATGCGAgtg                                                                                                 >  1:1946571/1‑150 (MQ=255)
                                                      aTTAGCACCACCATTACCACCACCATCACCGTTACCTGGGGTAATGGTGCGGGCTGACGCGTACAGTAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTg                                                                                              <  2:1429639/150‑1 (MQ=255)
                                                      aTTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGACGTAACAACCATGCGAGTGTTg                                                                                              >  1:381169/1‑150 (MQ=255)
                                                        tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                            <  1:453319/150‑1 (MQ=255)
                                                        tAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGaa                                                                                            >  1:1936404/1‑150 (MQ=255)
                                                         aGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTACGACCAAAGGTAACGAGGTAACAACCATCCCAGTGTTGaaa                                                                                           >  1:1867240/1‑149 (MQ=255)
                                                          gCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGt                                                                                          >  1:700487/1‑150 (MQ=255)
                                                             ccaccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTGTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcg                                                                                       >  1:530428/1‑150 (MQ=255)
                                                             ccaccaTTACCACCACCATCACCATTACCACAGGTAAAGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCAACGGACA‑GTGCGGGCTTTTTTTTTCGGCCAAAGGGAACGAGGTAACGACCAAGGGAGTGTTGAAGTTcg                                                                                       >  2:1608458/1‑150 (MQ=255)
                                                              caccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCAACTTACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGCGGTGACAACCATGCGAGTGTTGAAGTTcgg                                                                                      >  2:110401/1‑150 (MQ=255)
                                                               accaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggc                                                                                     <  2:34485/150‑1 (MQ=255)
                                                                ccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTcggcg                                                                                    <  2:1798848/150‑1 (MQ=255)
                                                                ccaTTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAAGCGAGTGTTGAAGTTcggcg                                                                                    >  2:851768/1‑150 (MQ=255)
                                                                     accaccaccATCACCATTACCACAGGTAACCGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTCACA‑GTGCGGGCATTTTTTTTCGACCAAAGGTAAAGACGTAACAACCAAGCGAGTGTTGAAGTTCGGCGGTACa                                                                               >  2:1737317/1‑150 (MQ=255)
                                                                        accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACCACCATGCGAGTGTTGAAGTTCGGCGGTACCTCa                                                                            >  2:568595/1‑150 (MQ=255)
                                                                        accaccaTCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCa                                                                            >  2:1922899/1‑150 (MQ=255)
                                                                            ccatcaccatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTgg                                                                        >  2:825577/1‑150 (MQ=255)
                                                                             catcaccatTACCCCAGGTAACGGTGCGGGCTGCCGCGTACAGGAAACACAGAAAAAAGCCCGCCCCAGACA‑CTGCGGGCATTTTTTTTCGACCAAAAGGAACGAGGTAAACACCATGCGAGTGTTGAAGTTCGGCGGTTCCTCAGTGGc                                                                       >  1:1301271/1‑150 (MQ=255)
                                                                                 accatTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCAAGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCaaaa                                                                   >  2:1092699/1‑149 (MQ=255)
                                                                                   catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                 >  2:653278/1‑150 (MQ=255)
                                                                                   catTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGc                                                                 <  1:1751305/150‑1 (MQ=255)
                                                                                      tACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGa                                                              >  2:554914/1‑150 (MQ=255)
                                                                                       aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGaa                                                             <  1:657518/150‑1 (MQ=255)
                                                                                       aCCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGAGGGCTTTTTTTATCTACCAAACGTAACGAGGTAACAACCATGCGAGTGATGAAGTTCGGCGGTACATCAGTGTCAAATGCAGaa                                                             >  1:322798/1‑150 (MQ=255)
                                                                                        ccACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAc                                                            <  1:838910/150‑1 (MQ=255)
                                                                                            aGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGttt                                                        >  2:761872/1‑150 (MQ=255)
                                                                                             ggTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGgttt                                                       >  2:1472199/1‑150 (MQ=255)
                                                                                              gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                      <  1:1370065/150‑1 (MQ=255)
                                                                                              gTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTc                                                      <  1:1366488/150‑1 (MQ=255)
                                                                                               tAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCt                                                     <  1:806905/150‑1 (MQ=255)
                                                                                                aaCGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTg                                                    <  2:1372547/150‑1 (MQ=255)
                                                                                                  cGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCg                                                  <  2:1443796/150‑1 (MQ=255)
                                                                                                  cGGTGCGGGCGGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCg                                                  <  1:1360706/150‑1 (MQ=255)
                                                                                                   ggTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCgt                                                 >  1:1730658/1‑150 (MQ=255)
                                                                                                      gCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGtt                                              >  2:555577/1‑150 (MQ=255)
                                                                                                       cGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTg                                             >  2:1830935/1‑150 (MQ=255)
                                                                                                         ggCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGcc                                           <  1:977986/150‑1 (MQ=255)
                                                                                                          gCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGCCA‑GTGAGGGCTTTTTTTTTAAACAAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCTGCGGTACATCAGTGGCAAATGAAGAACGTATTCTGCGTTTTGACg                                          >  1:486249/1‑150 (MQ=255)
                                                                                                              aCGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGatat                                      <  1:2002271/150‑1 (MQ=255)
                                                                                                               cgcgTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGTTAtt                                     >  2:450988/1‑150 (MQ=255)
                                                                                                                   tACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCTAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTgg                                 <  2:1266419/150‑1 (MQ=255)
                                                                                                                    aCAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGa                                >  2:165873/1‑150 (MQ=255)
                                                                                                                     cAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGCGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGAAATTCTGGaa                               >  1:1404750/1‑150 (MQ=255)
                                                                                                                       ggAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAg                             >  2:71702/1‑150 (MQ=255)
                                                                                                                       ggAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAAGTTTTCTGCGTGTTGCCGATATTCTGGAAAg                             <  1:1084969/150‑1 (MQ=255)
                                                                                                                        aaaaCACACAAAAAAGCCCGCACCCGACA‑GTGCGGGCTTTTTTTTTCGATCAAAGGTAACGAGGTAACAACCATGCGAGTGTGGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGc                            <  2:1036984/149‑1 (MQ=255)
                                                                                                                            cacaGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAAAGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATg                        >  2:283846/1‑150 (MQ=255)
                                                                                                                              caGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTAGAACAAAGGTAACGACTGAAAAACCAAGCGAGTGTTGAAGTTCTGCGGGACATAAGTGGCAAATGCAGAACGTTTTCACCGTGTTGCCGAGATTCTGGAAAGCAATGcc                      >  2:1905666/1‑150 (MQ=255)
                                                                                                                               aGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTGGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCca                     >  2:137579/1‑150 (MQ=255)
                                                                                                                                 aaaaaaGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCATGTTGCCGATATTCTGGAAAGCAATGCcagg                   >  2:1384425/1‑150 (MQ=255)
                                                                                                                                  aaaaaGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTGGCCGATATTCTGGAAAGCAATGCccggc                  >  2:286391/1‑150 (MQ=255)
                                                                                                                                  aaaaaGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGAGTTGAAGTTCGGCGGGACATCAGTGGCAAATGCAGAACGTTTTCGGCGGGTTGCCGCTATTCTGAAAACCAATGCcagac                  >  1:417419/1‑148 (MQ=255)
                                                                                                                                   aaaaGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggca                 >  1:1141556/1‑150 (MQ=255)
                                                                                                                                    aaaGCCCGCACCTGACA‑GTGCGGGC‑TTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg               >  1:1925185/1‑150 (MQ=255)
                                                                                                                                     aaGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg               >  1:1949193/1‑150 (MQ=255)
                                                                                                                                     aaGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCcaggcagg               <  1:182607/150‑1 (MQ=255)
                                                                                                                                       gCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg             <  1:121554/150‑1 (MQ=255)
                                                                                                                                       gCCCGCACCGGACA‑GTGCGGGCTTTTTTTTTATACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGAAGAACTTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcagggg             <  2:1818186/150‑1 (MQ=255)
                                                                                                                                         ccGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAggcaggggca           >  1:920365/1‑150 (MQ=255)
                                                                                                                                             aCCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACACCCATGCGAGTGTTGAAGTGCGGCTGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAACGCAATGCCAGGCAGGGGCAGGTg       <  2:1238396/150‑1 (MQ=255)
                                                                                                                                             aCCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg       <  2:1231863/150‑1 (MQ=255)
                                                                                                                                             aCCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTg       <  2:1038950/150‑1 (MQ=255)
                                                                                                                                                tGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGcc    >  1:1979041/1‑150 (MQ=255)
                                                                                                                                                 gACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGccc   >  1:1611523/1‑149 (MQ=255)
                                                                                                                                                  aCA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAc  <  1:1145600/150‑1 (MQ=255)
                                                                                                                                                  aCA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAc  <  2:970616/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GGCATAGCGCACAGACAGATAAAAATTACAGAGTACACAACATCCATGAAACGCATTAGCACCACCATTACCACCACCATCACCATTACCACAGGTAACGGTGCGGGCTGACGCGTACAGGAAACACAGAAAAAAGCCCGCACCTGACA‑GTGCGGGCTTTTTTTTTCGACCAAAGGTAACGAGGTAACAACCATGCGAGTGTTGAAGTTCGGCGGTACATCAGTGGCAAATGCAGAACGTTTTCTGCGTGTTGCCGATATTCTGGAAAGCAATGCCAGGCAGGGGCAGGTGGCCAC  >  NC_000913/145‑440

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: ‑