Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,214,516 G→C 5.3% F320L (TTC→TTG) ‡ 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,5160GC 5.3% 323.7 / 6.9 95F320L (TTC→TTG) ‡bluFblue light‑responsive regulator of BluR
Reads supporting (aligned to +/- strand):  ref base G (41/49);  new base C (3/2);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-01

CACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTT  >  NC_000913/1214367‑1214664
                                                                                                                                                     |                                                                                                                                                    
cACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTgg                                                                                                                                                      >  1:533201/1‑150 (MQ=255)
 aCTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGa                                                                                                                                                     <  2:1476505/150‑1 (MQ=255)
  ctctTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGAACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGaa                                                                                                                                                    >  1:1628406/1‑150 (MQ=255)
    cttctGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGc                                                                                                                                                  >  1:121977/1‑150 (MQ=255)
       ctGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTTCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCATGCTGGAAGCGTg                                                                                                                                               <  2:1410506/150‑1 (MQ=255)
       ctGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTg                                                                                                                                               <  1:978084/150‑1 (MQ=255)
        tGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGa                                                                                                                                              >  2:1317552/1‑150 (MQ=255)
        tGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGa                                                                                                                                              <  1:1289245/150‑1 (MQ=255)
         ggTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGAc                                                                                                                                             <  2:1505925/150‑1 (MQ=255)
          gtgtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACa                                                                                                                                            <  2:1209800/150‑1 (MQ=255)
          gtgtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCACGCTTCCAGCCTGACa                                                                                                                                            <  2:977308/150‑1 (MQ=255)
           tgtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAg                                                                                                                                           <  2:1593814/150‑1 (MQ=255)
            gtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAgg                                                                                                                                          >  1:1865269/1‑150 (MQ=255)
            gtTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAgg                                                                                                                                          >  1:1872033/1‑150 (MQ=255)
              ttccACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCACGCTTCCAGCCTGACAgaat                                                                                                                                        <  1:227353/148‑4 (MQ=255)
              tGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAg                                                                                                                                        <  1:530004/150‑1 (MQ=255)
              tGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAg                                                                                                                                        <  1:784585/150‑1 (MQ=255)
              tGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAg                                                                                                                                        <  2:232003/150‑1 (MQ=255)
               gCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGt                                                                                                                                       >  2:614523/1‑150 (MQ=255)
                ccACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTa                                                                                                                                      >  1:1685357/1‑150 (MQ=255)
                  aCGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAAc                                                                                                                                    <  2:533201/150‑1 (MQ=255)
                     cccATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACaaa                                                                                                                                 >  1:452920/1‑150 (MQ=255)
                      ccATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAAc                                                                                                                                <  1:1173078/150‑1 (MQ=255)
                      ccATAGCAAAAACTTAAATATAAATTGATGTAAAGCATTTTAAGATCGCATGAATAATTGCCTGCCGTGGACCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAAc                                                                                                                                >  2:896208/1‑150 (MQ=255)
                       cATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAAcc                                                                                                                               <  1:1396329/150‑1 (MQ=255)
                       cATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAAcc                                                                                                                               <  2:460402/150‑1 (MQ=255)
                       cATACCACTCTCTTGAATTTCAAGTGATGTACAGCAGTTTATGATCGCCTGAATACTTGCGGGCCGTGGCCCACTTTTATGCACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAAcc                                                                                                                               <  2:344320/150‑1 (MQ=255)
                                gACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAAcc                                                                                                                      >  2:493752/1‑150 (MQ=255)
                                   ttGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTgc                                                                                                                   <  2:581175/150‑1 (MQ=255)
                                   ttGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTgc                                                                                                                   <  1:1199278/150‑1 (MQ=255)
                                    tGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTTGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTgcg                                                                                                                  >  1:1014029/1‑150 (MQ=255)
                                    tGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTgcg                                                                                                                  <  1:675070/150‑1 (MQ=255)
                                    tGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTgcg                                                                                                                  <  1:1216693/150‑1 (MQ=255)
                                    tGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTgcg                                                                                                                  >  2:1761545/1‑150 (MQ=255)
                                    tGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTgcg                                                                                                                  >  2:820664/1‑150 (MQ=255)
                                      aaTTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGcc                                                                                                                <  2:1872033/150‑1 (MQ=255)
                                      aaTTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGcc                                                                                                                <  2:1865269/150‑1 (MQ=255)
                                        tttCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCaa                                                                                                              >  1:30036/1‑150 (MQ=255)
                                        tttCAAGTGATGTACAGCATTTTATGATCGACTGAATAATTGACTGCCGTGCCACACTTTTATGAACATTGCTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGAAAAACCTGCGACaa                                                                                                              >  2:679539/1‑150 (MQ=255)
                                          tCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAGCATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCaaaa                                                                                                            >  1:543358/1‑150 (MQ=255)
                                           cAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAAt                                                                                                           <  2:1457170/150‑1 (MQ=255)
                                           cAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAAt                                                                                                           <  1:1341062/150‑1 (MQ=255)
                                               tGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATc                                                                                                       <  1:1012969/150‑1 (MQ=255)
                                                   ttACAGCATTTTATGATCGCCTTAATAATTTCGTGCCGTGGCCCACTTTTATGAAGATTGGTAATCATTTCCTGACTGATTTTAATTCTGTCATGCTGGATGAGTGACAGGAGTAACAAACCAGCAAAACCTGAGCCAAAATTATCAAtt                                                                                                   <  2:806774/149‑1 (MQ=255)
                                                       aGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTa                                                                                               <  2:4360/150‑1 (MQ=255)
                                                        gCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTac                                                                                              >  2:805561/1‑150 (MQ=255)
                                                        gCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTac                                                                                              <  1:183243/150‑1 (MQ=255)
                                                               aTGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATa                                                                                       >  1:1345930/1‑150 (MQ=255)
                                                                   tCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACcag                                                                                   <  1:1434553/150‑1 (MQ=255)
                                                                       cTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACcag                                                                                   >  2:819731/1‑146 (MQ=255)
                                                                       cTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACcag                                                                                   <  1:819731/146‑1 (MQ=255)
                                                                        tGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGcc                                                                              <  2:1120206/150‑1 (MQ=255)
                                                                         gAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCt                                                                             >  2:827545/1‑150 (MQ=255)
                                                                             aatTGCCTGCCGTGGCCCACTTTTATGAGCATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCACAACCTGCGCCAAAATGATCAATTGCTCCACTGATACCAGCAGCCTTTAg                                                                         >  2:1619281/1‑150 (MQ=255)
                                                                              atTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACcagca                                                                                 >  1:1967243/1‑141 (MQ=255)
                                                                              atTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACcagca                                                                                 <  2:1967243/141‑1 (MQ=255)
                                                                              atTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGc                                                                        <  1:1766317/150‑1 (MQ=255)
                                                                               ttGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCg                                                                       <  2:1611649/150‑1 (MQ=255)
                                                                                tgcctgccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGTAGCGTGAAAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGa                                                                      <  1:1482149/150‑1 (MQ=255)
                                                                                   ctgccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGAttt                                                                   >  2:535601/1‑150 (MQ=255)
                                                                                    tgccGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTa                                                                  >  2:1552432/1‑150 (MQ=255)
                                                                                             ccACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGAtt                                                                    <  2:1478100/139‑1 (MQ=255)
                                                                                             ccACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGAtt                                                                    >  1:1478100/1‑139 (MQ=255)
                                                                                             ccACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCg                                                         <  1:1019371/150‑1 (MQ=255)
                                                                                             ccACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCg                                                         <  1:308824/150‑1 (MQ=255)
                                                                                             ccACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTAGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCg                                                         >  2:195499/1‑150 (MQ=255)
                                                                                              cACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTAAATCGCtaaag                                                        <  1:448666/150‑5 (MQ=255)
                                                                                               aCTTTTATGAACATTGGTATTCATTTCCTGACTGATTTTAATTGTGTCAGGCTGGAAGCGTGACAGGATTAACATACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACTCTGTTAGCAGCTGCCTTTAGCGATTAAATCGCtaaagg                                                       <  1:8620/150‑6 (MQ=255)
                                                                                               aCTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTAAATCGCtaaagg                                                       <  1:493752/150‑6 (MQ=255)
                                                                                                  tttATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCaaa                                                    >  2:1429010/1‑150 (MQ=255)
                                                                                                     aTGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACtc                                                 <  2:366288/150‑1 (MQ=255)
                                                                                                      tGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACtca                                                >  2:974225/1‑150 (MQ=255)
                                                                                                           aTTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCaaa                                           >  1:1618364/1‑150 (MQ=255)
                                                                                                                 aaTCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCgaga                                     >  2:52179/1‑150 (MQ=255)
                                                                                                                    cAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGatat                                  >  2:1925141/1‑150 (MQ=255)
                                                                                                                       ttCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGAc                               >  2:270749/1‑150 (MQ=255)
                                                                                                                          ctgactgaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTc                            <  1:1925141/150‑1 (MQ=255)
                                                                                                                              ctgaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACtt                        >  1:837638/1‑150 (MQ=255)
                                                                                                                                gaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                      <  1:1761545/150‑1 (MQ=255)
                                                                                                                                gaTTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTc                      <  2:1457136/150‑1 (MQ=255)
                                                                                                                                  ttttAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAg                    <  2:30036/150‑1 (MQ=255)
                                                                                                                                   tttAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGt                   >  2:1303426/1‑150 (MQ=255)
                                                                                                                                     tAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTaa                 >  1:1153503/1‑150 (MQ=255)
                                                                                                                                      actccTGTCACGCTTCCAGCCTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTaaa                >  1:1699582/5‑150 (MQ=255)
                                                                                                                                      actccTGTCACGCTTCCAGCCTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTaaa                >  1:627442/5‑150 (MQ=255)
                                                                                                                                        tccTGTCACGCTTCCAGCCTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAAtt              >  2:1627337/3‑150 (MQ=255)
                                                                                                                                          cTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCa            >  2:1969896/1‑150 (MQ=255)
                                                                                                                                          cTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCa            >  1:1990932/1‑150 (MQ=255)
                                                                                                                                           tGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCaa           >  2:1685110/1‑150 (MQ=255)
                                                                                                                                           tGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCaa           <  1:805561/150‑1 (MQ=255)
                                                                                                                                            gTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAAc          <  2:1685357/150‑1 (MQ=255)
                                                                                                                                             tCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACg         <  1:614523/150‑1 (MQ=255)
                                                                                                                                              cAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGCTATGACTTCACTTTCAGTAAATTCAACga        >  2:627669/1‑150 (MQ=255)
                                                                                                                                              cAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACga        >  1:1263213/1‑150 (MQ=255)
                                                                                                                                                ggCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACgatg      >  1:1205050/1‑150 (MQ=255)
                                                                                                                                                   tGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGAtt   <  2:1628406/150‑1 (MQ=255)
                                                                                                                                                    ggAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGAttt  <  2:543358/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CACTCTTCTGGTGTTGCCACGCCCATAGCACTGACTTGAATTTCAAGTGATGTACAGCATTTTATGATCGCCTGAATAATTGCCTGCCGTGGCCCACTTTTATGAACATTGGTAATCAATTCCTGACTGATTTTAATTCTGTCAGGCTGGAAGCGTGACAGGAGTAACAAACCAGCAAAACCTGCGCCAAAATGATCAATTGCTACACTGATACCAGCAGCCTTTAGCGATTTAATCGCTTCGGCAAACTCATCAAACCGAGATATGACTTCACTTTCAGTAAATTCAACGATGATTT  >  NC_000913/1214367‑1214664

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: