Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,392,498 C→G 5.9% G182R (GGC→CGC)  nuoM ← NADH:quinone oxidoreductase subunit M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,392,4980CG 5.9% 283.6 / 6.6 85G182R (GGC→CGC) nuoMNADH:quinone oxidoreductase subunit M
Reads supporting (aligned to +/- strand):  ref base C (38/41);  new base G (3/2);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

AAGGCGATGAAGAAGCCCAGCATTAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAA  >  NC_000913/2392350‑2392647
                                                                                                                                                    |                                                                                                                                                     
aacgCGATGAAGAAGCCCAGCATTAACAGGTATTCCACACCAGTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCa                                                                                                                                                      <  2:1526578/147‑1 (MQ=255)
      aTGAAGAAGCCCAGCATTAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGcc                                                                                                                                                <  2:879503/150‑1 (MQ=255)
       tGAAGAAGCCCAGCATTAACAGGTATTCCACACCACTGGACATTGGCGTAATAAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCt                                                                                                                                               >  2:213375/1‑150 (MQ=255)
          agaagCCCAGCATTAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTggg                                                                                                                                            <  2:832622/150‑1 (MQ=255)
          agaagCCCAGCATTAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTAGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTggg                                                                                                                                            >  1:1191176/1‑150 (MQ=255)
            aagCCCAGCATTAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGgtg                                                                                                                                          <  2:1514388/150‑1 (MQ=255)
             agCCCAGCATTAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGgtgt                                                                                                                                         >  2:1480053/1‑150 (MQ=255)
              gCCCAGCATTAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTa                                                                                                                                        >  1:1431033/1‑150 (MQ=255)
               cccAGCATTAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCGAGTGGCCTggtgatgt                                                                                                                                       <  2:1517261/150‑7 (MQ=255)
                ccAGCATTAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTaaa                                                                                                                                      >  1:843513/1‑150 (MQ=255)
                 cAGCATTAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAAt                                                                                                                                     <  1:62656/150‑1 (MQ=255)
                  aGCATTAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCATCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAATACCAGCGCCAGGATGTCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATg                                                                                                                                    <  1:1473161/150‑1 (MQ=255)
                  aGCATTAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATg                                                                                                                                    <  1:831807/150‑1 (MQ=255)
                      ttAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATgaaga                                                                                                                                <  2:1140835/150‑1 (MQ=255)
                      ttAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATgaaga                                                                                                                                <  2:728765/150‑1 (MQ=255)
                          cAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACtt                                                                                                                            >  2:1109097/1‑150 (MQ=255)
                              tATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGtt                                                                                                                        <  2:741553/150‑1 (MQ=255)
                                      accacTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGt                                                                                                                <  2:670415/150‑1 (MQ=255)
                                        cacTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGACGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGa                                                                                                              >  1:528320/1‑150 (MQ=255)
                                          cTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATa                                                                                                            <  1:547770/150‑1 (MQ=255)
                                            ggACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATAcg                                                                                                          >  2:275391/1‑150 (MQ=255)
                                               cATTGGCGTATTCAGCAGCGCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACGAGGCCACTCGCCTGTGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGt                                                                                                       <  2:1212236/150‑1 (MQ=255)
                                                       tATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGt                                                                                               <  1:568697/150‑1 (MQ=255)
                                                        aTTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTc                                                                                              <  1:952223/150‑1 (MQ=255)
                                                         ttCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCa                                                                                             >  1:860103/1‑150 (MQ=255)
                                                            agcagcTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCagag                                                                                          >  2:1477623/1‑150 (MQ=255)
                                                                gcTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCtt                                                                                      <  1:440490/150‑1 (MQ=255)
                                                                  tctTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTa                                                                                    <  1:775165/150‑1 (MQ=255)
                                                                   ctTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTAt                                                                                   >  2:229064/1‑150 (MQ=255)
                                                                   ctTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTAt                                                                                   <  1:474152/150‑1 (MQ=255)
                                                                     tCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGc                                                                                 >  2:950821/1‑150 (MQ=255)
                                                                       aTAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGccc                                                                               <  1:146339/150‑1 (MQ=255)
                                                                          gTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCca                                                                            >  1:1099964/1‑150 (MQ=255)
                                                                             gAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAg                                                                         >  2:513285/1‑150 (MQ=255)
                                                                             gAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAg                                                                         >  2:720634/1‑150 (MQ=255)
                                                                                   ccAGACGCCGGTCGCGTTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGAt                                                                   <  2:1358346/150‑1 (MQ=255)
                                                                                    cAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATc                                                                  <  1:275391/150‑1 (MQ=255)
                                                                                    cAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATc                                                                  <  2:1191176/150‑1 (MQ=255)
                                                                                     aGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCa                                                                 <  2:860103/150‑1 (MQ=255)
                                                                                       aCGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGCTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAgg                                                               <  1:563401/150‑1 (MQ=255)
                                                                                        cGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGa                                                              <  1:190819/150‑1 (MQ=255)
                                                                                        cGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGa                                                              >  1:1476255/1‑150 (MQ=255)
                                                                                        cGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGa                                                              >  1:1069311/1‑150 (MQ=255)
                                                                                         gCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGa                                                              <  1:417022/149‑1 (MQ=255)
                                                                                         gCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGa                                                              >  2:417022/1‑149 (MQ=255)
                                                                                          ccGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAg                                                            <  2:99618/150‑1 (MQ=255)
                                                                                          ccGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAg                                                            >  1:513522/1‑150 (MQ=255)
                                                                                            ggTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTa                                                          <  1:1044839/150‑1 (MQ=255)
                                                                                                gCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATc                                                      <  2:762787/150‑1 (MQ=255)
                                                                                                 cATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCg                                                     >  1:415863/1‑150 (MQ=255)
                                                                                                    tGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCa                                                  <  1:85410/150‑1 (MQ=255)
                                                                                                      tAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCAcc                                                >  2:41232/1‑150 (MQ=255)
                                                                                                      tAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCAcc                                                <  1:213375/150‑1 (MQ=255)
                                                                                                       aGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCa                                               >  2:70218/1‑150 (MQ=255)
                                                                                                       aGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCa                                               <  2:843513/150‑1 (MQ=255)
                                                                                                            aCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGcatc                                          >  1:1340034/1‑150 (MQ=255)
                                                                                                            aCAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGcatc                                          >  1:1244544/1‑150 (MQ=255)
                                                                                                             cAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGcatca                                         >  1:1255175/1‑150 (MQ=255)
                                                                                                             cAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGcatca                                         >  2:1340396/1‑150 (MQ=255)
                                                                                                              aaaaaCCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGcatcat                                        >  1:1190821/1‑150 (MQ=255)
                                                                                                                 aaCCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTc                                     >  1:1064376/1‑150 (MQ=255)
                                                                                                                  aCCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTcc                                    >  1:199952/1‑150 (MQ=255)
                                                                                                                     aGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCag                                 <  2:1190821/150‑1 (MQ=255)
                                                                                                                     aGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCag                                 <  1:225616/150‑1 (MQ=255)
                                                                                                                       cgcCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCagaa                               <  1:1109097/150‑1 (MQ=255)
                                                                                                                          cAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTcc                                    >  1:689993/1‑142 (MQ=255)
                                                                                                                          cAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTcc                                    <  2:689993/142‑1 (MQ=255)
                                                                                                                              aTGGCAATCAACATCACCAGGCCACTCGCCTGTGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGTGTTTTTATGCCGCCTGAGTGCTATCAGGTAGTACATCGTGACCAGTATCATTTGCCagaagtagatg                        <  2:1530702/150‑3 (MQ=255)
                                                                                                                               tcgCAATCAACATCACAAGGCAACTCCCCTGGGTATAAATGAATAACTTGCTTGCCGACGTAATACGAGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCagaagaagaaga                       >  2:1346481/3‑150 (MQ=255)
                                                                                                                               tGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCagaagaagaaga                       >  1:782504/1‑150 (MQ=255)
                                                                                                                                 gCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAAc                     <  1:456171/150‑1 (MQ=255)
                                                                                                                                 gCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCAGAAGAAGAAGAACa                    >  1:902994/1‑150 (MQ=255)
                                                                                                                                  cAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACa                    <  2:513522/150‑1 (MQ=255)
                                                                                                                                    aTCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAgg                  >  2:1397905/1‑150 (MQ=255)
                                                                                                                                       aaCATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAAc               >  2:422860/1‑150 (MQ=255)
                                                                                                                                           cacccAGGCGAGTGGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGt           >  1:282800/4‑150 (MQ=255)
                                                                                                                                            cACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCa                                               >  1:1229929/1‑113 (MQ=255)
                                                                                                                                            cACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCa                                               <  2:1229929/113‑1 (MQ=255)
                                                                                                                                            acccAGGCGAGTGGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCagaagaag                           >  1:173769/3‑133 (MQ=255)
                                                                                                                                            acccAGGCGAGTGGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCagaagaag                           <  2:173769/131‑1 (MQ=255)
                                                                                                                                             cccAGGCGAGTGGCCTGGGTGTAAATGAAGCACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGCGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCg         >  2:392650/2‑150 (MQ=255)
                                                                                                                                              ccAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGa        >  1:1151821/1‑150 (MQ=255)
                                                                                                                                              ccAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGa        >  1:778509/1‑150 (MQ=255)
                                                                                                                                               cAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGAt       >  1:151034/1‑150 (MQ=255)
                                                                                                                                                aGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATg      <  1:1480053/150‑1 (MQ=255)
                                                                                                                                                   ccACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCa   >  2:92518/1‑150 (MQ=255)
                                                                                                                                                   ccACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCa   >  1:597551/1‑150 (MQ=255)
                                                                                                                                                    cACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCaa  <  2:1255175/150‑1 (MQ=255)
                                                                                                                                                    cACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGCCACCAGCAACATTTCCCAGAAGAAGAAGAACAGGACCATGTAGATGGCaa  >  2:729760/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAGGCGATGAAGAAGCCCAGCATTAACAGGTATTCCACACCACTGGACATTGGCGTATTCAGCAGCTCTTCATAGTTGAAGGTCCAGACGCCGGTCGCATTGTAGTGAACAAAAACCAGCGCCAGGATGGCAATCAACATCACCAGGCCACTCGCCTGGGTGTAAATGAAGAACTTGGTTGCCGCCGTGATACGCGTTTTACCGTCAGAGGCTTTATGCCCCCACAGTGCGATCAGGAAGTACATCGGCACCAGCATCATTTCCCAGAAGAAGAAGAACAGGAACATGTCGATGGCAA  >  NC_000913/2392350‑2392647

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: