Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,214,511 G→C 8.6% P322R (CCT→CGT)  bluF ← blue light‑responsive regulator of BluR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,214,5110GC 8.6% 226.5 / 10.6 70P322R (CCT→CGT) bluFblue light‑responsive regulator of BluR
Reads supporting (aligned to +/- strand):  ref base G (27/37);  new base C (2/4);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

ACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATG  >  NC_000913/1214362‑1214660
                                                                                                                                                     |                                                                                                                                                     
acaTCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAg                                                                                                                                                       >  1:173866/1‑150 (MQ=255)
   tCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCACGCt                                                                                                                                                    >  2:787930/1‑150 (MQ=255)
        tcttctGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAg                                                                                                                                               <  1:1085349/150‑1 (MQ=255)
            ctGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTg                                                                                                                                           >  1:568449/1‑150 (MQ=255)
            ctGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTg                                                                                                                                           <  2:141007/150‑1 (MQ=255)
               gtgtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACa                                                                                                                                        <  1:1308539/150‑1 (MQ=255)
               gtgtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCACGCTTCCAGCCTGACa                                                                                                                                        <  2:50832/150‑1 (MQ=255)
               gtgtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCACGCTTCCAGCCTGACa                                                                                                                                        <  1:945424/150‑1 (MQ=255)
                 gtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCACGCTTCCAGCCTGACAga                                                                                                                                      <  1:601407/150‑2 (MQ=255)
                   tGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAg                                                                                                                                    <  2:156523/150‑1 (MQ=255)
                     ccACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTa                                                                                                                                  <  2:835276/150‑1 (MQ=255)
                       aCGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAAc                                                                                                                                <  2:43052/150‑1 (MQ=255)
                       aCGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAAc                                                                                                                                <  1:938047/150‑1 (MQ=255)
                        cGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTGACa                                                                                                                               <  2:1111744/150‑1 (MQ=255)
                        cGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACa                                                                                                                               <  2:124253/150‑1 (MQ=255)
                         gCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCACGCTTCCAGCCTGACAGAATTaaaat                                                                                                                              <  1:787930/150‑4 (MQ=255)
                           ccATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAAc                                                                                                                            <  1:1252177/150‑1 (MQ=255)
                                     gACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGt                                                                                                                                   <  2:1029850/133‑1 (MQ=255)
                                     gACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGt                                                                                                                                   >  1:1029850/1‑133 (MQ=255)
                                       cTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTg                                                                                                                >  2:655802/1‑150 (MQ=255)
                                          gAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTgaga                                                                                                             >  2:1278862/1‑147 (MQ=255)
                                            aTTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTgg                                                                                                                                                  <  2:33057/111‑1 (MQ=255)
                                            aTTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTgg                                                                                                                                                  >  1:33057/1‑111 (MQ=255)
                                                 aaGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATg                                                                                                      >  2:1245263/1‑150 (MQ=255)
                                                       tGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAAt                                                                                                >  1:1350522/1‑150 (MQ=255)
                                                       tGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAAt                                                                                                >  1:1062787/1‑150 (MQ=255)
                                                          aCAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGc                                                                                             <  1:576701/150‑1 (MQ=255)
                                                          aCAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGc                                                                                             <  2:1349337/150‑1 (MQ=255)
                                                           cAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCt                                                                                            >  2:329286/1‑150 (MQ=255)
                                                           cAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCt                                                                                            <  1:1119526/150‑1 (MQ=255)
                                                           cAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCt                                                                                            >  2:1338995/1‑150 (MQ=255)
                                                           cAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCt                                                                                            <  1:536337/150‑1 (MQ=255)
                                                           cAGCATATTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAAAATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCt                                                                                            <  2:1258300/150‑1 (MQ=255)
                                                             gCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTac                                                                                          <  1:329286/150‑1 (MQ=255)
                                                             gCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTac                                                                                          >  2:558989/1‑150 (MQ=255)
                                                                ttttATGATCGCCTGAATAATTGCCTGCCGTGGCCAACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACt                                                                                       >  1:449802/1‑150 (MQ=255)
                                                                    attaTCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATa                                                                                   <  1:620864/147‑1 (MQ=255)
                                                                         cGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACcagc                                                                              >  1:707944/1‑150 (MQ=255)
                                                                         cGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGAAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACcagc                                                                              >  2:636238/1‑150 (MQ=255)
                                                                             tGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGcc                                                                          <  2:701290/150‑1 (MQ=255)
                                                                             tGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGcc                                                                          <  2:1350522/150‑1 (MQ=255)
                                                                               aataatTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCtt                                                                        <  1:1338995/150‑1 (MQ=255)
                                                                                  aatTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAg                                                                     <  2:853678/150‑1 (MQ=255)
                                                                                   atTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGc                                                                    <  1:1121921/150‑1 (MQ=255)
                                                                                     tgcctgccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGa                                                                  <  1:312465/150‑1 (MQ=255)
                                                                                      gcctgccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGAt                                                                 >  2:1245281/1‑150 (MQ=255)
                                                                                      gcctgccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGAt                                                                 <  2:707944/150‑1 (MQ=255)
                                                                                      gcctgacGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGCAGTAACAAACCAGAAAAACCTGAGCCAAAATGATCAATTGCTACACTGATACAAGCAGCCTTTACCGAt                                                                 >  2:52916/1‑150 (MQ=255)
                                                                                           ccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAAt                                                            >  1:65721/1‑150 (MQ=255)
                                                                                             gTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCg                                                          >  2:296261/1‑150 (MQ=255)
                                                                                              tGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGc                                                         <  1:821933/150‑1 (MQ=255)
                                                                                               ggCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCt                                                        <  1:1245263/150‑1 (MQ=255)
                                                                                                 cccACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTc                                                      <  2:568449/150‑1 (MQ=255)
                                                                                                   cACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCgg                                                    >  1:141548/1‑150 (MQ=255)
                                                                                                    aCTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGc                                                   <  1:636238/150‑1 (MQ=255)
                                                                                                         taagAACATTGGTAATAAAATCCTAACTGATTATAATTCTGTAAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACt                                              >  2:1187760/4‑150 (MQ=255)
                                                                                                                   ggTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCg                                    <  2:895560/150‑1 (MQ=255)
                                                                                                                   ggTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCg                                    >  2:1380159/1‑150 (MQ=255)
                                                                                                                    gTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCga                                   >  1:146882/1‑150 (MQ=255)
                                                                                                                     tAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCgag                                  <  2:173866/150‑1 (MQ=255)
                                                                                                                          aaTTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATg                             <  1:96895/150‑1 (MQ=255)
                                                                                                                                    tgaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACttt                   >  1:1383389/1‑150 (MQ=255)
                                                                                                                                     gaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                  <  2:1422330/150‑1 (MQ=255)
                                                                                                                                     gaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                  <  1:1359668/150‑1 (MQ=255)
                                                                                                                                     gaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                  >  1:1109114/1‑150 (MQ=255)
                                                                                                                                      aTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCa                 >  1:962093/1‑150 (MQ=255)
                                                                                                                                       ttttAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAg                <  1:296261/150‑1 (MQ=255)
                                                                                                                                          tAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTaa             >  1:724075/1‑150 (MQ=255)
                                                                                                                                            ctccTGTCACGCTTCCAGCCTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAAt           >  1:656025/4‑150 (MQ=255)
                                                                                                                                            aTTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAAt           >  2:841138/1‑150 (MQ=255)
                                                                                                                                                 gTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAAc      <  1:655802/150‑1 (MQ=255)
                                                                                                                                                     ggCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCCCTTTCAGTAAATTCAACgatg  >  2:574009/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACATCCACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATG  >  NC_000913/1214362‑1214660

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: